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^i^This  paper  introduces  a new  approach  for  combining  the  integral  equation  and  high 
frequency  asymptotic  techniques,  e.g.,  the  geometrical  theory  of  diffraction.  The 
method  takes  advantage  of  th^  fact  that  the  Fourier  transform  of  the  unknown  surface 
current  distribution  Is  proportional  to  the  scattered  far  field.  A number  'of  asymp- 
totic methods  are  currently  available  that  provide  good  approximation  to  this  far  field 
In  a convenient  analytic  form  which  Is  useful  for  deriving  an  initial  estimate  of  the 
Fourier  transform  of  the  current  distribution. 

An  Iterative  scheme  Is  developed  for  systematically  Improving  the  Initial  form  of 
the  high  frequency  asymptotic  solution  by  manipulating  the  Integral  equation  In  the 
Fourier  transform  domain. 

A synthetlc-aperture-dlstrlbutlon  scheme  is  also  developed  In  which  the  approxi- 
mate scattered  far-fleld  pattern  obtained  by  asymptotic  techniques  is  Improved  by  sys- 
tematically correcting  the  scattered  field  distribution  on  an  aperture  erected  in  Jux- 
taposition with  the  obstacle.  The  Introduction  of  such  a planar  aperture  not  only  pro- 
vides an  additional  degree  of  freedom  In  performing  Improving  operations,  but  also 
renders  the  scheme  to  handle  n-dlmenslonal  geometries  by  (n  - 1) -dimensional  fast 
Fourier  transform  (FFT),  where  n - 2,  3,  and  circumvents  the  unwieldy  three-dimensional 
FFT,  making  It  a conceptually  simple  and  computationally  efficient  method.  J 

A salient  feature  of  the  method  Is  that  It  provides  convenient  validity  checks  of 
the  solutions  for  the  surface  current  distribution  and  the  scattered  far-fleld  pattern 
by  verifying  that  the  scattered  field  obtained  Indeed  satisfies  the  boundary  conditions 
at  the  surface  of  the  scatterer.  Another  Important  feature  of  the  method  Is  that  it 
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13. 

Diffractions  by  a strip,  a thin  plate,  a rectangular  cylinder,  and  a circular 
cylinder  are  presented  as  Illustrative  examples  that  demonstrate  the  usefulness  of  the 
approach  for  handling  a variety  of  electromagnetic  scattering  problems  in  the  resonance 
region  and  above.  Some  concluding  remarks  and  comparison  with  other  methods  are  also 
Included. 
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PREFACE 

A little  over  a century  ago»  Maxwell  put  together  the  fundamental 
laws  that  govern  all  the  electromagnetic  phenomei.  In  an  amazingly 
elegant  form,  l.e.,  the  set  of  equations  named  after  him.  Ever  since 
then,  the  major  contributions  in  the  field  of  electromagnetics  have  been 
solving  boundary  value  problems,  or  seeking  solutions  which  satisfy 
the  Maxwell's  equations  and  match  the  boundary  conditions  (including 
radiating  conditions  as  special  cases) in  specific  environments  under 
consideration. 

, I 

The  dyadic  Green  s functions  technique,  an  elegant  and  powerful 
method  for  solving  boundary-value  problems,  was  first  formulated  by 
Schwinger  In  the  early  1940' s.  However,  solutions  obtained  by  that 
technique  are  extremely  complicated  and  are  not  always  best  suited  for 
numerical  computations. 

In  the  late  1950' s,  Keller  conceived  the  concept  of  a diffraction-  I 

coefficient  approach  to  the  high-frequency  scattering,  and  was  able 

c 

1 

to  obtain  approximate  solutions  for  far-fleld  patterns  in  a very  simple 

and  straightforward  manner.  Although  his  method  falls  at  certain  aspect  i 

j 

angles  In  space.  It  has  been  applied  to  solve  many  practical  problems.  | 

In  the  1960's,  the  moment  method  became  popular  due  to  the  avail-  I 

j 

ability  of  the  large  scale  computer  systems.  However,  for  electrically  : 

large  scatterers,  the  moment  method  Is  limited  by  the  storage  of  a 

i 

computer. 

In  the  early  1970's,  a trend  of  combining  the  Keller's  geometrical 
theory  of  diffraction  with  the  moment  method  to  solve  high-frequency  i 

scattering  problems  was  Initiated.  However,  a general  approach 

for  combining  the  two  techniques  has  yet  to  be  developed.  | 


viil 


It  is  the  purpose  of  this  study  to  Introduce  a method  of  combining 
the  Integral  equation  and  asymptotic  techniques  for  solving  electro- 
magnetic scattering  problems.  The  method  Is  conceptually  simple  and 
computationally  efficient.  Comparisons  with  contemporary  combinational 
approaches  show  that  the  method,  when  fully  developed,  appears  to  be 
highly  promising  In  solving  practical  problems.  Of  course,  this  thesis 
Is  merely  a start,  and  there  Is  much  left  to  be  studied  and  investigated. 
It  Is  the  author's  hope  that  future  research  activities  along  the  same 
lines  as  described  herein  will  result  In  not  only  a feasible  but  also 
a practical  way  of  handling  real-world  high-frequency  scattering  problems. 
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ABSTRACT 


This  paper  Introduces  a new  approach  for  combining  the  Integral 
equation  and  high  frequency  asymptotic  techniques,  e.g.,  the  geometrical 
theory  of  diffraction.  The  method  takes  advantage  of  the  fact  that  the 
Fourier  transform  of  the  unknown  surface  current  distribution  is  pro- 
portional to  the  scattered  far  field.  A number  of  asymptotic  methods 
are  currently  available  that  provide  good  approximation  to  this  far 
field  in  a convenient  analytic  form  which  is  useful  for  deriving  an 
initial  estimate  of  the  Fourier  transform  of  the  current  distribution. 

An  iterative  scheme  is  developed  for  systematically  improving  the 
initial  form  of  the  high  frequency  asymptotic  solution  by  manipulating 
the  integral  equation  in  the  Fourier  transform  domain. 

A synthetlc-aperture-dlstrlbutlon  scheme  is  also  developed  in  which 
the  approximate  scattered  far-fleld  pattern  obtained  by  asymptotic  tech- 
niques is  improved  by  systematically  correcting  the  scattered  field 
distribution  on  an  aperture  erected  in  juxtaposition  with  the  obstacle. 
The  introduction  of  such  a planar  aperture  not  only  provides  an  additional 
degree  of  freedom  in  performing  improving  operations,  but  also  renders 
the  scheme  to  handle  n-dlmenslonal  geometries  by  (n  - 1) -dimensional 
fast  Fourier  transform  (FFT),  where  n - 2,3,  and  circumvents  the  un- 
wieldy three-dimensional  FFT,  making  it  a conceptually  simple  and  com- 
putationally efficient  method. 


A salient  feature  of  the  method  is  that  It  provides  convenient 
validity  checks  of  the  solutions  for  the  surface  current  distribution 
and  the  scattered  far-fleld  pattern  by  verifying  that  the  scattered 
field  obtained  Indeed  satisfies  the  boundary  conditions  at  the  surface 
of  the  scatterer.  Another  Important  feature  of  the  method  Is  that  It 
yields  both  the  Induced  surface  current  density  and  the  far  field. 

Diffractions  by  a strip,  a thin  plate,  a rectangular  cylinder,  and 
a circular  cylinder  are  presented  as  Illustrative  examples  that  demon- 
strate the  usefulness  of  the  approach  for  handling  a variety  of  electro- 
magnetic scattering  problems  In  the  resonance  region  and  above.  Some 
concluding  remarks  and  comparison  with  other  methods  are  also  Included. 
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1.  INTRODUCTION 

It  is  well-known  that  the  Integral  equation  methods  are  limited 
In  application  to  scatterers  whose  characteristic  dimensions  are  of 
the  order  of  one  wavelength  or  less.  On  the  other  hand,  the  high- 
frequency  asymptotic  techniques  can  be  reliably  used  only  when  the 
scatterer  Is  large  compared  to  the  wavelength,  and  neither  of  the  above 
two  methods  Is  suitable  In  the  resonance  region. 

This  thesis  Introduces  a new  hybrid  technique,^  based  on  a com- 
bination of  the  Integral  equation  and  asymptotic  methods,  that  Is  use- 
ful In  the  entire  frequency  range  encompassing  the  resonance  region 
and  above.  Another  Important  feature  of  the  method  Is  that  It  can  be 
used  to  check  and  Improve  the  accuracy  of  high-frequency  asymptotic 
solutions.  Such  an  accuracy  test  and  systematic  Improvement  of  the 
asymptotic  solution  are  often  needed,  but  no  reliable  methods  for  per- 
forming these  are  available  at  the  present  time. 

In  contrast  to  the  ray  optics  methods,  which  are  based  on  the  dif- 
fraction of  ray  fields  as  determined  by  the  local  properties  of  the  sur- 
face of  the  scatterer,  the  present  method  starts  with  the  Integral 
equation  formulation  Incorporating  the  boundary  conditions  on  the  entire 
surface  of  the  scatterer.  Conventionally,  th'*  solution  of  the  Integral 
equation  for  the  Induced  surface  current  Is  carried  out  by  matrix  methods 
[l  ]>  [2  ].  The  size  of  the  matrix  becomes  prohibitively  large  and  Its 
solution  extremely  time-consuming  when  the  characteristic  dimension  of 
the  scatterer  approaches  the  wavelength  of  the  Illuminating  field. 

^The  original  concepts  on  which  this  thesis  Is  based  were  described  at  the 
1975  URSI  Symposium  In  a paper  entitled,  "A  new  method  for  Improving  the 
GTD  solution  via  the  Integral  equation  formulation,"  by  R.  Mlttra  and 
W.  L.  Ko. 


The  approach  developed  In  this  thesis  circumvents  this  difficulty 
while  still  preserving  the  self-consistent  nature  of  the  Integral  equa- 
tion formulation  by  constructing  the  solution  In  the  Fourier  transform 
or  spectral  domain  rather  than  In  the  space  domain.  We  take  advantage 
of  the  facts  that  the  Fourier  transform  of  the  surface  current  distribu- 
tion is  directly  proportional  to  the  far  scattered  field  and  that  the 

asymptotic  methods  often  provide  a very  good  Initial  estimate  of  the 

2 

latter  quantity.  We  next  construct  an  iterative  solution  of  the  In- 

3 

tegral  equation  In  the  transform  domain  with  the  GTD  or  other  high- 
frequency  solution  as  the  zero-order  approximation. 

As  shown  In  the  thesis,  this  procedure  not  only  allows  us  to  Im- 
prove on  the  GTD  or  similar  solutions  but  also  provides  a convenient 
means  for  testing  the  satisfaction  of  the  boundary  conditions  on  the 
surface  of  the  scatterer.  Furthermore,  the  method  yields  not  only  the 
far  field  but  also  the  Induced  surface-current  distribution,  a feature 
not  readily  available  In  some  other  high-frequency  techniques. 

A detailed  derivation  of  the  Iterative  method  Is  presented  In 
Chapter  2.  A step-by-step  recipe  for  applying  the  Iterative  scheme 
to  construct  the  solution  of  the  Induced  surface  current  can  also  be 
found  In  Chapter  2. 

The  application  of  the  general  procedure  outlined  In  Chapter  2 Is 
Illustrated  by  two  examples:  the  two-dimensional  problem  of  a plane 

2 

A moment  method  solution  In  the  spectral  domain  has  also  been  devel- 
oped (see  [3  ]). 

3 

For  a comprehensive  discussion  of  the  Geometrical  Theory  of  Diffraction 
(GTD)  and  Its  application,  see  R.  G.  Kouyoumjlan,  Chapter  7 of  Numerical 
and  Asymptotic  Techniques  In  Electromagnetics.  R.  Mlttra  (Ed.), 

Springer  Verlag,  1975. 
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wave  diffraction  by  a strip  and  a three-dimensional  problem  of  a plane 
wave  diffraction  by  a thin  plate. 

In  Chapter  3,  the  iteration  procedure  is  applied  to  solve  the  strip 
problem.  This  problem  was  chosen  for  the  following  reasons:  It  is  shown 

that  when  the  angle  of  incidence  is  normal  or  near  normal,  the  GTD  solu- 
tion accurately  satisfies  the  boundary  condition  > 0 on  the  strip 
even  when  the  multiple  interaction  between  the  two  edges  of  the  strip 
is  neglected.  However,  it  is  found  that  when  the  angle  of  incidence  is 
near  grazing,  the  GTD  solution  is  quite  unsatisfactory,  while  the  it- 
erated solution  generated  by  the  hybrid  technique  does  display  the  cor- 
rect behavior. 

In  Chapter  4,  the  Iteration  procedure  is  applied  to  solve  the  plate 
problem.  This  problem  was  chosen  for  the  following  reasons:  In  the 

plate  problem,  the  difficulty  in  applying  GTD  to  this  geometry  stems 
from  the  fact  that  the  diffraction  coefficient  for  the  corners  of  the 
plate  is  not  known  and  neglecting  the  corner  effects  can  cause  substan- 
tial errors  in  the  resonance  frequency  region  where  the  plate  size  is  of 
the  order  of  one  wavelength  squared.  However,  the  present  technique, 
based  on  a combination  of  the  integral  equation  and  asymptotic  methods, 
does  allow  the  successful  computation  of  the  Induced  surface  currents 
on  the  plate. 

In  Chapter  5,  diffraction  by  planar  surface  solid  scatterers  with 
multiple  edges,  e.g.,  a rectangular  cylinder, is  considered.  A technique 
which  is  based  on  a representation  of  the  scattered  fields  in  terms  of 
the  spectrum  of  the  Induced  surface  current  on  the  scatterer  rather 
than  the  rays  emanating  from  it  is  discussed.  The  scattered  far-fleld 
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obtained  by  this  technique  compares  extremely  well  with  that  obtained 
from  the  moment  method  solution  for  thirty-two  unknowns.  Comparisons 
with  results  obtained  by  other  approaches  In  the  literature  [4  ] can 
also  be  found  In  Chapter  5. 

In  Chapter  6,  the  method  for  combining  Integral  equation  and  asymp- 
totic techniques  for  solving  electromagnetic  scattering  problems  of  a 
convexly  curved  surface  scatterer,  e.g.,  a circular  cylinder* Is  devel- 
oped. Detailed  descriptions  of  the  synthetlc-aperture-dlstrlbutlon 
approach  are  documented,  and  comparisons  with  other  methods  [4  ],  [5  ], 
that  have  been  developed  recently  for  solving  problems  of  a similar 
nature  are  also  Included  In  Chapter  6.  The  Idea  Is  to  transfer  Informa- 
tion characterizing  the  curved  surface  onto  a planar  aperture  where 
operations  for  Improving  the  high-frequency  asymptotic  solution  are 
carried  out.  A salient  feature  of  the  method  Is  that  the  far  field  Is 
obtained  by  Fourier  transforming  the  aperture-field  distribution;  hence, 
It  can  handle  problems  of  an  n-dlmenslonal  obstacle  by  an  (n  - 1)- 
dlmenslonal  fast  Fourier  transform  (FFT),  where  n ■ 2,3.  Therefore,  the 
unwieldy  three-dimensional  FFT  Is  circumvented;  hence, the  method  Is 
numerically  efficient.  Moreover,  the  method  Is  conceptually  simpler 
than  other  conventional  methods  because  operations  for  Improvement  are 
carried  out  on  a plane  rather  than  on  a curved  surface. 

In  Chapter  7,  concluding  remarks  are  made  and  future  research 
activities  along  similar  lines  based  on  the  present  approach  are  recom- 
mended. 

In  this  method,  the  fast  Fourier  transform  [6  ] Is  a necessary 
tool;  hence,  a brief  discussion  on  FFT  Is  presented  In  Appendix  A. 
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2.  FORMULATION  OF  AM  ITERATION  METHOD 
2.1  Derivation  of  the  Method 

We  begin  our  analysis  with  the  electric-field  Integral  equation 
[7  ] for  a perfectly  conducting  scatterer.  The  equation  may  be  sym- 
bolically written  as 

(G  * j)^  » -eJ  , r e S (2.1) 

where  J(r')  Is  the  Induced  surface  current  density  we  are  attempting 
to  determine,  and  the  subscript  t signifies  tangential  to  the  surface  S; 

Is  the  tangential  component  of  the  Incident  electric  field  E^  on  the 
surface  S of  the  scatterer;  and  G is  the  free  space  Green's  Dyadic,  which 
yields  the  scattered  electric  field  when  operating  on  J.  In  anticipation  of 
Fourier  transforming  (2.1),  we  extend  it  over  all  space  by  first  defining 
a truncation-projection  operator  0 as  follows:  i 

For  any  vector  A(r),  which  Is  a vector  function  of  position  r,  ' 

6(A)*/a  6(r  - r )dr  , r € S (2.2a)  ■ 

U S S ^ 

t 

where  5 (r  - r ) Is  the  Dirac  delta  function  and  the  subscript  t sig- 
8 

nlfles  tangential  to  the  surface  S.  Also  Its  complementary  operator 
6 Is  defined  as 

0(A)  - A - 0(5) 

We  can  then  rewrite  (2.1)  as 

G * j « 0(-E^)  + e(G  * (0J))  , for  all  space 


t 


t 


6 

As  Indicated  above,  (2.3)  is  valid  at  all  observation  points  whether 
on  or  off  the  surface  S.  Note  that  the  Integral  equation  (2.1)  is  em- 
bedded in  (2.3)  and  that  we  have  made  use  of  the  obvious  identity  6J  ■ J. 

Next  we  Fourier  transform  (2.3)  by  introducing  the  transform  re- 
lationships 

OP  • • 

F(k)  - / F(?)  d^  - F(F(?)1  (2.4a) 

«Q0 

and 


00  — 

F(?)  - (-2^)^  / F(k)  dk  - F"^(F(k)]  (2.4b) 

>00 

with  *■  on  top  denoting  the  transformed  quantities. 

The  transformed  version  of  (2.3)  reads 

GJ  - -Ej,  + F (2.5) 

where  F - F[6(G  * (6J))]  and  is  the  transform  of  the  tangential  com- 
ponent of  the  incident  field  truncated  on  S.  Note  that  the  convolution 
operation  in  (2.3)  is  transformed  into  an  algebraic  product  upon  Fourier 
transformation. 

A formal  solution  to  (2.5)  can  now  be  written 

j - g“^(-Ej  + F)  . (2.6) 

Equation  (2.6)  says  that  if  we  had  available  the  Fourier  transform  of 
the  scattered  electric  field,  we  could  construct  the  solution  for  the 
Induced  surface  current  density  in  the  transform  domain  by  adding  it 
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to  -E^,  which  is  known,  and  by  performing  an  algebraic  division  repre- 
sented by  G~^.  In  practice,  of  course,  F is  not  known  and  must  be  solved 
for  along  with  J if  (2.6)  is  to  be  used  in  the  form  as  shown.  However, 
Instead  of  using  this  form,  we  proceed  to  derive  an  iterated  form  of 
the  equation  as  shown  below: 


■j(n+l)  _ g"^(-Ej  + 


(2.7) 


which  indicates  the  (n  + l)th  approximation  of  J from  the  nth  approxi- 
mation for  F.  We  next  show  how  F^”^  Itself  can  be  derived  from 
To  this  end,  we  use  the  identity 

f - F[F"^[g5]  - 6(F‘^[g!])]  (2.8) 


which  may  be  verified  by  writing  (2.8)  as 
F - F[G  * J - e(-E^)] 


(2.9) 


and  using  (2.3)  to  get 
F - F[0(G  * (0J))] 


(2.10) 


which,  of  course,  is  the  definition  of  F.  We  can  now  use  (2.8)  to  de- 
rive the  nth  approximation  F^"^  of  F from  the  nth  approximation  of  J, 
l.e.,  3^*^^  The  relationship  is  written  as 


- F[F"^[g5^"^]  - 0(F“^[5j^"^])]  . (2.11) 

The  desired  iteration  relating  and  may  now  be  written. 

Using  (2.7)  and  (2.11) 


5(nfl) 


+ F[F‘^[g5^"^]  - 0(F”^[gJ^"^])]] 


(2.12) 
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2.2  Recipe  for  Applying  the  Method 

The  step-by-step  procedure  for  constructing  the  solution  of  the 
transformed  surface  current  J will  now  be  given: 

1.  Begin  with  an  estimate  of  which  Is  the  Fourier  trans- 

form of  the  Induced  surface  current,  or  equivalently,  the 
scattered  far  field  within  a known  multiplicative  constant. 

Typically,  the  Initial  approximation  for  J,  viz.,  J can  be  ob- 

tained as  follows: 

(a)  Estimate  F,  the  Fourier  transform  of  the  scattered  field,  F, 

4 

outside  the  scatterer,  using  GTD  or  other  asymptotic  solutions. 

(b)  Subtract  E^,  the  Fourier  transform  of  the  tangential  component 
of  the  Incident  electric  field  truncated  to  the  surface  of  the 
scatterer. 

(c)  Multiply  the  result  of  Step  (b)  by  5”^.  Note  that  G ^ is 
known  and  the  operation  Is  algebraic. 

(d)  Take  the  inverse  Fourier  transform  of  the  result  of  Step  (c) 
and  truncate  It  to  the  surface  of  the  scatterer  to  obtain 

the  initial  approximation  for  J. 

2.  Multiply  by  G,  the  known  transform  of  the  Green's  Dyadic. 

Note  this  Involves  algebraic  multiplication  and  not  the  usual 
time-consuming  convolution  operation. 


^Note  that  GTD  (Keller's)  solutions  may  either  have  singularities  or  may 
be  In  error  near  shadow  and  reflection  boundaries  or  at  caustics,  and  the 
Uniform  Theory  of  Diffraction  (UTD)  [8  ] and  the  Uniform  Asymptotic  Theory 
(UAT)  [9  ] break  down  at  caustics.  The  Spectral  Theory  of  Diffraction  (STD) 
[10],  on  the  other  hand,  is  uniform  for  all  observation  angles.  The  cri- 
terion for  choosing  any  of  these  asymptotic  forms  of  solution  is  conve- 
nience of  computation  for  desired  accuracy.  For  a comparative  evaluation 
of  the  accuracy  of  the  GTD,  UTD,  UAT  and  STD,  the  reader  is  referred  to  [11]. 
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3.  Take  the  inverse  Fourier  transform  of  the  product  using 

both  visible  and  invisible  ranges. 

4.  Apply  the  truncation-projection  operator  9 to  F which 

gives  the  approximation  to  the  tangential  component  of  the 
scattered  electric  field  on  the  surface  S.  The  accuracy  of 
the  solution  can  be  conveniently  checked  at  this  point  by  verl- 
fying  the  satisfaction  of  the  boundary  condition  by  the  tangential 
component  of  E®,  viz.,  {E®  ■ mentioned  in  Chapter  1, 

this  is  an  Important  feature  of  the  method. 

5.  Subtract  0(F  from  the  total  F already  eval- 

uated. 

6.  Take  the  Fourier  transform  of  the  difference  obtained  In  Step  5. 

7.  Subtract  E^,  the  Fourier  transform  of  the  tangential  component 
of  the  incident  electric  field  truncated  on  the  surface,  from 
the  result  in  Step  6. 

8.  Multiply  the  result  obtained  in  Step  7 by  G Note  that  G ^ 

is  also  known  and  the  operation  is  again  algebraic  as  in  Step  2. 

The  result  thus  obtained  is  which  is  the  first  iteration 

of  the  scattered  far  field. 

9.  Take  the  Inverse  Fourier  transform  of  obtained  in  Step  8 
and  evaluate  it  on  S to  get  the  desired  Induced  surface  current 

on  the  scatterer.  In  other  words,  perform  the  operation  I 

; 

6(F  For  an  exact  solution,  this  operation  is  redundant,  \ 

{ 

since  J ■ 0J,  and  hence,  0(F  ^[F[0J]])  ■ 00J  ■ J.  However,  the  I 

i 

Fourier  inversion  of  an  nth  approximate  solution  will  not  I 

give  rise  to  a current  distribution  that  is  nonzero  except  on  S.  ] 

i 

This  step  provides  a test  for  the  accuracy  and  for  the  'i 

\ 

I 

j 
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convergence  of  the  approximate  solution  by  comparing  the  approx- 
imate 5^°^  with  F[e(F"^[J^^^l)]. 

10.  Take  F[0(F  to  derive  an  improved  approximation  for 

-fi) 

11.  Repeat  as  necessary  using,  for  instance,  the  improved  ^ 

from  Step  10  in  the  Iteration  Equation  (2.12)  to  generate  the 

-(2) 

next  higher-order  approximation  J . 

Before  closing  this  chapter  we  should  point  out  that  Galerkln's 
method  applied  in  the  transformed  domain  also  provides  an  alternate,  and 
in  some  cases  the  more  desirable,  approach  for  deriving  the  solution  to 
the  transformed  Integral  equation.  etc.,  as  well  as  other  suit- 

able functions  may  be  employed  as  a basis  set  for  this  purpose. 

In  the  following  two  chapters,  we  show,  in  some  detail,  the  appli- 
cation of  the  procedure  just  described  to  a two-dimensional  and  a three- 
dimensional  scattering  problem. 
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3.  DIFFRACTION  BY  AN  INFINITE  STRIP 


3.1  Introduction 

In  the  last  chapter,  we  presented  a general  iteration  method  for  ob- 
taining solution  of  the  Integral  equation  In  the  transform  domain  with  the 
GTD  or  other  high-frequency  solution  as  the  zero-order  approximation.  This 
Iteration  method  not  only  allows  us  to  Improve  on  the  GTO  or  similar  solu- 
tions but  also  provides  a convenient  means  for  testing  the  satisfaction 
of  the  boundary  conditions  on  the  surface  of  the  scatterer.  Furthermore, 
the  method  yields  not  only  the  far-fleld  but  also  the  Induced  surface- 
current  distribution,  a feature  not  readily  available  In  some  other 
high-frequency  techniques. 

The  application  of  the  general  procedure  outlined  In  Chapter  2 Is 
Illustrated  In  this  chapter  by  using  it  to  solve  the  two-dimensional 
problem  of  a plane  wave  diffraction  by  a strip.  This  problem  was  chosen 
for  the  following  reasons:  It  Is  shown  that  when  the  angle  of  Incidence 

Is  normal  or  near  normal,  the  GTD  solution  accurately  satisfies  the 
boundary  condition  ■ 0 on  the  strip  even  when  the  multiple  inter- 

action between  the  two  edgei^  of  the  strip  Is  neglected.  However,  It 
Is  found  that  when  the  angle  of  Incidence  Is  near  grazing,  the  GTD  solu- 
tion Is  quite  unsatisfactory,  while  the  Iterated  solution  generated 
by  the  hybrid  technique  does  display  the  correct  behavior. 

3.2  Geometry  of  the  Strip  Problem 

The  geometry  of  the  electromagnetic  scattering  problem  Involving 
a perfectly  conducting  Infinite  strip  of  zero  thickness  Illuminated 
by  a uniform  plane  wave,  whose  electric  Intensity  vector  Is  oriented 
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parallel  to  the  edges  of  the  strip,  Is  depicted  in  Figure  3.1.  For 
convenience  of  analysis,  an  arbitrary  incident  wave  can  always  be  decom- 
posed Into  two  components  with  respect  to  the  z-axls,  namely,  TM 

Z 

(E-wave)  and  TE  (H-wave).  In  the  following  discussion  we  consider  the 
z 

E-wave  case  only;  the  H-wave  case  can  be  solved  In  a similar  manner 
by  considering  - zH^. 

The  Incident  field  Is  given  by 

. -lk(xcos<|i«+y8i*'<l'n) 

eJ(P,«  - e 0 0 (3.1) 

z 

where  the  e time  dependence  Is  understood. 

3.3  Iteration  Method  Applied  to  the  Strip  Problem 

The  integral  equation  formulation  [12]  for  the  problem  at  hand 
takes  the  form 

j a 

-E^(x)  - / J (x’)G(x  - x')  dx*  , X e [-a, a]  (3.2) 

z * z 

-a 

where  J (x')  Is  the  algebraic  sum  of  the  Induced  surface  current  den- 
z 

sltles  on  the  top  and  the  bottom  surfaces  of  the  thin  strip.  The  kernel 
G is  the  two-dimensional  free-space  Green's  function  given  by 

G(x  - x')  - HQ^^(kQ|x  - x'|)  (3.3) 

where  Is  the  Hankel's  function  of  the  first  kind  of  order  zero, 

kg  - 2ir/X  Is  the  free-space  propagation  constant.  Note  that  (3.2)  is 
the  conventional  Integral  equation  which  equates  the  Integral  repre- 
sentation of  the  tangential  component  of  the  scattered  E-fleld  radiated 


13 


Figure  3.1.  Diffraction  by  a strip  Illuminated  by  an  E-wave. 


HALF  -e'  half 
PLANE  PUNE  F, 

« Z 


a a 

STRIP 


(a)  (b) 

Figure  3.2.  F.(x)  can  be  approximated  by  the  GTD  solution  to  the  half- 
plane  problem  (a)  shown  on  the  left-hand  side;  F2(x)  can 
be  approximated  by  the  GTD  solution  to  the  half-plane 
problem  (b)  shown  on  the  right-hand  side. 
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by  the  Induced  surface  current  density  to  the  negative  of  the  tangential 
component  of  the  Incident  E-fleld  on  the  surface  of  the  perfectly  con- 
ducting scatterer  as  required  by  the  satisfaction  of  the  boundary  con- 
dition. Hence,  (3.2)  Is  valid  on  the  strip  only. 

An  extended  Integral  equation  that  Is  valid  for  all  x can  be  ob- 
tained by  Including  the  scattered  fields  outside  the  strip  as  well. 

If  the  scattered  field  on  the  Interval  (-»,-a)  Is  designated  by  F^(x) 
and  the  scattered  field  on  the  Interval  (a,<»)  Is  designated  by  F2(x), 
then  the  extended  form  of  (3.2)  becomes  [13] 

a . 

/ J (x')G(x  - X*)  dx'  - 0(-E  (x))  + F,  (x)  + F,(x)  (3.4) 

‘ Z Z 1 z 

-a 


where  6 Is  defined  Iti  (2.2a). 

Since  the  Fourier  transform  of  the  Induced  surface  current  density 
can  be  related  to  the  far  field,  (3.4)  Is  Fourier  transformed  to  give 


J^(o)G(o) 


e(-Ep(o)  + Fj^(a)  + F2(a) 


(3.5) 


where  ~ on  top  Indicates  the  Fourier  transform  pair  defined  In  (2.4) 
which  simplifies  In  the  present  one-dlmenslonal  problem  to 


00 

F(o)  - / F(x)  e"^“*  dx 

•>00 


(3.6a) 


and 

00 

F(x)  " •^  / F(o)  da 

—00 


(3.6b) 
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The  Fourier  transform  of  the  two-dimensional  Green's  function  In 
(3.5)  takes  the  form 


G(a) 


(3.7) 


Note  that  (3.5)  Is  an  algebraic  equation  In  the  spectral  domain  In  con- 
trast to  the  convolution  form  of  the  Integral  equation  (3.4)  In  the 
spatial  domain.  The  reason  for  working  In  the  spectral  domain  will  be- 
come clear  when  the  method  of  solution  for  (3.5)  Is  developed.  Following 
the  procedure  discussed  In  Chapter  2 and  In  terms  of  the  notations  Intro- 
duced In  the  present  problem,  we  proceed  as  follows: 

1.  Obtain  J^^^(a),  the  Initial  approximation  of  the  Fourier  transform 
of  the  Induced  surface  current  density,  or  equivalently,  the  scattered  far 
field  within  a known  multiplicative  constant,  as  follows: 

1.1.  Find  the  expressions  for  the  first  estimate  of  F^^^(a)  + F^^^a). 
Note  that  GTD  may  be  used  to  get  closed-form  expressions  for 
F^^^(a)  and  F^^^  (o)  since  Fj^(x)  and  F2(x)  can  be  approximated 
by  the  GTD  solutions  to  the  two  half-plane  problems  as  shown 
In  Figure  3.2  (page  13).  The  expressions  for  F^®^  and  F^^^ 
as  obtained  from  GTD  read 


(a) 


" lkacos(fiQ 
2 ® 


/sit  sin 


Q 

2 J (a  + 


laa 

e 

k cos  (|)Q)*''a  + k 


la(a+kcos(|)Q) 

le 

(o  + k cos  ^q) 


(3.8a) 


and 


I'.'l  . 'fv 
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(o) 


^ " -ikacos(>Q 
2 ® 


-laa 

e 

(a  + k cos  ♦q)*'V'  - 


a 


le 


-iaCotfkcosi^Q) 


(a  + k cos  ^q) 


(3.8b) 


Note  that  these  expressions  are  free  of  singularities  for  all  o. 

1.2.  Solve  for  the  initial  approximation  of  J_(a) , (a) , by 

z z 

carrying  out  the  operations  shown  below; 

^0(-eS(c»)  + F^°^(a)  + 


z 


/ 

- F 

6 

r 

G(a) 


2.  Use  (2.12)  to  further  Improve  the  solution  as  necessary. 

The  check  for  satisfaction  of  the  integral  equation  can  be  applied 
very  simply  by  computing  J(a)G(a),  taking  its  inverse  Fourier  transform, 
and  verifying  how  well  it  approaches  -E^  on  the  surface  of  the  scatterer. 


3.4  Numerical  Results  and  Discussions 

3.4.1  Normal  incidence 

Figure  3.3  shows  the  calculated  induced  surface  current  density 
distribution  on  the  strip  with  ka  - 4 (1.3X  wide)  for  normal  incidence. 
Note  that  the  current  density  becomes  large  at  the  edges,  as  it  should 
for  E-wave  incidence,  although  no  specific  condition  was  enforced  at 
the  edges,  nor  any  special  care  exercised.  Note  also  that  the  approx- 
imate current  is  confined  essentially  on  the  surface  of  the  strip  and 
extends  very  little  outside  of  this  surface.  Thus,  the  solution  in  this 
case  is  very  close  to  the  true  solution  and  this  is  easily  verified  by 
truncating  the  current  density,  computing  the  scattered  field  it  radiates 
on  the  strip,  and  verifying  that  the  scattered  field  is  Indeed  very 
nearly  equal  to  -E^. 
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Figure  3.4  depicts  the  result  for  ka  ■ 40,  l.e.,  a 13X  strip. 

Note  that  the  peak  in  the  center  is  no  longer  present  and  the  current 
there  approaches  that  given  by  the  physical  optics  approxlnatlon.  There 
are  now  more  oscillations,  however,  and  the  current  density  has  a sharp 
dip  before  rising  to  infinity  at  the  edges. 

Figure  3.5  displays  the  moment  method  applied  in  the  spectral  do- 
main solution  [ 3 ] and  the  comparison  with  the  one  obtained  here  is 
quite  favorable. 

Figure  3.6  exhibits  the  satisfaction  of  the  boundary  condition 
after  one  iteration.  As  mentioned  before,  such  a.  test  is  not  available 
in  the  conventional  GTD  approach. 

3.4.2  Near  grazinR  incidence 

Let  us  next  turn  to  the  Interesting  case  of  a near  grazing  inci- 
dence where  the  zero-order  current  density  has  a long  tall  extending 
beyond  the  edge  of  the  strip  (see  Figure  3.7).  This  result  is  to  be 
expected  since  the  two  half-plane  GTD  solutions  used  in  the  zero-order 
approximation  represent  a poor  approximation  for  the  Induced  current 
for  shallow  incidence  angles.  If  this  tall  is  truncated,  the  remain- 
ing portion  of  the  current  density  on  the  strip  produces  a scattered 
field  on  the  surface  of  the  strip  which  is  significantly  different  from 
-E^,  as  may  be  seen  from  Figure  3.8. 

Figure  3.9  shows  the  effect  of  one  iteration  on  the  zero-order 
GTD  solution  shown  in  Figure  3.7.  Note  that  the  current  density  is 
significantly  altered  in  the  neighborhood  of  the  shadowed  edge  demonstra- 
ting the  fact  that  even  with  a relatively  poor  initial  guess,  the  con- 
vergence is  quite  rapid  in  this  case. 


INDUCED  CURRENT  (AMP/M) 


20 


NORMALIZED  DISTANCE  D/A 

Figure  3.5.  Moment  method  (applied  in  the  spectral  domain)  solution 
of  the  magnitude  of  the  Induced  surface  current  density 
distribution  normalized  to  1/Z  on  the  strip  of  ka  - 50. 
(15.92  X wide),  - 90°.  ° 


Figure  3.6.  Magnitude  of  the  scattered  E-field  evaluated  on  the 
strip  of  ka  ■ 40.,  ■ 90°  (one  iteration). 
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To  see  that  this  Is  Indeed  an  improved  solution,  the  truncated 
portion  of  It  Is  used  to  calculate  the  scattered  field.  It  Is  observed 
that  the  satisfaction  of  the  boundary  condition  has  been  Improved  as 
shown  In  Figure  3.10. 

To  verify  the  convergence  of  the  solution  numerically,  one  more 
Iteration  Is  performed  and  the  result  Is  depicted  In  Figure  3.11.  Note 
that  the  shape  of  the  surface  current  density  does  not  change  much  which 
Indicates  a settling  down  of  the  solution  has  occurred.  Also,  note 
that  the  tall  extending  outside  of  the  atrip  has  been  reduced  to  an  In- 
significant quantity,  which,  when  truncated,  will  produce  little  effect 
on  the  scattered  field  on  the  surface  of  the  strip. 

To  further  validate  the  solution,  the  moment  method  solution  [ 3 ] 
of  the  same  problem  with  slightly  different  parameters  Is  shown  In  Figure 
3.12  for  a comparison.  Again,  the  agreement  Is  good.  However,  In  terms 
of  computational  efficiency,  the  present  method  Is  far  superior  to  the 
moment-method  solution  for  the  accuracy  realized. 

3.5  Summary 

Before  closing  this  chapter.  It  is  worthwhile  to  recapitulate  the 
main  points  of  the  approach  discussed.  The  strip  problem  has  been  solved 
by  a combination  of  the  Integral  equation  and  asymptotic  high-frequency 
techniques.  Formulation  of  the  Integral  equation  In  the  Fourier  trans- 
form domain  allows  one  to  conveniently  obtain  the  zero-order  approximation 
to  the  transformed  unknown  surface  current  density  from  the  solution  of 
two  half-plane  problems. 
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Figure  3,10. 


Magnitude  of  the  scattered  E-fleld  evaluated  on  the 
strip  of  ka  - 40.,  - 10°  (one  iteration). 
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NORMALIZED  DISTANCE  D/A 

Figure  3.12.  Moment  method  (applied  in  the  spectral  domain)  solution 
of  the  magnitude  of  the  Induced  surface  current  density 
distribution  normalized  to  1/Z_  on  the  strip  of 
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Higher-order  solutions  have  been  obtained  via  the  Iteration  steps 
outlined  above  and  the  numerical  convergence  has  been  demonstrated. 

The  iteration  process  generates  the  proper  edge  singularities  even  when 
they  are  not  present  in  the  original  approximation,  e.g.,  physical  optics. 
However,  additional  iterations  are  necessary  in  that  case.  Validity 
of  the  solution  has  been  substantiated  by  numerically  verifying  the 
satisfaction  of  the  boundary  condition. 
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4.  DIFFRACTION  BY  A FINITE  THIN  PLATE 


4.1  Introduction 

Having  Illustrated  the  usefulness  of  the  hybrid  technique  for  com- 
bining the  Integral  equation  and  GTD  techniques  relevant  to  a two- 
dimensional  scatterer,  viz.,  the  strip,  we  now  turn  to  the  more  general 
three-dimensional  problem,  a thin  rectangular  plate  Illuminated  by  a 
plane  wave.  This  problem  was  chosen  for  the  following  reasons:  In 

the  plate  problem,  the  difficulty  In  applying  GTD  to  this  geometry  stems 
from  the  fact  that  the  diffraction  coefficient  for  the  corners  of  the 
plate  Is  not  known  and  neglecting  the  corner  effects  can  cause  sub- 
stantial errors  In  the  resonance  frequency  region  where  the  plate  size 
Is  of  the  order  of  one  wavelength  squared. 

Before  we  discuss  the  present  approach  to  solve  the  plate  problem, 
we  would  like  to  mention  two  other  approaches  [14],  [4  ] that  are  based 
on  a combination  of  asymptotic  and  Integral  equation  techniques.  The 
one  developed  by  Thiele  [14]  decomposes  a given  problem  Into  two  parts, 
one  of  which  Is  handled  using  the  GTD  method  and  the  other  using  the 

moment  method.  For  the  case  of  a wire  antenna  on  a finite  ground  plane, 

the  effect  of  the  edge  diffraction  from  the  ground  plane  Is  evaluated 
using  GTD  and  the  result  Is  subsequently  used  to  augment  the  Impedance 
matrix  of  the  monopole  antenna  over  an  Infinite  ground  plane.  Although 
the  method  works  rather  well  when  GTD  results  are  accurately  known  for 

the  ground  plane  problem,  e.g.,  a ground  plane  of  circular  shape,  no 

convenient  means  Is  available  for  Improving  the  solution  when  there  are 
corners  In  the  plane  that  contribute  substantially  to  the  scattered  field. 
The  latter  situation  arises  when  the  ground  plane  Is  of  rectangular  shape 
and  Is  not  large  compared  to  the  wavelength,  or  when  the  antenna  Is 
mounted  close  to  one  or  more  of  the  edges. 
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The  second  method  developed  by  Burnside  [4  ] tends  to  rectify  the 
situation  alluded  to  above  by  solving  for  surface  currents  via  the  moment 
method  in  the  regions  where  the  GTD  solution  is  not  accurate,  and  by 
using  asymptotic  forms  for  the  surface  currents  in  regions  where  a good 
approximation  for  these  currents  can  be  employed.  However,  this  method 
cannot  be  conveniently  applied  to  either  the  strip  problem  with  grazing 
incidence,  or  to  the  large  plate  problem  discussed  in  this  thesis.  For 
the  strip  problem,  the  GTD  solution  is  quite  Inaccurate  when  the  incident 
angle  of  the  Illuminating  wave  is  near  grazing.  For  the  plate  problem, 
the  current  does  not  settle  down  to  known  asymptotic  form  in  the  center 
region  of  the  plate  until  it  is  at  least  three  to  four  wavelengths  squared. 
The  moment  method  is  incapable  of  handling  the  number  of  unknowns  required 
to  accurately  solve  for  the  current  distribution  on  plate  sizes  that 
are  larger  than  2X  squared. 

In  the  following  sections,  we  present  a new  approach  using  the 
iteration  method  developed  in  Chapter  2 to  demonstrate  that  a successful 
computation  of  the  Induced  surface  currents  on  the  plate  is  not  only 
feasible,  but  also  numerically  efficient. 

4.2  Iteration  Method  Applied  to  the  Plate  Problem 

For  the  sake  of  simplicity  we  consider  only  the  case  of  an  x-polarlzed 
uniform  plane  wave  which  is  normally  incident  on  a square  thin  plate. 

The  geometry  of  the  problem  is  depicted  in  Figure  4.1,  where  the  plate 
is  located  in  the  z - 0 plane. 

Using  classical  electromagnetic  theory,  the  following  coupled 

Integro-differentlal  equations  for  the  current  components  J and  J are 

X y 

readily  obtained: 
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Figure  4. 


Figure  4, 


. Diffraction  by  a finite  rectangular  thin  plate  Illuminated 
by  a normally  Incident  plane  wave  with  polarization  as  shown. 


. Regions  In  the  z ■ 0 plane  In  each  of  which  the  zero-order 
approximation  of  the  scattered  field  Is  obtained  according 
to  Table  4.1. 


Ay(x,y) 


and 


0 
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(4.1a) 


(4.1b) 


where  x e (-a, a),  y e (-b,b),  and  z - 0.  A and  A in  (4.1a)  and  (4.1b) 

X y 

are  the  x-  and  the  y-componenta  of  the  magnetic  vector  potential,  re- 
spectively, containing  the  unknowns  J ,J  implicitly. 

X y 

Since  the  convolution  of  the  induced  surface  current  density  with 
the  free-space  Green's  function  gives  the  magnetic  vector  potential,  we 
have  the  expressions  which  are  yalld  in  the  z ■ 0 plane 


Ax(x.y)  * * G(x,y) 


(4.2a) 


and 


Ay(x,y)  - Jy(x,y)  * G(x,y) 


(4.2b) 


where  * denotes  the  convolution  operation  and  the  free-space  Green's 
function  is  given  by 


G(x  - x',y  - y') 


1 exp(ikr) 
4tt  r 


(4.3) 


where 

r - - x')^  + (y  - y')2  , z,z'  - 0 . 
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Note  that  (4.1a),  (4.1b)  are  conventional  integro-dlfferentlal  equations 

which  are  valid  on  the  plate  only.  To  obtain  an  extended  form  of  these 

equations,  additional  unknown  functions  F (x,y)  and  F (x,y)  will  be  in- 

X y 

troduced.  The  domain  of  these  functions  is  the  region  complementary 
to  the  plate  in  the  z ■ 0 plane.  Hence,  the  extended  form  of  (4.1a), 
(4.1b)  can  be  written  as 


2 

— r + k^lA^(x,y)  + 
i3x 


— Ay(x,y)  - -itoe^ 


(-6(Ej^(x,y)) 

+ 0(F^(x,y))j 


and 


Ay(x,y) 


A (x,y) 
3x3y  X 


-iueQ0 (Fy(x,y)) 


(4.4a) 


(4.4b) 


where  9 and  0 are  operators  defined  In  (2.2a)  and  (2.2b),  respectively. 
Note  that  (4.4a),  (4.4b)  are  valid  for  the  entire  z ■ 0 plane.  These 
equations  are  now  Fourier  transformed  to  obtain 


(a^  - k^)J  (o,e)G(o,B)  + aPJ  (o,e)G(o,0) 
X y 


- “ (-0(ES(a,3)  + e(F  )(o,3))  (4.5a) 

4q  X X 


and 


ikjj  ^ 

(3  - k^)J  (o,3)G(a,3)  + o3J^(a,B)G(a,3)  9(F  )(a,B) 

y X y 


(4.5b) 
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where  ■>  on  top  Indicates  the  Fourier  transform  as  defined  In  (3.6a) 
with  transform  variables  (a, 6)  corresponding  to  (x,y),  respectively, 


and  Zq  > /Uq/Cq  Is  the  free-space  Impedance.  In  writing  (4.5a)  and 
(4.5b),  (4.2a)  and  (4.2b)  have  been  utilized.  The  Fourier  transform 
of  the  free-space  Green's  function,  specialized  to  the  z - 0 plane  as 
given  In  (4.3),  Is 


G(o,6)  - f • 


1.2  2 

7^0  - a - 6 


(4.6) 


Observe  that  (4.5a),  (4.5b)  are  two  algebraic  equations  In  the  transform 

domain  as  opposed  to  the  two  Integro-dlfferentlal  equations  In  (4.4a) 

and  (4.4b).  It  Is  a simple  step  to  derive  the  zero-order  solutions  of 

J (a, 6)  and  J (a, 3)  once  the  estimates  of  6(F  )(a,3)  and  6(F  )(a,3)  are 
X y X y 

available.  One  merely  solves  the  two  coupled  algebraic  equations  for 

these  two  unknowns  J (o,6)  and  J (o,6).  For  the  present  case,  6(F  )(x,y) 

X y y 

is  zero  due  to  the  particular  choice  of  x-polarlzed  normal  incident  plane- 
wave  illumination.  With  this  In  mind,  the  first-order  solutions  of  the 
transformed  surface  current  density  can  be  expressed  as 


■ [T  '"'J — ^ 


2 2 
^ -^0 


*^0^0 


1-9 (E^  (a, 8)  + 0(F^)(o,8)] 


(4.7a) 


and 


Jy(a,6) 


-2 


1? 2" 

7^0 " “ ■ 


r tts 

7 


0(Ej)(a,8)  + e(F^)(cx,6)] 


(4.7b) 


i 


i 
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J (a,  6)  and  J (a, 6)  are  then  inverse  Fourier-transformed  and  truncated 
X y 

to  obtain  the  induced  surface  current  densities  on  the  plate. 

Next,  we  estimate  the  zero-order  approximation  to  @(F^)(x,y')  using 
the  GTD  solutions  to  four  pertinent  half-plane  problems.  The  z ■ 0 plane 
containing  the  plate  has  been  divided  into  regions  as  shown  in  Figure  4.2 
(page  32),  where  the  hatched  region  is  occupied  by  the  plate  and  the  scattered 
field  in  this  region  must  be  equal  to  the  -6(E^(x,y))  to  satisfy  the 
boundary  condition.  The  rest  of  the  z 0 plane  hds^been  designated 
by  digital  numbers  acid  Roman  numerals,  and  the  manner  in  which  the 
scattered  fields  in  these  various  regions  are  obtained  is  concisely 
tabulated  in  Tabl<;  4.1. 


TABLE  4.1 

THE  ZERO-ORDER  APPROXIMATION  OF  6(F  (x,y))  IN  VARIOUS  REGIONS 
EXTERNAL  TO  THE  PLATE  OBTAINED  VIA  THE  USE  OF  GTD  METHOD 


Region 

Scattered  Field 

Remarks 

I.  II 

E (x,z  - 0) 

Derived  from  H-wave 

X 

strip  GTD  solution 

III,  IV 

E„(y.z  ■ 0) 

Derived  from  E-wave 

strip  GTD  solution 

1,  2, 

0 

First  estimate 

3,  4 

In  deriving  the  zero-order  approximation,  the  eoattered  fields  in 
regions  1,  2,  3 and  4 are  neglected  although  these  fields  are  nonzero 
in  higher-order  approximations. 
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The  zero-order  approximation  to  the  scattered  fields  Is  computed 
In  Regions  III  and  IV,  by  starting  with  the  E-wave  GTD  solution  for 
the  strip,  and  truncating  It  so  that  It  is  nonzero  only  In  these  regions. 
In  partlcul^.r. 


and 


+ 1 


2 sln(aa) 

6 J a 


1 


2 aln(aa) 

C( 


(4.8a) 


(4.8b) 


In  Regions  I and  II,  the  H-wave  GTD  solution  for  a strip  Is  used  to 

obtain  H (x,z  » 0),  and  then  E (x,z  - 0)  Is  constructed  from  Maxwell's 
y x 

equations.  The  resulting  solution  Is  again  truncated  so  that  It  Is 
nonzero  only  In  the  appropriate  regions.  In  particular. 


-i/k  + a e^®**  . 

7^ r + 


iJv.^  - 2 sln(Bb) 


“J 


e 


and 
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Having  completed  the  estimation  of  the  zero-order  approximation 
to  the  scattered  field  6(F^(x,y))  external  to  the  plate,  we  now  proceed 
to  solve  for  the  Induced  surface  current  on  the  plate.  To  this  end  we 
return  to  (4.7a)  and  (4.7b)  and  substitute  the  Fourier  transform  of 
0(F^(x,y))  and  compute  and  J^.  The  desired  Induced  surface  current 
densities  In  the  space  domain  are  then  obtained  by  Inverse  Fourier 
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transformation  and  truncation.  If  necessary,  the  Iteration  scheme  dis- 
cussed In  the  previous  sections  can  be  followed  to  obtain  higher-order 
solutions.  Convergence  of  the  solution  can  be  checked  by  performing 
one  more  Iteration  and  checking  to  see  whether  the  solution  has 
"settled  down."  Validity  of  the  solution  can  be  assured  by  computing 
the  scattered  field  on  the  plate  using  the  solution  of  the  surface  cur- 
rent just  obtained  to  see  how  well  the  boundary  condition  Is  satisfied 
on  the  surface  of  the  plate.  It  should  be  ''lear  now  that  all  of  these 
steps  follow  exactly  the  same  line  as  In  the  case  of  the  strip  problem 
discussed  In  Chapter  3. 

4.3  Numerical  Results  and  Discussions 

The  numerical  result  for  the  dominant  x-component  of  the  surface 
current  density  for  a one-wavelength  squared  plate  Is  shown  In  Figure 
4.3a.  Note  that  the  surface  current  density,  which  goes  to  zero  at  the 
two  edges  perpendicular  to  the  Incident  electric  Intensity  vector,  tends 
to  grow  without  bound  at  the  other  two  edges  parallel  to  the  Incident 
electric  Intensity  vector,  although  no  special  edge  condition  has  been 
enforced  to  derive  this  behavior.  To  see  that  edge  behavior  better,  a 
90°  rotation  of  the  surface  current  In  Figure  4.3a  Is  shown  In  Figure 
4.3b.  It  Is  clearly  seen  from  Figure  4.3b  that  the  cross  section  of 
the  x-component  of  the  surface  current  density  at  x » 0 closely  re- 
sembles the  surface  current  on  the  strip  plotted  In  Figure  3.3. 

Figures  4.4a,  4.4b,  and  4.4c  exhibit  the  change  In  the  behavior 
of  the  current  distribution  both  In  the  middle  of  the  plate  and  at  the 
corners  as  the  plate  size  Is  progressively  Increased. 


Figure  A. 4b.  Magnitude  of  the  donlnant  x-coaponent  of  the  surface  cur 
density  on  a 3X  x 3A  plate  (ka  - 9.43);  plate  region: 

yE(-l,l),  noraal  Incidence  with  x-polarlzatlon. 
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Figure  4.4c.  Magnitude  of  the  dominant  x— component  of  the  surface  current 
density  on  a 4X  x 4X  plate  (ka  = 12.6);  plate  region: 

ye (-1,1),  normal  incidence  with  x-polarization. 
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The  corner  singularities,  which  are  highly  localized,  appear  to  be  | 

present  only  for  relatively  small  sized  plates  where  corner-corner  inter-  I 

actions  may  play  a significant  role.  It  has  not  been  possible  to  verify  ^ 

I 

the  correctness  of  the  behavior  of  the  solution  at  the  corners  since  no 

f 

reliable  analytical  or  umerlcal  solutions  for  the  comer  problem  are  | 

available  at  this  time.  | 

Figures  4.5a  and  4.5b  show  the  cross-polarized  component  of  the  | 

surface  current  density  on  the  plate.  This  current  density  goes  to  zero 
at  the  line  of  symmetry  in  the  middle  of  the  plate  and  has  a tendency 
to  grow  without  bound  at  the  edges.  The  results  for  the  one-wavelength 
squared  plate  have  been  checked  by  moment-method  solutions  and  the  agree- 
ment is  good.  For  such  an  electrically  small  plate,  results  are  avail- 
able for  comparison  in  the  literature  [15],  [16].  However,  for  elec- 
trically large  plates  the  matrix  size  becomes  prohibitively  large  when 
the  conventional  moment  method  type  of  approach  is  used.  In  contrast, 
the  accuracy  and  the  convergence  of  the  solution  Improve  even  further 
for  a large  scatterer.  It  should  also  be  mentioned  that  the  number  of 
grid  points  at  which  the  current  density  of  the  plate  has  been  evaluated 
is  2048.  Such  fine  details  of  the  current  behavior  would  also  be  im- 
practical to  obtain  using  the  moment  method. 
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Magnitude  of  the  cross-polarized  component  of  the  sur- 
face current  density  on  a lA  x IX  plate  (ka  ■ 3.14); 
plate  region:  xe(-l,l),  ye (-1,1),  normal  incidence  with 
x-polarlzatlon. 
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Figure  4.5b.  Magnitude  of  the  cross-polarized  component  of  the  sur- 
face current  density  on  a 3X  x 3X  plate  (ka  « 9.43); 
plate  region:  xe(-l,l),  ye (-1,1),  normal  incidence  with 
x-polarlzatlon. 


5.  A NEW  LOOK  AT  THE  SCATTERING  OF  A PLANE 
WAVE  BY  A RECTANGULAR  CYLINDER 
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5.1  Introduction 

The  problem  of  scattering  of  electromagnetic  waves  by  a rectangular 
cylinder  has  been  Investigated  by  numerous  scientific  researchers  In 
the  past;  yet  Its  flavor  has  not  diminished  but  quite  on  the  contrary, 

It  stimulates  a strong  Interest  In  further  pursuing  this  problem  lately 
as  evidenced  by  the  number  of  papers  appeared  in  the  literature  [4],  [17], 
[18],  [19],  [20],  [21]  that  deal  with  this  subject. 

Due  to  the  two-dimensional  nature  of  the  problem,  the  conventional 
moment  method  [17]  can  be  applied  In  this  case  to  an  electrically  large 
cylinder  when  Its  symmetric  properties  are  also  exploited.  Therefore, 

It  provides  a convenient  validity  check  for  other  methods.  Indeed, 
Burnside  et  al.  [4]  and  Mautz  et  at,  [18]  have  claimed  their  different 
approaches  valid  only  by  a comparison  of  results  with  the  conventional 
moment  method. 

Both  of  the  aforementioned  papers  made  use  of  the  so-called  uniform 
geometrical  theory  of  diffraction  developed  by  Kouyoumjlan  et  al.  [8]  In 
which  the  classical  wedge  diffraction  coefficient  put  forth  by  Keller  has 
been  modified  using  multiplicative  factors  in  such  a manner  that  one  of 
these  factors  goes  to  zero  as  the  Keller's  diffraction  coefficient  goes 
to  Infinity  at  one  of  the  shadow  or  reflection  boundaries  so  that  the 
product  remains  finite.  The  resulting  diffraction  coefficient  is  Indeed 
applicable  at  all  aspect  angles.  It  is  conceivable  that  an  additive 
term  which  goes  to  Infinity  as  the  Keller's  diffraction  coefficient  goes 
to  Infinity  but  with  opposite  sign  will  work  equally  well  as  does  the 
multiplicative-factor  type  of  remedy.  Indeed,  this  addltlve-term  type 


of  remedy  has  been  followed  by  Lewis  et  al.  [22].  Both  of  these  remedies 
multiplicative  and  additive  alike,  still  have  room  for  further  develop- 
ment, and  it  Is  difficult  to  Judge  which  one  is  superior.  Furthermore, 
both  of  these  theories  'are  based  on  an  Ansatz  that  has  no  physical  Justl- 

I 

flcatlon  at  all  to  starts  with. 

In  this  chapter,  we  tackle  the  scattering  problem  using  a completely 
different  technique  which  Is  based  on  a representation  of  the  scattered 
fields  In  terms  of  the  spectrum  of  the  Induced  surface  current  on  the 
scatterer  rather  than  the  rays  emanating  from  It.  This  spectral  domain 
interpretation  of  high-frequency  diffraction  phenomena  has  been  docu- 
mented In  detail  In  a recent  paper  by  Mlttra,  Rahmat-Samll  and  Ko  [10]. 

Ko  and  Mlttra  [23]  also  developed  a method  based  on  the  spectral  domain 
concept  for  combining  the  asymptotic  high-frequency  technique  and  the 
Integral  equation  formulation.  The  method  has  been  applied  successfully 
to  Infinitesimally  thin  scattering  objects,  e.g.,  an  infinite  perfectly 
conducting  strip  and  a finite  perfectly  conducting  square  plate  [23]. 

To  further  develop  the  method,  we  apply  It  to  the  perfectly  conducting 
square  cylinder. 

This  particular  geometry  has  been  chosen  because  it  serves  as  a good 
example  for  demonstrating  that  the  method  can  be  extended  to  treat  with 
ease  scattering  objects  with  finite  thickness  and  planar  surfaces. 
Moreover,  the  validity  of  the  method  can  be  easily  established  by  com- 
paring with  well-established  solutions  derived  by  entirely  different 
approaches  available  In  the  literature.  We  show  that  starting  out  with 
Keller's  wedge  diffraction  coefficient,  only  one  iteration  gives  a far- 
fleld  pattern  that  compares  extremely  well  with  that  obtained  from  the 
moment-method  solution  for  thirty-two  unknowns. 
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To  emphasize  the  simplicity  and  efficiency  of  the  present  approach. 

It  is  worth  mentioning  that  the  GTD  result  deviates  significantly  In 
those  directions  which  coincide  with  the  orientations  of  the  surfaces  of 
the  square  cylinder  even  when  multiple-edge  Interactions  are  Included 
In  a self-consistent  manner  by  Introducing  eight  unknown  diffraction  co- 
efficients. Only  when  one  further  modifies  the  GTD  method  with  the 
uniform  diffraction  coefficient  and  combines  It  with  the  moment  method 
using  a total  of  twenty-four  unknowns  will  the  result  compare  well  with 
the  moment-method  solution  for  thirty-two  unknowns. 

The  present  approach  shows  how  the  Keller's  GTD  solution  can  be  Im- 
proved In  a straightforward  manner  based  on  a physical  argument  rather 
than  a remedy  from  an  Ansatz.  Thus,  It  gives  Insight  Into  the  mechanism 
involved  in  problems  of  high-frequency  scattering  from  scatterers  with 
planar  facets. 

5.2  Formulation 

The  geometry  of  the  electromagnetic  scattering  problem  involving  a 
perfectly  conducting  Infinite  rectangular  cylinder  of  cross  section  2a  x 2b 
Illuminated  by  a uniform  plane  wave,  \diose  magnetic  Intensity  vector  Is 
oriented  parallel  to  the  edges  of  the  cylinder.  Is  depicted  In  Figure  5.1. 
For  convenience  of  analysis,  an  arbitrary  incident  wave  can  always  be  de- 
composed Into  two  components  with  respect  to  the  z-axls,  namely,  TM 

z 

(E-wave)  and  TE^  (H-wave) . In  the  following  discussion  we  consider  the 

H-wave  case  only;  the  E-wave  case  can  be  solved  In  a similar  manner  by 

-1  ~ 1 
considering  E ■ zEq. 

In  the  H-wave  case,  the  Incident  H-field  is  given  by 

^ -ik(xcos(|)Q+ysin(JiQ) 

HgCp.'t')  ■ e 


(5.1) 
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where  the  e time  dependence  Is  understood.  The  geometrical  optics 
reflected  field  can  be  derived  easily  from  (5.1)  once  the  geometrical 
configuration  Is  known.  Hence,  the  geometrical  optics  part  of  the  total 
field  will  not  be  discussed  on  further. 

We  now  turn  to  the  diffracted  field  which  results  from  the  sharp 
edges  on  the  scattering  cylinder.  As  alluded  to  In  Section  5.1,  we  make 
use  of  the  Keller  wedge  diffraction  coefficient  to  find  the  diffracted 
fields.  To  make  this  discussion  as  self-contained  as  possible  and  to 
Introduce  the  notations  that  we  have  adopted  In  this  chapter,  a brief 
review  of  the  wedge  diffraction  coefficient  Is  In  order. 

The  geometry  of  a perfectly  conducting  wedge  Immersed  In  a uniform 
H-wave  In  the  canonical  wedge  diffraction  problem  Is  shown  In  Figure  5.2. 
According  to  Keller's  geometrical  theory  of  diffraction,  the  diffracted 
field  by  the  edge  of  the  wedge  can  be  constructed  by  the  following  formula: 


H 


d 


z 


H-* 


(at  the  edge  of  the  wedge) 


(5.2) 


where 


and  ~ means  "asymptotically  equal  to."  In  the  above  formula,  nir  denotes 
the  exterior  region  while  (2  - n)ir  denotes  the  wedge  angle.  Note  that  n 
Is  not  necessarily  an  Integer;  In  the  rectangular  cylinder  case.  It  will 
take  on  the  value  3/2.  The  subscript  w attached  to  the  angle  of  Incidence 
(||q  and  the  angle  of  diffraction  4i  serves  as  a reminder  that  these  angles 
are  used  In  the  canonical  problem  and  their  senses  are  defined  as  shown 
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In  Figure  5.2.  This  fact  Is  worth  emphasizing  since  It  Is  crucial  In 
constructing  the  diffracted  field  when  several  wedges  are  Involved  In 
a scattering  problem. 

As  shown  In  Figure  5.3,  the  major  contributions  to  the  diffracted 
field  are  from  wedges  1,  2,  and  3 since  the  three  edges  of  these  wedges 
are  being  Illuminated  by  the  Incident  field  whereas  wedge  4 Is  In  the 
dark  and  hence  It  can  be  Ignored  In  the  zeroth-order  approximation  to 
the  diffracted  field.  With  this  In  mind,  the  diffracted  far  field  can 
be  written  In  a concise  fashion  as  follows: 


(5.4) 
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where  ({),.  ■ ir/4.  In  the  above  equation,  subscripts  1,  2,  and  3 on 

u z 

designate  wedges  1,  2,  and  3,  respectively.  Each  of  the  H^'s  from 

z 

the  Individual  wedges  Is  of  the  form  given  by  (5.2)  and  (5.3)  with 
proper  values  for  ({i^  and  substituted  and  proper  phase  shifts  taken 
care  of  due  to  the  transfer  of  the  phase  center  located  at  the  edge 
of  the  wedge  In  each  of  the  Individual  canonical  problems  to  the  common 
phase  center  located  at  the  origin  of  the  rectangular  coordinate  system 
shown  In  Figure  5.3. 

5.2.1  Pole  singularities  In  the  diffraction  coefficients 
It  Is  well-known  that  the  Keller's  diffraction  coefficient  (5.3) 
as  Introduced  In  the  last  section  will  not  give  the  correct  value  for 
the  diffracted  field  at  certain  directions,  namely,  at  the  shadow  and 
the  reflection  boundaries  predicted  by  geometrical  optics.  One  can  easily 
see  from  (5.3)  that  the  denominator  of  either  one  of  the  two  terms  In 
the  square  brackets  vanishes  at  one  of  these  shadow  and  reflection 
boundaries  and  gives  Infinite  field  values  at  those  directions.  We 
know  that  the  field  actually  remains  finite  and  varies  continuously 
across  the  shadow  and  reflection  boundaries.  Hence,  we  cannot  apply 
Keller's  diffraction  coefficient  near  the  shadow  and  reflection  bound- 
aries. It  was  precisely  for  this  reason  that  some  work  had  been  done 
In  trying  to  construct  a uniform  diffraction  coefficient.  As  mentioned 
earlier  In  the  Introduction  of  this  chapter,  both  multiplicative-  and 
addltlve-type  approaches  have  been  pursued  and  have  shown  some  success. 

In  the  following  discussions  we  show  that  It  Is  not  necessary  to  use 
the  aforementioned  "uniform  diffraction  coefficient"  to  solve  the  rec- 
tangular cylinder  problem  under  consideration. 
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{ To  start,  we  make  direct  use  of  Keller's  diffraction  coefficient 

; as  given  In  (5.3)  and  substitute  It  In  (5.2)  to  find  the  wedge  dlf- 

I 

fracted  field.  Observing  that  the  diffraction  coefficient  contains 
: pole  singularities  which  give  nonphysical  fields  at  the  shadow  and 

reflection  boundaries,  we  will  first  subtract  out  the  pole  singularities. 

5 

The  reason  for  subtracting  out  the  pole  singularities  becomes  clear 

i 

In  the  next  section  where  these  singularities  In  the  diffraction  co- 
efficient are  shown  to  correspond  exactly  to  those  In  the  physical 
optics  current  that  exists  on  the  seml-lnfinlte  surface  of  the  wedge 
In  the  canonical  problem  but  should  not  have  existed  In  the  rectangular 
cylinder  case  since  the  surfaces  of  the  rectangular  cylinder  are  finite 
In  width.  With  this  In  mind,  the  expressions  for  the  zeroth-order  approx- 
imation to  the  diffracted  field  can  be  written  explicitly  as  follows: 

For  wedge  1 (see  Figure  5.3), 
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where 


0 1 <(>0  — "f  0 £ <(*  f.  2^ 


Note  that  the  second  term  In  the  square  brackets  Is  the  pole  singu- 
larity of  the  first  term  at  the  shadow  boundary  and  the  fourth  term 
is  the  pole  singularity  of  the  third  term  at  the  reflection  boundary. 
For  wedge  2 (see  Figure  5.3), 
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Note  that  no  singularities  are  contained  in  the  first  term  in  the  given 
angular  domain  of  4).  The  third  term  in  the  square  brackets  is  the 
pole  singularity  of  the  second  term  at  the  reflection  boundary  for 
0 £ 4>  <,  TT  and  the  fourth  term  is  that  for  £ iji  £ 2ir.  Also,  note 
that  there  are  no  shadow  boundaries  in  this  case. 

For  wedge  3 (see  Figure  5.3), 


(at  edge  of  wedge  3) 
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Note  that  the  second  term  In  the  square  brackets  Is  the  pole  singularity 
of  the  first  term  at  the  shadow  boundary  and  the  fourth  term  Is  that 
of  the  third  term  at  the  reflection  boundary.  For  future  reference, 
all  the  pole  singularities  contained  in  the  diffraction  coefficient, 
which  have  been  discussed  In  detail  above  and  shown  explicitly  In 
(5.5a),  (5.5b)  and  (5.5c),  are  tabulated  In  Table  5.1. 

5.2.2  Pole  singularities  in  the  physical  optics  currents 

It  was  discussed  in  the  last  section  how  the  diffracted  far  field 
can  be  constructed  from  the  diffraction  coefficients.  In  applying 
these  diffraction  coefficients,  one  must  recognize  the  fact  that  they 
were  derived  from  solving  a canonical  problem  which  Involved  an  In- 
finite perfectly  conducting  wedge,  or  more  descriptively,  two  seml- 
Inflnlte  perfectly  conducting  half  planes  joining  each  other  at  an 
angle.  Therefore,  the  solution  of  the  canonical  problem  predicts  a 
physical  optics  current  existing  on  the  Illuminated  seml-lnflnlte 
surfaces  of  the  wedge.  Since  the  support  of  this  physical  optics 
current  Is  infinite.  It  gives  rise  to  infinite  fields  at  certain 
directions  In  space.  In  the  following  discussions,  we  show  that  all 
the  pole  singularities  contained  In  the  wedge  diffraction  coefficient 
are  precisely  due  to  these  physical  optics  currents  with  infinite 
supports. 

To  tackle  this  problem  one  step  at  a time.  Figure  5.3  has  been 
redrawn  in  more  detail  in  Figures  5.4a,  5.4b,  Figures  5.5a,  5.5b,  and 
Figures  5.6a,  5.6b.  These  figures  exhibit  the  physical  optics  currents  on  the 
Illuminated  surfaces  of  the  Illuminated  wedges  of  the  rectangular  cylinder  as 
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TABLE  5.1 


POLE  SINGULARITIES  AT  SHADOW  AND  REFLECTION  BOUNDARIES 
THAT  HAVE  BEEN  SUBTRACTED  FROM  THE  PERTINENT  DIFFRACTION  COEFFICIENTS 
(S.  B. : SHADOW  BOUNDARY,  R.  B. : REJECTION  BOUNDARY) 


Wedge 

Pole  at 

Pole  singularity 

1 

S.  B. 

ik  (-acos  (|)Q+b8  in(^Q  ) 
e 

({>  - (ii  + (|)q) 

1 

R.  B. 

ik(-aco8(j)Q-b8in^Q) 

e 

4>  - (ti  - <l>o) 

2 

R.  B.  1 

lkacos4>Q  -ikbsin(j>Q 
-e  e 

4>  - ir  + (J.Q 

2 

R.  B.  2 

-lkacos(|)Q  lkbsin(|),^ 
e e 

“ 2ir  + 

3 

R.  B. 

“ikacos(J)Q  -lkb8in(^Q 
-e  e 

4»  ~ 2ir  + (|)q 

3 

S.  B. 

lkacos(|iQ  -lkbsln4)Q 
-e  e 

(J)  - (tt  + ♦q) 

i 


well  as  the  shadow  and  reflection  boundaries  created  by  the  Illuminating 


plane  Incident  wave  as  predicted  by  classical  geometrical  optics. 

J.O 

Referring  to  Figure  5.4a,  the  far  field  due  to  the  physical  optics 

P 0 

current  on  semi- Infinite  surface  A of  wedge  1 can  be  obtained 

using  classical  electromagnetic  theory.  First  of  all,  the  subscripts 
and  superscripts  attached  to  J and  H must  be  understood  properly. 

The  first  subscript  Indicates  the  vector  component  of  the  quantity; 
the  second  subscript  designates  the  surface  on  which  the  current  exists; 
and  the  third  subscript  depicts  the  Illuminated  wedge  under  scrutlnl- 
zatlon.  The  superscript  reveals  that  the  quantity  under  consideration 
Is  obtained  from  physical  optics.  In  classical  electromagnetic  theory, 

I 

P 0 • P aO* 

the  far  field  H^^^^'and  the  surface  electric  current  radiating  It 

In  free  space  are  related  In  the  following  manner. 


1- 

„P.O.  e^^P  e ^ , . -,P.O. 

“zAl  ^ “^xAl 


(5. 


where 


I I 


-lk(xco.*„+b.ln*  ) 

2e  5(y  - b)  e e ^ 


dxdy 


-a 


which  Is  the  Fourier  transform  of  the  surface  electric  current  density 

P.O.  -1 

J^^'obtalned  by  physical  optics,  l.e.,  2h  x H . The  factor  of  2 Is 


present  because  the  original  problem,  which  Involved  an  electric  cur- 
rent radiating  In  the  presence  of  a perfectly  conducting  surface  on 
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which  It  resides,  Is  equivalent  to  the  problem  of  twice  the  electric 
current  radiating  In  free  space.  If  we  Ignore  all  of  the  effects  at 
Infinity,  which  can  be  done  by  Introducing  an  Infinitesimal  loss  In 
the  medium,  then  the  far  field  due  to  the  physical  optics  current  on 
surface  A of  wedge  1 can  be  expressed  as: 


1-^ 

Ikp  ^4  . 

-Ik  ^ jg  1 ^ - 2 sin  <|)  - 


- Ikb  sln(j>Q-lkbsln(|>  la(kcos(^  ^ ^ ^ 

1 e * 

l(k  cos  (jijj  + k cos  (j)) 


It  Is  very  Interesting  to  observe  that  the  denominator  In  (5.7)  vanishes 
at  ((i  IT  + which  are  precisely  the  reflection  and  shadow  boundaries 
as  shown  In  Figure  5.4b. 

Knowing  that  the  field  produced  by  the  physical  optics  current 
does  possess  pole  singularities  at  the  shadow  and  reflection  boundaries, 
let  us  check  the  singular  behavior  as  ())  approaches  To  this 

end,  the  denominator  will  be  expanded  Into  a Taylor's  series  In 
[<ji  - (it  + series,  only  those  terms  up  to  the  first 

power  of  (|)  will  be  kept  while  all  higher  power  terms  are  discarded 
since  their  contributions  are  small  compared  with  the  first  power  term 
In  the  series  as  approaches  Infinitely  close  to  tt  + i^q.  Putting 
this  In  terms  of  mathematical  language,  we  have 

H**:°*^. 
zAl  di4-(b 

y yp 

(5.8) 

where 


i 

t 

i 


i 


Through  straightforward  algebraic  manipulations,  we  have 


p Q -lkbsln<j>Q+lkacos(|iQ  gilv 

'^zAl  (JrHT+(J)Q  ® 
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‘1 

/2nk 


-lkb8ln(ir+(j)Q)+lkaco8  (ir+(j)Q) 


(|)  - (ir  T (|)q) 


(5.9) 


At  thl8  point,  one  may  go  back  to  (5.5a)  or  Table  5.1  to  aee  that  theae 
are  Indeed  exactly  the  pole  elngularltles  at  the  ahadow  and  reflection 
boundarlea  of  the  wedge  diffraction  coefficient 

Referring  to  Figure  5.5a  and  proceeding  along  similar  lines,  one 
can  obtain 


„P-0.  . 

"zA2  “ 


2 sin  (j>  e 


- Ikb  sln({iQ“lkb8ln(Ji  -la(  kcoa  i^Q+kcos  (p  )" 


-l(k  cos  ({iq  + k cos  p) 


(5.10) 


and 


1- 

„P-0. 

“zB2 


-lkacosiJ)Q-lkaco8(}i  -lb  (ksln(j)Q+ksln^) 
-2  cos  (t  e e 


l(k  sin  ({iq  + k sin  (}>) 


(5.11) 


p.O.  P.O 

Similarly,  the  singular  behavior  of  H^^2*  ^zB2*  shown  to 

be  as  follows: 


p Q -lkbsln<|iQ-lkacos(frQ  ^Ikp 

**zA2  ® /p 


‘I 
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TStHc 


-lkbsln(ir+(|)Q)-lkacos  (it+<^q) 


- (if  ^ <l»o) 


(5.12) 
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and 


(5.13) 


P.O. 

It  Is  very  satisfying  to  observe  that  the  pole  singularity  of  at 

* . V - Is  exactly  the  pole  singularity  at  the  reflection  boundary 

off  surface  A of  the  diffraction  coefficient  In  (5.5b)i  and  the 

pole  singularity  of  at  ♦ - 2x  - singularity 

at  the  reflection  boundary  off  surface  B of  tbs  diffraction  coefficient 

D In  (5.5b)  (see  Figure  5.5b).  These  pole  singularities  are  also 

P 0 • 

listed  in  Table  5.1.  The  pole  singularity  of  at  ♦ - ir  + (|)q 

cancels  the  pole  singularity  of  H^-°‘at  ^ exactly.  It  is 

also  interesting  to  note  that  the  direction  d)  - tt  + is  inside 

wedge  2. 

Finally,  referring  to  Figure  5.6a  and  proceeding  as  before,  one 
obtains 
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,P.O. 
zB3 


ikp 


2 cos  d e 


-Ikacost^Q-ikacos^  ib  (ksin(JiQ+ksin(|)) 


i(k  sin  d'n  + 


(5.14) 
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P.Ot 

The  singular  behavior  of  Is  studied  similarly  and  the  result  Is 

shown  below. 


-lkacos((iQ+lkbsln<jiQ  ^Ikp 


(5.15) 


Once  again,  by  referring  to  Figure  5.6b  and  (5.5c),  one  can  verify  that 

P.O. 

these  pole  singularities  In  are  Indeed  those  at  the  shadow  and 

reflection  boundaries  of  the  diffraction  coefficient  D.  .. 

hj 

To  recapitulate,  the  pole  singularities  contained  In  the  diffraction 
coefficient  at  the  shadow  and  reflection  boundaries  correspond  exactly 
to  the  pole  singularities  contained  In  the  physical  optics  currents 
with  Infinite  supports.  This  demonstration  explains  why  we  have  sub- 
tracted out  the  pole  singularities  In  the  diffraction  coefficient  before 
we  applied  It  to  compute  the  diffracted  far  field.  The  main  reason 
Is  that  we  are  applying  these  diffraction  coefficients  to  construct  the 
far  field  diffracted  by  a finite  two-dimensional  scatterer  where  there 
Is,  In  fact,  lack  of  all  these  physical  optics  currents  with  Infinite 
supports  which  would  have  existed  on  the  Illuminated  surfaces  of  the 
seml-lnflnlte  wedge  In  the  canonical  problem  used  to  Identify  these 
wedge  diffraction  coefficients.  After  these  pole  singularities  have 
been  subtracted  from  the  diffraction  coefficients,  the  resulting  expres- 
sions for  the  diffraction  coefficients  can  be  used  directly  In  the 
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construction  of  the  diffracted  far  field.  This  procedure  has  been  shown 
explicitly  In  (5.5a),  (5.5b),  and  (5.5c). 

The  superposition  of  ^z3*  in  manner 

Indicated  by  (5.4),  gives  the  diffracted  far  field  from  the  rectangular 
cylinder.  The  far-fleld  pattern  computed  by  using  (|)q  ■ ■|■,a  ■ b ■ IX, 

In  (5.4)  Is  depicted  In  Figure  5.7.  This  plot  shows  that  there  are 
no  sporadic  variations  near  the  shadow  and  reflection  boundaries  at 
225°,  135°,  and  315°.  This  pattern  has  already  been  In  very  good 
agreement  with  that  obtained  by  using  the  moment  method,  as  shown  In 
Figure  5.11.  By  good  agreement,  we  mean  that  all  the  locations  of 
the  peaks  and  nulls  are  close  to  the  right  positions  and  the  levels  are 
more  or  less  on  the  right  track.  Except  for  the  obvious  discontinuities 
In  the  pattern  at  the  angles  of  0°,  90°,  180°,  and  270°,  which  cor- 
respond to  the  directions  In  which  the  four  surfaces  of  the  rectangular 
cylinder  are  oriented.  It  Is  really  remarkable  how  good  a far-fleld 
pattern  a zeroth-order  GTD  solution  can  give!  A simple  physical  In- 
terpretation for  the  existence  of  these  discontinuities  can  be  found; 
and  their  elimination  Is  discussed  In  the  following  sections. 

Before  closing  this  section.  It  Is  noteworthy  that  due  to  the 
complete  symmetry  of  the  rectangular  cylinder,  the  pole  singularities  In 
the  wedge  diffraction  coefficients  need  not  be  subtracted  out  explicitly 
as  done  In  (5.5a),  (5.5b),  and  (5.5c)  before  we  use  them  to  compute 
the  far  field.  In  other  words,  the  wedge  diffraction  coefficient  as 
given  In  Its  original  form  shown  In  (5.3)  could  have  been  used  directly 
to  compute  the  far-fleld  pattern;  and  Indeed  this  Is  done  and  the  result 
Is  plotted  In  Figure  5.8.  It  Is  Identical  to  Figure  5.7.  The  automatic 
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Figure  5.7.  Diffracted  far-field  pattern  of 
cylinder  obtained  by  using  (5.4 
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pole  cancellations  happened  In  this  completely  symmetric  case  are 
shown  schematically  In  Figure  5.3  where  pole  singularities  at  shadow  or 
reflection  boundaries  Indicated  by  like  arrow  pairs  annihilate  each 
other.  That  this  Is  Indeed  the  case  can  be  verified  by  checking  (5.5a), 
(5.5b),  and  (5.5c)  for  the  pole  singularities  enclosed  In  the  square 
brackets  as  well  as  paying  attention  to  the  phase  factor  outside  the 
square  brackets.  While  It  Is  admittedly  a happy  coincidence  that  the 
Keller's  diffraction  coefficients  can  be  used  directly  In  the  rectangular 
cylinder  problem.  In  more  general  cases  the  procedure  of  subtracting 
pole  singularities  explicitly  from  the  diffraction  coefficient,  as 
done  In  (5.5a),  (5.5b),  and  (5.5c),  should  be  followed  In  order  to  get 
meaningful  results  near  the  shadow  and  the  reflection  boundaries. 

5.2.3  Discontinuities  In  the  far-fleld  pattern  at  <|)  « 0,  , tt,  and 

In  the  previous  sections  the  problem  of  diffraction  by  a rectangular 
cylinder  was  solved  by  using  Keller's  wedge  diffraction  coefficients 
minus  the  pole  singularities  In  constructing  the  diffracted  far  field. 
Although  the  far-fleld  pattern  varies  smoothly  and  remains  finite  across 
the  shadow  and  reflection  boundaries,  there  are  noticeable  discontinuities 
In  the  pattern  at  4)  > 0,  ir,  and  These  directions  are  those  in 
which  the  surfaces  of  the  rectangular  cylinder  are  oriented.  In  the 
next  section,  a physical  Interpretation  of  the  existence  of  these  dis- 
continuities Is  given.  Before  doing  that,  these  discontinuities  are 
studied  and  the  results  are  displayed  In  full  detail  In  this  section. 

Let  us  refer  to  (5. A).  The  discontinuity  In  the  far  field  at 
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Following  the  same  lines,  the  discontinuity  in  the  far  field  at 
(J)  - is 


<(♦  • f -)  - - i +) 
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(5.16a) 


(5.16b) 


(5.16c) 


(5.16d) 
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(5.17a) 
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In  the  above  equation, 


7 -)  - lira 
^ ^ e-K)  ^ 


2 - 


and 


- f +) 


lim  ((^  - I + e) 
e-H)  ^ ^ 


where  e > 0 is  any  small  positive  real  number.  Similarly,  the  dis- 
continuity In  the  far-field  pattern  at  ^ * ir  is 
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(5.17b) 


(5.17c) 


(5.18a) 


(5.18b) 
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(5.19c) 


3ir 


Finally,  the  discontinuity  in  the  far-fleld  pattern  at  ■ y is 
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To  summarize,  the  discontinuities  in  the  far  field  at  - 0, 
y IT,  and  y have  been  studied  in  full  detail  and  the  results  documented 
in  this  section.  The  fact  that  the  GTD  solution  of  the  far  field  con- 
tains discontinuities  is  well-knovm.  In  the  next  section,  the  cause 
of  the  existence  of  these  discontinuities  is  investigated  and  a simple 
physical  interpretation  is  given  which,  in  turn,  provides  a clue  to 


their  elimination. 
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5.2.4  Physical  Interpretation  of  the  existence  of  the  dlscontl- 
nuities  In  the  far-field  pattern  at  <|»  ■ 0,  tt,  and  ^ 

In  the  last  section,  the  discontinuities  In  the  far  field  at 
7T  3ir 

())  >■  0,  ir,  and  — have  been  studied  and  documented.  The  fact  that 
the  GTD  solution  of  the  far  field  contains  discontinuities  Is  well- 
known.  The  cause  of  the  existence  of  these  discontinuities  Is  Investi- 
gated and  a simple  physical  Interpretation  Is  given  which.  In  turn, 

provides  a clue  to  their  elimination. 

Note  that  In  (5.4),  the  far  field  Is  constructed  by  using 

H^2*  ^z3'  which  has  Its  own  angular  domain  of  defini- 
tion. More  explicitly,  Is  defined  on  the  angular  range  0 £ ^ ^ 

(see  Figure  5.9a);  outside  of  this  angular  range,  simply  assumed 

to  be  zero.  Likewise,  defined  on  0 £ (j>  5 it  and  ^ 4*  2it 

d IT 

(see  Figure  5.9b);  and  is  on  0 £ (ji  £ y and  it  £ (Ji  £ 2tt  (see  Figure 

5.9c).  It  would  have  been  a valid  assumption  In  the  canonical  problem 
since  fields  would  have  existed  only  in  the  exterior  region  of  the  wedge. 

In  the  rectangular  cylinder  problem.  It  should  be  realized  that  the 
fields  do  not  terminate  abruptly  at  the  geometrical  planes  represented 
by  dashed  lines  In  Figures  5.9a,  5.9b,  and  5.9c.  These  dashed  planes  would 
have  been  occupied  by  the  seml-lnflnlte  wedge  surfaces  In  the  canonical  prob- 
lems. Physically,  the  fields  should  be  continuous  across  these  dashed  planes 
In  space.  Therefore,  In  constructing  the  far  field  as  shown  In  (5.4), 
we  have  effectively  created  a discontinuity  at  each  of  these  dashed 
planes.  The  presence  of  these  discontinuities  Is  solely  artificial 
because  they  are  created  by  our  abrupt  truncation  of  the  fields  to 
the  regions  In  space  corresponding  to  the  exterior  regions  of  the  wedges 
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Figure  5.9a.  3 


contribute  to  the  discontinuities  In  the 
far  field  at  ({>  - (j>  - 0,  respectively. 
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SURFACE  B 


Figure  5.9c.  J^3t  JyB3  contribute  to  the  discontinuities  in  the 
far  field  at  4 - ir,  <|i  “ 7,  respectively. 
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used  In  the  canonical  problems.  Now  that  the  cause  of  the  existence  of 
the  discontinuities  In  the  far-fleld  pattern  at  (|»  ■ 0,  ir,  and 
has  been  found,  we  proceed  to  give  them  a physical  Interpretation  which 
In  turn  provides  a clue  to  their  elimination. 

In  the  following  discussion,  a typical  case  Is  studied  and  the 
solution  to  this  case  Is  derived  In  a step-by-step  fashion.  While 
the  results  of  other  cases  will  be  exhibited,  their  derivations  are 
left  to  the  Interested  reader. 

The  typical  case  as  shovm  In  Figure  5.9a  Is  studied.  Let  us  con- 
centrate on  the  extension  of  surface  A.  This  extension  Is  represented 
by  a dashed  line  from  x*atox»»aty*bin  the  direction  ^ ■ 0. 

On  one  side  of  this  plane,  there  exists  as  shown  In  Figure  5.9a, 
on  the  other  side  there  Is  zero  field.  By  classical  electromagnetic 
theory,  we  conclude  that  there  Is  an  electric  surface  current  density 
existing  In  the  dashed  plane.  Such  a current  radiates  In  free  space 
to  support  this  discontinuity  In  the  H-fleld.  In  compact  mathematical 
language,  we  write 
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(5.21) 


where  h Is  the  unit  normal  vector  to  the  surface  A,  the  subscript  A 
refers  to  the  extension  of  surface  A,  the  subscript  1 designates  the 


wedge  under  consideration,  the  superscript  d means  diffraction,  and  the 


rest  should  be  self-explanatory.  Note  that  a factor  of  2 Is  not  present 

d 


In  (5.21)  because  this  Is  simply  the  discontinuity  In  free  space 

of  the  H-fleld  and  there  Is  no  backing  of  perfect  conductor  In  this 
situation. 

Referring  to  (5.5a)  and  keeping  In  mind  that  the  original  dif- 
fraction coefficient  Is  perfectly  valid  In  the  direction  ^ ■ 0,  which 
Inqplles  that  the  procedure  of  subtracting  the  pole  singularities  need 
not  concern  us  here,  we  can  write 
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Note  that  In  the  above  equation,  an  approximation  of  p by  x has  been 
made  and  a constant  value  zero  for  the  variable  ^ has  been  assumed. 
The  far-fleld  expression  as  given  In  (5.5a)  can  be  legitimately  used 
to  compute  the  H-fleld  on  the  dashed  plane  because  any  observation 
points  In  that  plane  will  be  at  least  at  a distance  of  2a  > 2A  away 
from  the  edge  of  wedge  1. 


The  H-field  generated  by  is 
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where  the  superscript  Ex  stands  for  "excess"  to  remind  the  reader  that 
this  field  has  been  created  artificially  by  abruptly  truncating  the 
diffracted  fields,  therefore,  It  should  not  have  existed;  hence,  the 
word  "excess."  (5.23)  Is  given  by 
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By  letting  x(k  - a)  ■ t^,  the  Integral  with  respect  to  x can 

be  transformed  Into  an  Integral  with  respect  to  t which  may  be  Identi- 
fied with  the  well-tabulated  Fresnel  Integrals,  l.e.. 
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where  ■ 0)  is  given  in  (5.22b). 


Now  we  have  the  expression  for  the  field  generated  by  as  shown 

in  (5.27).  The  next  step  is  to  examine  the  discontinuity  introduced 
at  (|)  “ 0 by  . To  this  end,  all  we  need  to  pay  attention  to  is 
the  following  limit: 
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Using  (5.28)  and  (5.27),  the  discontinuity  introduced  at  (|>  ■ 0 by 
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It  is  Indeed  enlightening  to  see  that  this  discontinuity  is  exactly 
the  discontinuity  given  in  (5.16e).  The  study  of  the  typical  case 
is  now  completed.  For  the  cases  depicted  in  Figures  5.9a,  5.9b,  and 
5.9c,  we  simply  list  the  results  as  follows: 

The  H-fleld  generated  by  In  Figure  5.9a  is 
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The  H-field  generated  by  Figure  5.9b  is 
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The  H-field  generated  by  j“g2  Figure  5.9b  is 
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The  H-field  generated  by  Figure  5.9c  is 
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The  H-fleld  generated  by  3°^^  Figure  5.9c  is 
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To  examine  the  discontinuity  introduced  at  4)  * ^ by 
(5.33)  to  obtain 
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Again,  this  discontinuity  is  identical  to  that  given  in  (5.17c). 

^(i£x  ^d£x 

To  examine  the  discontinuity  at  4>  ■ it,  we  use  both  h“2  and 
given  by  (5.31)  and  (5.34),  respectively,  since  both  J^2 
tribute  to  the  discontinuity.  Hence,  the  discontinuity  at  i|)  ■ ir  is 
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which  is  identical  to  the  discontinuity  in  (5.19c). 
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which  Is  identical  to  the  discontinuity  in  (5.20c). 

Before  closing  this  section,  it  is  worthwhile  to  recapitulate 
the  main  points  discussed.  The  discontinuities  at  “ 0,  -j,  ir,  and 
in  the  far-fleld  pattern  obtained  by  the  GTD  technique  have  been 
shown  to  be  caused  by  the  abrupt  truncation  of  the  diffracted  field  which 
is  effectively  equivalent  to  artificially  introducing  sheets  of  elec- 
tric surface  current  in  free  space.  When  these  current  sheets  radiate 
strongly  in  the  directions  ^ 0,  it,  and  they  produce  discon- 
tinuities corresponding  exactly  to  those  in  the  GTD  far-fleld  pattern. 
This  discovery  provides  us  a clue  to  further  improve  the  GTD  far-field 
pattern.  One  simply  subtracts  out  all  these  artificially  created 
"excess"  fields  from  the  GTD  far-fleld  pattern  to  get  rid  of  these 
discontinuities.  The  Improved  GTD  far-field  pattern  is  discussed  in 
the  next  section  along  with  some  comparisons  with  the  conventional 
moment  method  solution  and  results  obtained  by  a GTD-moment  method 
hybrid  technique  [4  ] and  a self-consistent  method  [24]. 

5.3  Improved  Far-Field  Pattern  and  Comparison  with  Results  Obtained 
by  Other  Approaches 

In  the  last  section,  Keller's  wedge  diffraction  coefficient  was 
used  to  construct  the  far  field  scattered  by  a rectangular  cylinder. 

It  was  shown  that  when  the  pole  singularities  contained  in  the  physical 
optics  current  with  infinite  support  on  the  seml-lnflnlte  wedge  surface 
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were  subtracted  from  the  wedge  diffraction  coefficient,  the  far  field 
remained  finite  and  varied  smoothly  across  the  geometrical  optics 
shadow  and  reflection  boundaries.  Other  chan  the  noticeable  discontinuities 
In  the  directions  In  which  the  surfaces  of  the  rectangular  cylinder  are 
oriented,  that  far-fleld  pattern  was  In  fair  agreement  with  the  con- 
ventional moment  method  solution. 

Then  It  was  found  that  these  discontinuities  were  caused  by  abrupt 
truncations  of  the  diffracted  field  to  the  regions  In  space  which  cor- 
respond to  ^iie  exterior  regions  In  the  wedge  canonical  problems.  These 
abrupt  truncations  were  demonstrated  to  be  effectively  equivalent  to 
artificial  Introductions  of  seml-lnflnlte  electric  surface  current 
sheets  In  free  space.  An  In-depth  study  of  these  current  sheets  con- 
sequently showed  that  they  radiated  strongly  and  Introduced  Identical 
discontinuities  In  those  directions  corresponding  to  the  discontinuities 
In  the  GTD  far  field.  This  discovery  provides  a clue  which  leads  to 
further  Improvement  of  the  GTD  far-fleld  pattern. 


The  elimination  of  these  discontinuities  can  be  achieved  by  sub- 
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(5.31),  (5.32),  (5.33),  and  (5.34),  correspondingly,  from  the  GTD  H 


given  In  (5.4).  The  Improved  GTD  far-fleld  pattern  so  obtained  Is 


displayed  In  Figure  5.10.  It  Is  very  Interesting  to  observe  that  not 

TT  j 3tT 


only  the  discontinuities  at  ^ 0,  y,  ir,  and  y-  disappear  completely, 

but  also  that  nulls,  which  were  not  present  In  the  GTD  far-fleld  pattern 


In  Figure  5.7,  are  showing  up  In  the  neighborhood  of  > tt  and 
3ir 


(|)  ■ -j-.  The  disappearance  of  the  discontinuities  Is  to  be  expected  since 
these  discontinuities  have  been  examined  carefully  In  the  last  section 


270® 


Figure  5.10. 


Improved  scattered  far-field  pattern  of  the  rectangular 
cylinder;  4>q  - tt/4,  a - b = IX. 


(see  (5.29),  (5.35),  (5.36),  and  (5.37))  and  were  shown  to  be  identical 
to  those  in  the  GTD  far-fleld  pattern  In  Figure  5.7. 

Finally,  results  obtained  by  using  (a)  conventional  moment  method  [17] 
with  32  unknowns,  (b)  hybrid  moment  method  and  GTD  technique  by 
Burnside  et  al.  [4  ] with  24  unknowns,  and  (c)  self-cons  Latent  method  [24] 
with  8 unknowns  are  exhibited  In  Figure  5.11.  If  one  pays  attention 
to  the  Insertion  In  Figure  5.11  showing  the  geometry  of  their  problem, 
one  will  realize  that  their  Incident  field  Is  from  a direction  exactly 
opposite  to  the  one  used  In  this  study,  hence,  their  far-fleld  pattern 
Is  opposite  to  that  In  Figure  5.10.  Note  that  In  Figure  5.11,  the  solu- 
tion obtained  by  the  self-consistent  method,  which  apparently  Includes 
all  higher-order  multiple-edge  Interactions  In  Its  formulation,  still  has 
the  discontinuities  at  ^ - 0,  90°,  180°,  and  270°.  While  the  conventional 
moment  method  solution  and  the  moment  method-GTD  hybrid  solution  deviate 
only  slightly  In  the  backscatterlng  direction,  the  present  method  solu- 
tion as  shown  In  Figure  5.10  Is  In  good  agreement  with  them. 
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Figure  5.11.  Scattered  far-fleld  patterns  obtained 

by  other  approaches.  i 
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5.4  Accuracy  Check 

In  the  previous  sections,  a method  has  been  discussed  for  obtaining 
an  Improved  scattered  far-£leld  pattern  of  a perfectly  conducting  rec- 
tangular cylinder  Illuminated  by  a plane  wave  whose  magnetic  Intensity 
vector  Is  parallel  to  the  axis  of  the  cylinder.  The  resulting  Improved 
scattered  far-fleld  pattern  has  been  compared  with  results  obtained  by 
other  approaches  In  the  literature.  In  particular,  results  obtained  by 
using  (a)  conventional  moment  method  [17]  with  32  unknowns;  (b)  hybrid 
moment  method  and  GTD  technique  by  Burnside  et  al.  [4]  with  24  unknowns, 
and  (c)  self-consistent  method  [24]  with  8 unknowns  have  been  used  for 
comparisons . 

Out  of  these  three  different  approaches,  only  the  conventional 
moment  method  has  Qonsistently  taken  Into  account  In  its  formulation 
the  boundary  condition,  requiring  the  tangential  components  of  the  total 
electric  field  be  vanished  on  the  surface  of  the  rectangular  cylinder. 

In  Burnside's  hybrid  moment  method  and  GTD  technique,  the  boundary  con- 
dition Is  not  enforced  on  the  entire  surface  of  the  scattering  object, 
as  In  the  case  of  the  conventional  moment  method.  Instead,  point  matching 
Is  applied  at  the  midpoint  of  each  of  the  pulse  current  segments  located 
near  the  diffracting  edges  of  the  wedges  of  the  rectangular  cylinder, 
and  at  some  arbitrary  points  In  the  GTD  regions  of  the  surface  of  the 
rectangular  cylinder  In  order  to  obtain  an  equal  number  of  equations  as 
the  Involved  unknowns  Including  the  current  samples  and  the  diffraction 
constants.  The  final  result  Is  somewhat  dependent  on  the  locations  of 
the  arbitrarily  chosen  additional  matching  points  In  the  GTD  region  on 
the  surface  of  the  scatterer.  Therefore,  It  takes  a Judicious  choice 


f 

i 


96 


of  the  locations  of  the  aforementioned  matching  points  and  quite  often,  In 
such  a situation,  more  matching  points  than  unknowns  are  needed  such 
that  some  averaging  schemes  can  be  applied  to  obtain  meaningful  results. 

In  the  self-consistent  method,  no  considerations  are  given  to  the  satis- 
faction of  the  boundary  condition. 

Although  the  comparisons  of  the  scattered  far-fleld  pattern  with 
those  obtained  by  the  aforementioned  methods  are  very  favorable.  Indicating 
the  solutions  are  In  good  agreement  with  the  true  solution,  an  Independent 
accuracy  check  must  be  applied  to  further  validate  the  approximate  solu- 
tion, especially  when  there  are  no  available  results  In  the  literature 
to  compared  with.  Such  an  accuracy  check  Is  often  needed  but  Is  not 
readily  available  In  the  high-frequency  asymptotic  techniques. 

In  this  section,  we  present  a method  which  allows  us  to  calculate 
the  tangential  components  of  the  scattered  electric  field  on  the  surface 
of  the  scatterer  using  the  approximate  scattered  far-fleld  pattern.  The 
accuracy  of  the  solution  can  then  be  checked  simply  by  observing  whether 
these  tangential  components  of  the  scattered  electric  field  on  the  sur- 
face of  the  scatterer  are  equal  to  the  negative  of  the  tangential  components 
of  the  Incident  electric  field  on  the  surface  of  the  scatterer.  If  the 
outcome  of  such  an  observation  Is  favorable,  then  the  solution  Is  good; 
otherwise,  further  Improvement  Is  needed. 

5.4.1  Method  of  computation 

In  classical  electromagnetic  theory.  It  Is  well-known  that  the 
scattered  magnetic  far  field  Is  related  to  the  Fourier  transform  of 
the  scattered  electric  field  on  an  aperture  In  a relatively  simple  manner. 
Consider  the  aperture  shown  In  Figure  5.12.  This  aperture  Is  a plane 
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at  y ■■  b,  containing  the  surface  of  the  rectangular  cylinder  as  Its 
central  portion.  The  scattered  far-fleld  pattern  enclosed  In  the  hemi- 
sphere corresponding  to  this  aperture  as  shown  In  Figure  5.12  Is  used 
to  compute  the  Fourier  transform  of  the  scattered  electric  field  on  the 
aperture.  In  particular,  for  the  H-wave  case  considered  In  this  chapter, 
the  Fourier  transform  of  the  tangential  component  of  the  scattered  elec- 
tric aperture  field  can  be  written  as 
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(5.38) 


where 

Ikp+lf 

Zq  Is  the  free  space  characteristic  Impedance 

T ((|i)  Is  that  portion  of  the  scattered  magnetic  far  field  pattern 
enclosed  In  the  proper  hemisphere  corresponding  to  the  aperture 
under  consideration,  and 
p and  ifi  are  the  polar  coordinates. 

The  phase  factor  In  (5.38)  Is  to  assure  that  the  aperture  at  y ■ b,  but 
not  at  y > 0,  Is  under  consideration.  The  tangential  scattered  electric 
aperture  field  In  the  space  domain  can  be  obtained  readily  by  an  Inverse 
Fourier  transform,  viz.. 


- b)  -=  F"^[E^(<^(a))]  . 


(5.39) 
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The  tangential  scattered  electric  aperture  field  in  (5.39)  is  truncated 
to  the  surface  of  the  scatterer»  i.e.,  x e [-‘a, a]  and  an  observation  can 
be  made  to  judge  whether  the  tangential  scattered  electrical  field  on 
the  surface  of  the  rectangular  cylinder  is  Indeed  equal  to  the  tangential 
Incident  electric  field. 

The  hemisphere  shown  in  Figure  5.12  Is  for  the  aperture  at  y ■ b. 
Similar  hemispheres  and  apertures  corresponding  to  the  other  surfaces 
of  the  rectangular  cylinder  can  be  established  and  the  tangential  com- 
ponents of  the  scattered  electric  field  on  these  surfaces  can  be  com- 
puted In  a similar  fashion. 

5.4.2  Results  and  comments 

In  the  previous  sections,  an  Improved  scattered  magnetic  far-fleld 
pattern  has  been  obtained  (see  Figure  5.10).  The  accuracy  of  that  pattern 
can  be  verified  by  computing  the  tangential  component  of  the  scattered 
electric  field  on  the  surface  of  the  rectangular  cylinder  using  the 
method  outlined  In  Subsection  5.4.1.  The  resulting  tangential  component 
of  the  scattered  electric  field  on  the  surface;  x t [-a, a],  y ■■  b.  Is 
shown  In  Figures  5.13a  and  5.13b.  These  curves  were  obtained  by  using  that 
portion  of  the  Improved  pattern  (Figure  5.10)  enclosed  In  the  hemisphere 
schematically  Indicated  In  Figure  5.12,  and  by  using  the  one-dlmenslonal 
FFT  for  the  Inverse  Fourier  transform  operation. 

The  magnitude  of  the  tangential  scattered  E-fleld  shown  In  Figure 
5.13a  oscillates  around  the  constant  value  266.58,  which  Is  the  magnitude 
of  the  tangential  component  of  the  incident  E-fleld.  More  precisely, 
for  normalized  H-wave  incidence,  |h^|  - lt|E^|  ■ 377,  and  the  Incident 
angle  4>q  - 45°,  therefore,  the  magnitude  of  the  tangential  component  of 
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Figure  5.13a.  Magnitude  of  the  scattered  E-field  on  the  aperture 
shown  in  Figure  5.12,  truncated  to  the  surface  of 
the  rectangular  cylinder. 
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Figure  5.13b.  Phase  of  the  scattered  E-£leld  on  the  aperture  shown 
in  Figure  5.12,  truncated  to  the  surface  of  the 
rectangular  cylinder. 
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the  Incident  E-£leld  Is  |e^  cos  45° | « 266.58.  The  phase  of  the  tangential 
scattered  E-fleld  shown  in  Figure  5.13b  Is  varying  linearly  across  the 
surface  of  the  rectangular  cylinder  and  can  be  readily  verified  to  be  of 
a phase  difference  of  tt  or  180°  from  the  phase  of  the  tangential  Incident 
E-fleld  on  the  surface.  Similar  behaviors  of  the  tangential  scattered 
E-fleld  on  other  surfaces  of  the  rectangular  cylinder  can  be  obtained 
and  are  not  repeated  here.  These  observations  confirm  that  the  Improved 
far-fleld  pattern  obtained  In  the  last  section  Is  indeed  a very  good 
approximation  to  the  true  scattered  field.  Of  course,  the  comparisons 
with  results  obtained  by  other  approaches  as  done  In  the  previous  section 
not  only  further  validate  the  approximate  solution,  but  also  demonstrate 
the  effectiveness  of  the  accuracy  checking  method  developed  In  this 
section. 

5.5  Summary 

In  this  chapter,  the  scattered  far-fleld  pattern  of  a perfectly 
conducting  rectangular  cylinder  Illuminated  by  a plane  H-wave  has  been 
obtained  by  a zeroth-order  GTD  approximation  and  the  result  Improved  by 
a straightforwardly  physical  interpretation  of  the  existence  of  the 
discontinuities  in  the  zeroth-order  GTD  far-field  pattern. 

The  Improved  scattered  far-field  pattern  has  been  compared  with 
results  obtained  by  other  different  approaches  and  the  similarities  and 
differences  between  these  results  have  been  discussed.  Generally  speaking 
the  Improved  pattern  obtained  by  the  present  approach  is  in  good  agreement 
with  the  conventional  moment-method  solution  with  32  unknowns. 

An  accuracy  checking  method  has  also  been  presented  so  that  an  in- 
dependent check  on  the  satisfaction  of  the  boundary  condition  on  the 


surface  of  the  rectangular  cylinder  can  be  performed.  The  results  of 
such  an  accuracy  check  are  quite  favorable,  which  demonstrates  that  the 
Improved  far-fleld  pattern  obtained  by  the  present  approach  Is  Indeed  a 
close  approximation  to  the  true  solution.  The  merit  of  such  an  Inde- 
pendent accuracy  checking  scheme  Is  that  the  approximate  solution  ob- 
tained can  be  validated  without  making  comparisons  with  other  methods. 
This  accuracy  check  Is  much  desired  especially  when  there  are  new  re- 
sults obtained  and  no  available  Information  In  the  literature  that  can 
be  trustworthy  to  compare  with. 
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6.  A SYNTHETIC-APERTURE-DISTRIBUTION  APPROACH  TO 
THE  HIGH-FREQUENCY  ELECTROMAGNETIC  SCATTERING 
OF  OBSTACLES  WITH  CONVEXLY  CURVED  SURFACE 


6.1  Introduction 

In  the  previous  chapters,  we  have  discussed  new  approaches  for 
solving  the  high-frequency  electromagnetic  scattering  problems  Involving 
obstacles  with  planar  surfaces  and  sharp  edges.  For  all  of  these  problems, 
the  zeroth-order  approximation  to  the  scattered  far-fleld  can  be  con- 
structed from  Keller's  wedge-dlf fraction  coefficient  (Including  the 
edge-dlf fraction  coefficient  as  a special  case)  In  a relatively  simple 
manner.  The  Improvement  of  the  zeroth-order  solution  can  be  achieved 
In  a systematic  way,  as  demonstrated  In  the  previous  chapters,  where 
the  Infinite  perfectly  conducting  thin  strip,  the  finite  perfectly  con- 
ducting rectangular  thin  plate,  and  the  perfectly  conducting  rectangular 
cylinder  have  been  used  as  Illustrative  examples. 

In  reality,  most  of  the  scattering  objects  possess  a smoothly 
curved  surface  Instead  of  planar  facets  and  sharp  edges.  A circular 
perfectly  conducting  cylinder  Is  the  simplest  geometry  of  this  category 
of  smooth-surfaced  obstacles.  The  scattering  by  a circular  cylinder  Is 
among  the  few  electromagnetic  scattering  problems  that  can  be  solved 
by  the  method  of  separation  of  variables  In  partial  differential  equa- 
tions. From  a mathematician's  point  of  view,  the  scattering  by  a cir- 
cular cylinder  Is  solved  when  the  scattered  field  Is  expressed  In  terms 
of  an  Infinite  series  of  eigenfunction  expansion  Involving  transcendental 
functions.  However,  such  an  eigenfunction  solution  usually  converges 
slowly  and  many  terms  In  the  Infinite  series  expansion  must  be  Included 
before  a "settled-down"  solution  can  be  obtained.  For  high-frequency 


scattering,  the  obstacle  is  normally  large  In  terms  of  wavelengths. 

When  It  comes  to  numerical  computations,  such  an  eigenfunction  expansion 
Is  no  more  than  a mathematical  elegance  but  Is  totally  useless  In 
practicality. 

For  many  decades,  sclentlflc-mlnded  people  have  stared  at  this 
exact  solution  which  Is  expressed  In  such  an  elegant  manner  and  yet  so 
formidable  to  get  numerical  results  out  of  It.  Before  the  advent  of 
the  latest  generation  of  large  computer  systems,  many  of  the  former 
graduate  students  have  spent  a good  deal  of  their  time  working  on  re- 
search projects  Involving  numerical  evaluations  of  the  exact  solutions 
with  the  aid  of  desk  calculators.  The  Introduction  of  the  first  large 
computer  system  was  expected  to  release  them  of  these  tedious  calcu- 
lations. However,  they  soon  found  out  that  the  evaluation  of  the  exact 
eigenfunction  expansion  solutions  on  the  computer  was  extremely  time- 
consuming.  Furthermore,  they  discovered  that  the  problem  of  numerical 
Instabilities  would  occur  when  Increasingly  higher-order  terms  were 
Included  successively  In  the  computation.  Therefore,  In  high-frequency 
scattering  where  the  higher-order  terms  must  be  Included  In  the  series 
to  obtain  a convergent  solution,  the  numerical  evaluation  of  the  exact 
solution  Is  still  hopeless. 

Even  before  their  disenchantment  with  the  large  scale  computer 
systems,  people  were  so  frustrated  by  the  numerical  evaluation  of  the 
exact  solution  that  they  started  deriving  asymptotic  formulas  for  both 
high-frequency  and  low-frequency  approximations  to  the  exact  solution. 
Most  of  the  high-frequency  asymptotic  formulas  may  be  obtained  either 
directly,  e.g.,  by  the  Luneburg-Kllne  method  [25],  or  by  asymptotic 
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evaluations  of  contour-integral  representations  of  the  exact  solution# 
Saddle-point  Integration  and  the  stationary-phase  method  are  most  often 
used  to  obtain  an  asymptotic  series  If  the  Integral  representation  of 
the  exact  solution  Is  the  starting  point  of  the  derivation. 

For  high-frequency  scattering.  In  most  of  the  cases,  the  first 
term  of  the  asymptotic-series  expansion  gives  an  extremely  good  approx- 
imation to  the  exact  solution.  However,  for  the  case  In  which  the  first 
term  does  not  give  a good  approximation  to  the  exact  solution.  Including 
the  next  higher-order  term  In  the  asymptotic  series  may  or  may  not  be 
able  to  Improve  the  approximate  solution.  This  Is  a well-known  property 
of  the  asymptotic  series, because  In  most  situations,  the  asymptotic 
series  Is  a divergent  Instead  of  a convergent  series. 

To  date,  a method  for  systematically  Improving  the  asymptotic 
solution  Is  still  the  most  desirable.  In  spite  of  all  the  efforts 
that  have  been  concentrated  on  this  solution  In  the  two  decades  since 
Keller's  first  conception  of  the  diffraction  coefficient  approach,  no 
systematic  way  exists  to  Improve  the  hlgh-frequency-asymptotlc  solution 
when  necessary.  Many  theories,  which  are  Improved  versions  of  the 
original  Keller's  geometrical  theory  of  diffraction  (GTD)  [26]  have  been 
proposed,  among  these  theories  are  the  uniform  theory  of  diffraction 
(UTD)  [8],  the  uniform  asymptotic  theory  (UAT)  [9],  and  the  spectral 
theory  of  diffraction  (STD)  [1(9 . Nevertheless,  the  computational  ef- 
ficiency of  these  theories  when  applied  to  practical  situations  still 
remains  to  be  seen. 

Two  of  the  state-of-the-art  methods  for  solving  the  high-frequency 
scattering  problems  Involving  perfectly  conducting  convex  cylindrical 
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structures  have  been  presented,  one  has  been  published  In  the  literature 
[4j,  while  the  other  has  been  presented  at  the  USNC/URSI  1975  annual 
meeting  [5  ].  The  first  method  uses  a combination  of  the  geometrical 
theory  of  diffraction  and  the  moment  method  to  solve  a circular  cylinder 
scattering  problem,  wherein  the  cylindrical  surface  is  divided  into 
three  regions,  namely,  the  physical  optics  region  in  the  well-illuminated 
side,  the  GTD  region  in  the  deep  shadow,  and  the  pulse  region  in  the 
transition  region  between  the  lit  and  the  dark  sides.  The  surface 
current  is  first  obtained,  then  the  far-fleld  can  presumably  (Burnside 
et  dl.  did  not  show  the  far-field  pattern)  be  obtained  by  a numerical 
double-integration  of  the  surface  current.  There  are  three  weaknesses 
in  Burnside's  method: 

(1)  The  composite  surface  current  computed  by  using  different  tech- 
niques in  different  regions  on  the  cylindrical  surface  as  in- 
dicated above  is  not  smoothly  connected  at  the  boundaries  of 

each  region.  These  discontinuities  are  clearly  shown  in  Figure  6.20; 

(2)  The  far-fleld  pattern  must  be  computed  by  a numerical  double- 
integration  of  the  surface  current.  Such  a computation  may  be 
time  consuming;  and 

(3)  Furthermore,  the  discontinuities  at  the  high  end  in  the  surface 
current  (see  Figure  6.20)  may  introduce  extraneous  high-level  side- 
lobes  in  the  far-fleld  pattern  (see  Section  6.4  for  a demonstration). 
The  second  method  used  an  N-slded  polygon  to  model  the  circular 

cylinder.  Wang  and  Richmond's  analysis  is  based  on  the  self-consistent 
geometrical  theory  of  diffraction  [24].  Two  cylindrical  waves  with 
unknown  amplitudes  are  assumed  to  travel  in  opposite  directions 
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(traverse  to  the  axis  of  the  cylinder)  on  each  facet  of  the  polygon. 

The  boundary  conditions  for  the  corners  are  applied  to  set  up  a matrix 
equation  for  2N  unknowns  (the  amplitudes  associated  with  the  traveling 
cylindrical  waves).  Grout's  method  Is  used  to  solve  the  matrix  equation. 
After  the  amplitudes  for  the  traveling  waves  have  been  determined, 
the  far  field  Is  obtained  by  the  self-consistent  method.  There  are 
three  weaknesses  In  Wang  and  Richmond's  method: 

(1)  Sharp  edges  are  Introduced  artificially  In  modeling  the  smooth 
circular  cylindrical  surface  by  a polygonal  cylinder.  In  other 
words,  the  nature  of  the  problem  has  been  changed  substantially; 

(2)  In  order  to  accurately  model  an  electrically  large  circular  cylinder. 
It  Is  conceivable  that  N,  the  number  of  polygonal  facets  used, 

must  be  large;  otherwise,  the  problem  of  a polygonal  Instead  of  a 
circular  cylinder  Is  being  solved.  Hence,  the  number  of  unknowns, 

2N,  In  the  matrix  equation  can  be  large  and  the  computation  be- 
comes, again,  time  consuming;  and 

(3)  In  addition,  the  self-consistent  method,  as  pointed  out  In  the 
last  chapter,  may  not  be  able  to  get  rid  of  the  extraneous  vari- 
ations In  the  far-fleld  pattern  In  the  directions  coincident  with 
the  orientations  of  the  Individual  facets  of  the  modeling  poly- 
gonal cylinder. 

It  Is  also  well-known  In  the  high-frequency  scattering  that  asymp- 
totic evaluation  of  the  Integral  representation  of  the  scattered  field 
results  in  an  Infinite  series  whose  terms  can  be  Interpreted  as  highly 
attenuated,  surface-mode,  traveling-waves  around  the  surface  of  the 
obstacle.  The  diffraction  coefficients  and  attenuation  constants  In 


the  Keller's  geometrical  theory  of  diffraction  for  smoothly  curved 
surfaces  can  be  identified  from  the  terms  in  the  series.  Such  traveling 
waves  propagate  around  the  obstacle  at  a slightly  slower  speed  than 
that  of  light  in  free  space,  and  continuously  shed  energy  in  directions 
tangential  to  the  surface  as  they  travel  along;  hence,  they  are  highly 
attenuated.  These  waves  are  known  as  "creeping  waves."  It  was  once  of 
great  hope  that  the  creeping-wave  contributions  would  complete  the  pic- 
ture of  the  geometrical-optics  solution  of  scattering  problems  by  pro- 
viding a means  of  calculating  the  diffracted  field  in  the  deep-shadow 
region.  Conceptually,  the  creeping  wave  theory  is  attractive;  but, 
actually,  due  to  the  hlgh-attenuatlon  property  of  the  creeping  waves, 
'those  waves  that  "crept"  around  the  obstacle  more  than  once  contribute 
little  to  the  total  diffracted  field.  Therefore,  if  the  first-order 
creeping  wave  failed  to  Improve  the  geometrical-optics  solution  to  a 
certain  extent,  including  higher-order  creeping  waves  seldom  Improves 
the  solution  further.  Besides,  creeping  waves  are  too  cony>llcated  to 
compute,  hence,  are  not  suited  for  programming  on  computers. 

All  of  the  above  discussions  point  to  the  fact  that  a conceptually 
simple,  computationally  efficient  method  of  solving  high-frequency 
electromagnetic  scattering  problems  involving  curved-surface  scatterers 
is  still  much  desired. 

In  this  chapter,  we  Introduce  a new  approach  to  the  high-frequency 
electromagnetic  scattering  of  obstacles  with  a convexly  curved  surface. 
The  idea  is  to  transfer  information  characterizing  the  curved  surface 
onto  a planar  aperture  where  maneuvers  for  improving  the  geometrical- 
optics  solution  are  carried  out.  In  conventlonp'i.  strategies,  maneuvers 
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for  Improvement  are  carried  out  either  directly  on  the  surface  current 
or  directly  on  the  far  field.  The  Introduction  of  a planar  aperture  In 
the  present  approach  Is  a marked  difference  from  all  of  the  conventional 
methods.  A salient  feature  of  the  method  Is  that  the  far  field  Is 
obtained  by  Fourier  transforming  the  aperture-field  distribution;  hence, 
It  can  handle  problems  Involving  an  n-dlmenslonal  obstacle  by  an 
(n  - 1) -dimensional  fast  Fourier  transform  (FFT) , where  n 2,3. 
Therefore,  the  unwieldy  three-dimensional  FFT  Is  circumvented,  hence, 
the  method  Is  numerically  efficient.  Another  Important  feature  of  the 
method  Is  that  the  aperture-field  distribution  Is  slowly  varying,  l.e., 
not  rapidly  oscillatory.  In  magnitude  and  In  phase.  Hence,  maneuvers 
for  Improvement  can  be  achieved  In  a relatively  straightforward  manner. 
Moreover,  the  method  Is  conceptually  simpler  than  other  conventional 
methods  in  the  literature,  because  maneuvers  for  improvement  are  car- 
ried out  on  a plane  rather  than  on  a curved  surface. 

The  best  way  to  explain  the  method  Is  to  apply  it  to  solve  the 
scattering  problem  of  a plane  wave  by  a perfectly  conducting  circular 
cylinder,  because  the  principle  Involved  can  be  Illustrated  without 
the  unnecessary  complexities. 

To  make  this  chapter  as  self-contained  as  possible,  we  briefly 
derive  the  exact  solution  In  Section  6.2  and  the  geometrical-optics 
solution  in  Section  6.3.  Synthesizing  the  approximate  aperture-field 
distribution  Is  demonstrated  In  Section  6.4.  Computation  of  surface 
current  Is  carried  out  In  Section  6.5.  An  accuracy  check  of  the  ap- 
proximate solution  Is  discussed  in  Section  6.6.  A summary  for  the 
chapter  can  be  found  In  Section  6.7. 
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6.2  Exact  Solution 

In  this  section,  we  briefly  discuss  the  exact  solution  to  the  prob- 
lem of  scattering  of  a plane  wave  by  a perfectly  conducting  circular 
cylinder.  The  geometry  of  the  problem  Is  shown  In  Figure  6.1.  A per- 
fectly conducting  circular  cylinder  located  at  the  origin  of  the  co- 
ordinate system  Is  Illuminated  by  a uniform  plane  wave  whose  electric 
Intensity  vector  Is  parallel  to  the  axis  of  the  cylinder. 

Using  the  method  of  separation  of  variables  [27],  the  scattered 
field  can  be  written  as 

E®(p,<j>)  - I (1)"  a^H^^^(kp)  e^”^  (6.1) 

n— 00 

where  H^^^(kp)  Is  the  Hankel  function  of  the  first  kind,  and  the  a ’s 
n n 

are  unknown  coefficients.  The  time  convention  Is  understood. 

To  determine  the  we  evaluate  (6.1)  at  p ••  a,  the  radius  of  the 

cylinder,  l.e.,  at  the  surface  of  the  cylinder  and  then  apply  the 
boundary  condition.  To  do  this,  we  first  expand  the  Incident  electric 
field  Into  a series,  namely, 

00 

E^Cp.-fr)  - I (1)"  J^(kp)  e^"*^  (6.2) 

n— 00 

where  J^(kp)  Is  the  Bessel  function  of  the  first  kind.  The  boundary 
condition  requires  that  the  total  tangential  electric  field  vanish  on 
the  surface  of  a perfectly  conducting  obstacle.  Therefore, 

E®(a,(|.)  + E^(a,(^)  - 0 . 

2 2 


(6.3) 


113 


Using  (6.1)  and  (6.2)  In  (6.3),  we  obtain 


a 


-J  (ka) 
_ n 

n 


(6.4) 


Substituting  (6.4)  in  (6.1),  we  have 


“ J (ka)  » . . 

E®(p,<|.)  - - I (i)"nif) 

* n— » H^^^ka)  " 

n 


(6.5) 


The  far-fleld  pattern  of  the  scattered  field  is  of  Interest.  At 
large  distances  from  the  cylinder,  we  can  use  the  asymptotic  formula 
for  namely. 


to*.  Ji£  * 


Ikp 


(6.6) 


and  (6.5)  becomes 


r-T“  ir  “ J (ka)  . . 

' ^ ikp  „ ^ n 


"z<P'*>^7.kp®  ^ ® ^L«H<^>(ka) 


(6.7) 


The  magnitude  of  the  ratio  of  the  scattered  field  to  the  incident  field 
is  therefore 


J irkp 


n«.®  (to) 

n 


(6.8) 


This  is  the  scattered  far-fleld  pattern.  For  ka  ■ 1,  5,  10,  and  6ir, 
the  far-fleld  patterns  are  shown  in  Figures  6.2a,  6.2b,  6.2c,  and  6. 2d. 
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Experience  shows  that  the  number  of  terms  one  should  take  In  (6.8) 
to  compute  the  far-fleld  pattern  Is  of  the  order  of  ka,  which  means 
that  the  convergence  of  the  series  is  extremely  slow  for  an  electrically 
large  cylinder  and  the  computation  becomes  time-consuming.  For  high- 
frequency  scattering,  one  must  resort  to  approximation  solutions  to 
obtain  meaningful  numerical  results. 

The  surface  current  on  the  cylinder  Is  also  of  Interest.  It  can 
be  obtained  from  Maxwell's  equations,  l.e.. 


3E 


p*a 


-iojy  3p 


p-a 


(6.9) 


i s 

where  E « E + E . Using  (6.2)  and  (6.5),  and  simplifying  the  result 
z z z 

by  the  Wronsklan  of  Bessel's  equation,  we  obtain 


,n 


trkaZr,  ^ u(l) 
0 n«-“>  H 

n 


Iniji 

(ka) 


(6.10) 


where  Zq  is  the  free-space  characteristic  Impedance.  The  surface  cur- 
rent for  the  case  of  a ■ 3X  is  shown  in  Figures  6.3a  and  6.3b. 

6.3  Geometrical  Optics  Solution 
6.3.1  Fields  of  a ray 

A point  source  radiating  In  all  directions  Is  Indicated  In  Figure  6.4. 
The  energy  between  a set  of  such  rays  Is  assumed  u>>  remain  constant. 

Four  of  these  rays,  as  shown  In  Figure  6.5,  constitute  a tube  of  rays. 

The  energy  passing  through  the  cross-sectional  areas  and  A2  must 
be  a constant  under  the  assumption  of  no  fringing.  A linear  relation 


SURFACE  CURRENT  PHASE 


Figure  6.3b.  Phase  of  the  total  Induced  surface  current  on  a circular  cylinder 
with  radius  a = 3X. 
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exists  between  the  length  of  a side  and  Its  position  p or  (p  + 1) 
shown  In  Figure  6.5.  Thus,  the  ratio  of  the  areas  Is 


2 

£_ 


(P  + 


(6.11) 


A field  quantity  at  the  position  (p  + 1)  may  be  written  In  terms 
of  Its  value  at  the  position  p . The  field  quantity  considered  here  Is 
a scalar.  A vector  field  can  be  resolved  Into  two  component  fields, 
one  polarized  In  and  the  other  normal  to  the  plane  of  Incidence,  and 
eAch  conq>onent  field  treated  as  a scalar  field.  It  should  be  empha- 
sized that  this  Is  exactly  the  physical  scattering  process  that  occurs 
at  a boundary  between  two  different  media.  The  field  quantity  Is  then 


p + 1 


(6.12) 


where  U represents  the  scalar  E-  or  H-fleld.  , 

One  may  also  consider  a more  general  case  Involving  an  astigmatic 
ray  tube  shown  In  Figure  6.6.  In  this  case,  the  scalar  field  quantity 
Is  given  by 


“2  ■ 


l)(p2  + il) 


(6.13) 


The  electrical  phase  of  the  ray  tube  Is  given  by  the  spatial  phase 
delay  factor  e^^^,  when  the  time  convention  e Is  understood.  This, 


of  course,  assumes  that  the  phase  reference  position  coincides  with 
the  amplitude  reference.  Now  one  may  write  an  expression  for  the  fields 
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of  this  astigmatic  tube  of  rays  In  the  form 


^2  " V >/(p^  + 0(p2  + 


lk£ 


1) 


(6.14) 


There  are  other  phase  changes  that  must  be  carefully  preserved. 


Observe  that  the  quantity  under  the  radical  becomes  negative  as  I be- 

,0  IT 


comes  less  than  -p^.  This  Introduces  a phase  Jump  of  -90  or  - radians. 
This  point  Is  called  a caustic.  The  fields  at  the  caustic  point  cannot 
be  determined  by  geometrical  optics  but  the  correct  fields  can  be  found 


on  either  side  of  the  caustic.  If  1 becomes  less  than  p^,  then  the 


phase  shift  of  180  or  ir  radians  occurs. 

The  general  geometrical  optics  expression  for  the  field  may  be 
written  In  the  form 


U(l)  - UqO  FU)  e 


(6.15) 


where  Uq  Is  a reference  amplitude  at  £ - 0; 


(|)q  Is  the  reference  phase  at  £ » 0; 


I ^1^2 

F(£)  - J-^ — +"i)  (p — +~[)  spatial  attenuation  factor; 


ik£ 


and  e Is  the  spatial  phase  delay  factor.  In  other  words, 


ua)  - Iteferenc.  field)  ettenuetlo!:]  fepatlal  phase 

factor  delay  factor 


(6.16) 
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6.3.2  Reflection  at  a curved  surface  boundary  separating  two 
different  media 

In  this  section  we  discuss  a simple  technique  for  obtaining  the 
reflected  field  from  a curved  surface  which  Is  a boundary  between  two 
different  media.  The  cylindrical  or  two-dimensional  geometry  Is  con- 
sidered simply  to  avoid  any  unnecessary  cosq>lexltles.  Other  surfaces 
can  be  treated  In  the  same  manner. 

The  source  Is  a line  source  as  shown  In  Figures  6.7a  and  6.7b.  A plane 
wave  is  obtained  simply  by  allowing  to  go  to  infinity.  The  reflected 
rays  now  appear  to  be  coming  from  a virtual  source  located  a distance 
Inside  the  surface.  All  factors  In  (6.16)  should  be  written  In 
terms  of  distance  from  the  element  of  the  surface  as  Illustrated  by 

The  quantity  r^  is  the  radius  of  curvature  of  the  surface. 

Now  we  proceed  to  determine  the  reflected  fields.  From  (6.15), 
the  reflected  field  Is  obtained  and  takes  the  form 


U 


r 


iklj^ 

e 


(6.17) 


where  R Is  the  reflection  coefficient  and  Uq  Is  the  reference  field  at 
the  surface. 

The  only  remaining  task  is  to  evaluate  in  (6.17).  From  trigo- 
nometry, the  angle  (see  Figure  6.7a)  Is  expressed  as  follows: 

= TT  - a - (ir  - 6q) 

= 6q  - a . (6.18) 


Similarly, 
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(a) 

Figure  6.7a.  Coordinate  system. 


(b) 


Figure  6.7b.  Reflection  from  a curved  surface. 
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(6.19) 


The  element  of  arc-length  Is  equal  to  r^(Aa)  (see  Figure  6.7b)  and  It 
is  also  equal  to  (ACj^)  Aq/cos  0q.  Therefore, 

(45  jl,  (48,  - 4.)  to 

“ (JOQ  0^  “ cog  9^ 


(6.20) 


where  the  last  equality  Is  obtained  by  means  of  (6.18)  and  A denotes 
the  differential  In  the  usual  sense.  By  the  same  token, 

(6.21) 


r^(Ao) 


(AC2> 

cos  6n 


(A6q  + Aa) 
cos  6/v 


where  the  last  equality  Is  obtained  with  the  aid  of  (6.19).  Eliminating 
Aa  from  (6.20)  and  (6.21)  yields 


Pi 


2 

r^  cos  0Q 


(6.22) 


Substituting  (6.22)  In  (6.17)  completes  the  solution  for  the  reflected 
field  from  the  curved  surface. 

Thus,  the  expression  for  the  field  reflected  by  a two-dimensional 
curved  surface  has  been  obtained  In  a relatively  simple  manner.  The 
validity  of  this  expression  has  been  checked  by  comparing  solutions  of 
this  type  with  exact  solutions  for  spheres  and  cylinders  [28].  For 
further  discussions  on  the  subject  of  electromagnetic  reflection  from  a 
perfectly  conducting  surface,  the  Interested  reader  Is  referred  to  the 
literature  [29]. 
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6.3.3  Scattering  of  a plane  wave  by  a circular  cylinder 
In  this  section  we  apply  the  general  formulas  derived  In  the  last 
section  for  reflection  at  an  arbitrarily  convexly-curved  boundary-surface 
to  the  special  case  of  a circular  cylinder.  The  geometry  of  the  scatter- 
ing problem  Is  shown  In  Figure  6.8.  A plane  wave  polarized  parallel  to  the 
axis  of  the  circular  cylinder  located  at  the  origin  of  the  coordinate 
system  Is  Incident  from  the  negative  x-dlrectlon,  viz., 


=1,  ,v  ..  lkpcos(|> 
E (p,(ji)  - ze 


(6.23) 


-lujt 


where  the  time  convention  e Is  understood. 

To  obtain  the  geometrical  optics  reflected  field,  we  use  (6.17). 
Since  this  Is  a perfectly  conducting  circular  cylinder,  the  reflection 
coefficient,  R,  Is  -1.  The  reference  field,  Uq,  Is  the  Incident  field 
at  the  point  of  reflection,  P^,  on  the  circular  cylinder.  From  (6.23), 


Uq  “ E^(p  - a,<p  - O 


IkacosC 


(6.24) 


The  spatial  attenuation  factor  Is  determined  when  the  caustic  position, 
p^.  Is  known.  To  find  p^,  we  use  (6.22).  For  plane  wave  Incidence, 

-*■  «>.  The  radius  of  curvature  In  this  case  Is  simply  a constant,  l.e., 
the  radius  of  the  circular  cylinder,  a.  Hence,  the  caustic  position 
Is  given  by 


a cos  6, 


(6.25) 
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Note  that  the  angle  of  reflection,  6^,  in  (6.25)  Is  a function  of 
l.e.,  It  depends  on  the  location  of  the  point  of  reflection  on  the 
cylindrical  surface.  The  quantity  In  (6.17)  Is  the  distance  between 
the  point  of  reflection,  P^(a,C),  and  the  point  of  observation,  r(p,i|)). 

This  distance  Is  denoted  by  £ In  Figure  6.8.  To  obtain  an  expression 
for  the  scattered  field  at  an  observation  point,  r(p,4>),  we  must  first 
find  a relation  between  the  parameters  (6q,£,C)  and  the  coordinates 
(p,(^).  However,  there  Is  no  simple  relation  that  can  be  found  to  ex- 
press 0Q,  £,  and  C explicitly  In  terms  of  p and  4>.  Therefore,  In 
practice,  one  has  to  resort  to  numerical  methods.  ' There  are  basically 
two  routes  that  one  can  follow,  namely,  either  searching  for  the  point 
of  reflection  when  a point  of  observation  Is  given,  or  Interpolating 
for  the  field  when  a point  of  reflection  Is  chosen. 

Nevertheless,  a simple  closed-form  expression  can  be  found  for  the 
far  field  (p  -^  ») . In  particular,  we  use  the  following  approximations 
when  p >>  a > X; 

p - £ In  the  amplitude  (6.26a) 

p <^  £ + a cos  In  the  phase  (6.26b) 

®0  " 1^  2 0 £ (j>  £ 2tt  (6.26c)  . 

and 

C = for  0 < $ < 2w  . (6.26d) 

Using  (6.24),  (6.25),  and  (6.26)  In  (6.17),  and  after  some  algebraic 
manipulations,  we  obtain 
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where 


0 _<()»£  2ir 
p >>  a > X 

and  the  time  convention  e la  understood.  The  scattered  far-£leld 
pattern  obtained  from  (6.27)  for  the  case  of  a - 3X  Is  shown  In  Figure 
6.9.  For  further  discussions  on  the  geometrical  optics  reflected  fields, 
the  Interested  reader  Is  referred  to  the  literature  [30]. 

6.4  Synthesizing  the  Approximate  Aperture  Field  Distribution 

In  the  last  section  we  discussed  the  geometrical  optics  solution 
to  the  problem  of  scattering  of  a plane  wave  by  a perfectly  conducting 
circular  cylinder.  The  geometrical  optics  scattered  far-fleld  pattern 
Is  In  excellent  agreement  with  the  exact  (eigenfunction  series  expansion) 
solution  In  the  hemisphere  containing  the  lit  surface  of  the  circular 
cylinder.  However,  the  geometrical  optics  solution  falls  completely  In 
the  forward-scattering  direction  because  It  predicts  a zero  field  In- 
stead of  the  maximum  of  the  main  beam  In  that  direction.  This  result 
does  not  come  as  a surprise  since  In  a geometrical  optics  solution  only 
the  Illuminated  portion  of  the  surface  of  an  obstacle  Is  taken  Into  ac- 
count, whereas  the  shadowed  portion  of  the  surface  of  the  obstacle  Is 
completely  Ignored  In  the  solution.  As  a matter  of  fact,  the  shadowed 
side  of  the  scatterer  could  have  been  of  any  arbitrary  shape  and  the 
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geometrical  optics  solution  would  have  still  predicted  the  same  scattered 
£ar-£leld  pattern  as  long  as  the  Illuminated  portion  o£  the  sur£ace 
were  cylindrical. 

In  spite  o£  the  serious  discrepancy  In  the  £orward-scatterlng  di- 
rection, the  geometrical  optics  £ar-£leld  solution  Is,  nonetheless.  In 
good  agreement  with  the  exact  solution  £or  a wide  range  o£  angles  ex- 
tending £rom  the  lit  hemisphere  deep  Into  the  dark  hemisphere.  Only 
within  a narrow  range  of  angles  In  the  Immediate  neighborhood  o£  the 
£orward-scatterlng  direction  does  the  geometrical  optics  £ar-£leld  start 
to  deviate  £rom  the  exact  £ar-£leld  pattern.  This  observation  leads  to 
the  conclusion  that  the  geometrical  optics  solution  xs  a very  good  ap- 
proximation to  the  exact  solution  except  In  the  neighborhood  o£  the 
£orward-scatterlng  direction. 

It  Is  our  conjecture  here  that  the  main  beam  In  the  £orward- 
scatterlng  direction  may  be  obtained  £rom  an  aperture-£leld  distribution. 
1£  an  ln£lnlte,  planar  aperture  cutting  the  entire  space  Into  two  hall- 
spaces  Is  erected  on  the  £ar-slde  o£  the  cylinder  when  we  view  It  £rom 
the  source  side,  we  may  stop  all  the  geometrical  optics  rays  at  this 
aperture  and  calculate  the  aperture-£leld  distribution.  For  observa- 
tion points  beyond  this  aperture,  we  simply  use  Huygens'  principle  to 
£lnd  the  £lelds.  Theoretically,  this  procedure  o£  stopping  the  rays 
at  the  aperture  and  then  letting  the  aperture-£leld  propagate  beyond 
It  should  not  be  dl££erent  £rom  the  original  ray  £lelds  which  propagate 
continuously  undisturbed.  Hence,  we  should  obtain  nothing  more  than 
the  original  geometrical  optics  £ar-£leld  pattern.  However,  1£  the 
aperture  Is  In  juxtaposition  with  the  cylinder,  then  that  portion  o£ 


the  aperture  field  In  the  vicinity  of  the  cylinder  must  be  completely 
different  from  the  ray  fields.  As  a matter  of  fact,  the  geometrical 
optics  ray  fields  are  no  longer  meaningful  In  that  region.  Consequently, 
an  option  Is  opened  here  for  modifying  the  aperture-field  distribution 
so  that  an  Improved  far-fleld  pattern  can  be  obtained.  Methods  of 
synthesizing  the  aperture-field  distribution  are  discussed  In  the  follow- 
ing subsections. 

The  format  of  presentation  of  the  methods  is  as  follows.  In  each 
of  the  subsections,  a method  of  synthesizing  the  rture-fleld  distri- 
bution Is  Illustrated  and  the  resulting  far-fleld  pattern  displayed. 

Each  subsection  Is  concluded  by  a discussion  of  the  method  and  a hint 
for  a better  synthesis.  These  methods  for  synthesizing  the  aperture- 
field  distribution  are  arranged  In  the  order  of  Increasing  degree  of 
sophistication  and  accuracy.  The  first  two  methods  presented  below 
give  totally  unsatisfactory  far-fleld  patterns;  however,  they  serve 
as  good  examples  to  Indicate  some  of  the  pitfalls  that  one  may  run  Into, 
hence,  their  Inclusion  here  Is  well  justified. 

6.4.1  Method  I 

All  scientific  progress  starts  from  observations.  It  Is  no  ex- 
ception In  the  development  of  this  method.  The  exact  scattered  electric 
aperture  field  In  the  shadow  region  cast  by  the  geometrical-optics  rays 
Is  depicted  In  Figure  6.10.  Note  that  the  width  of  this  shadow  Is 
simply  the  diameter  of  the  circular  cylinder.  In  this  particular  ex- 
ample, the  radius  of  the  circular  cylinder  Is  thrice  the  wavelength  of 
the  Incident  plane  wave  which  Is  polarized  In  the  z-dlrectlon  and  has 


a magnitude  of  1.  Due  to  symmetry » only  half  of  the  shadow-region  scat- 
tered aperture-field  Is  plotted. 

The  magnitude  curve  has  a constant  value  1 at  the  center  of  the 
shadow  and  hardly  varies  until  It  Is  about  two  wavelengths  away  from  the 
center,  at  which  point  the  magnitude  starts  to  decrease.  The  overshoot 
In  the  magnitude  curve  Is  highly  exaggerated  because  the  scale  runs  from 
0.7  to  1.1  only.  The  value  of  the  magnitude  of  the  exact  scattered 
aperture  field  at  the  edge  of  the  geometrical-optics  shadow  Is  approx- 
imately 0.73.  Note  that  the  phase  variation  throughout  the  entire 
shadow  region  Is  nothing  more  than  attributing  a negative  sign  to  the 
magnitude.  This  fact  does  not  come  as  a surprise,  because.  In  the 
deep-shadow  region,  the  total  electric  field  should  be  zero;  hence,  the 
scattered  E-fleld  In  that  region  should  be  just  the  negative  of  the 
Incident  E-fleld. 

Having  observed  this  fact  about  the  scattered  electric  aperture 
field,  a first  attempt  to  Improve  the  geometrical-optics  far-fleld  Is 
naturally  to  assume  an  aperture  distribution  which  Is  the  negative  of  the 
Incident  E-fleld  In  the  geometrical-optics  shadow  region  and  Is  zero 
outside  this  shadow  region.  Alternatively,  the  scattered  electric 
aperture  field  outside  the  shadow  region  could  have  been  assumed  the 
value  of  the  geometrical-optics  scattered  E-fleld.  However,  the  process 
of  stopping  the  geometrical-optics  rays  at  the  aperture  and  then  letting 
them  radiate  again  to  generate  the  far-fleld  Is  not  much  different  from 
the  far-fleld  obtained  directly  from  the  geometrical-optics  rays  prop- 
agating to  the  far-fleld  zone  undisturbed.  The  only  significant  dif- 
ference Is  an  additional  contribution  generated  by  the  shadow-region 
aperture-field  to  the  far-fleld.  The  scattered  far-fleld  pattern  Is 
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obtained  by  superimposing  the  geometrical-optics  far-fleld  and  the 
additional  far-fleld  computed  by  Fourier  transforming  the  assumed 
scattered  aperture  E-fleld  in  the  shadow  region. 

The  far-fleld  pattern  obtained  by  this  method  is  shown  as  a dotted 
curve  in  Figure  6.11.  Recall  that  the  Fourier  transform  of  the  rectangular 
pulse  aperture  distribution  results  In  a sin  x/x,  or  the  sine  function, 
type  of  far-fleld  pattern.  Therefore,  the  main  beam  In  the  forward- 
scattering direction  shows  up.  A comparison  of  Figure  6.2,  Figure  6.9, 
and  Figure  6.11  manifests  the  remarkable  ln^irovement  In  the  geometrical- 
optics  far-fleld  pattern  by  the  additional  far-fleld  generated  by  the 
deep-shadow  aperture-distribution  which  Is  precisely  what  was  Ignored 
In  the  geometrical-optics  solution.  The  deep  nulls  In  the  neighborhood 
of  the  forward-scattering  direction  are  due  to  the  phase  Interferences 
between  the  geometrical-optics  solution  and  the  Fourier  transform  of  the 
aperture-distribution.  Moreover,  the  relatively  pronounced  oscillations 
in  the  far-fleld  pattern  are  caused  by  the  discontinuities  in  the  aper- 
ture-field distribution.  In  the  next  subsection  we  discuss  another  method 
of  synthesizing  the  aperture-field  distribution  so  that  the  extraneous  oscil- 
lations can  be  suppressed. 

6.4.2  Method  II 

The  method  discussed  In  the  last  subsection  substantially  Improved 
the  geometrical-optics  solution  of  the  scattered  far-fleld  pattern. 

However,  extraneous  oscillations  occurred  in  the  far-fleld  pattern  and 
further  Improvement  is  necessary.  It  is  well-known  in  the  field  of 
antenna  design  that  a tapered  aperture  field  produces  a low  side-lobe 
radiation-pattern.  Therefore,  the  second  attempt  to  improve  the 
geometric-optics  solution  is  to  synthesize  a smoothly  tapered  aperture- 
distribution  of  the  scattered  E-fleld.  To  this  end,  consider  a 
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sinusoidal  roll-off  function  of  the  form 

f (y)  - - f + I cos  7T  (6.28) 

where  A and  B are  parameters.  Note  that  the  attractive  property  of 
such  a sinusoidal  roll-off  function  Is  that  It  serves  as  a very  good 
Interpolating  curve  between  the  constant  values  -1  and  zero.  Not  only 
the  values,  but  also  the  slopes  of  the  function  being  Interpolated  are 
matched  at  the  two  end  points  of  the  Interpolation  Interval  designated  by 
[A,B].  In  a sense,  such  a curve  possesses  the  properties  of  a spline. 

The  parameters  A and  B are  chosen  to  correspond  to  the  width  of 
the  penumbra-region,  or  better-known  as  the  "transition-region"  on  the 
aperture.  To  date,  no  reliable  techniques  are  available  In  the  liter- 
ature to  determine  the  width  of  the  transition-region  to  a considerable 
degree  of  accuracy.  In  the  literature.  It  has  always  been  assumed  that 
the  transition-region  Is  located  symmetrically  about  the  boundary  of 
the  geometrical  shadow.  However,  such  an  assumption  Is  not  very  accurate  be- 
cause from  the  exact  aperture -fie Id  distribution,  the  transition-region 
Is  observed  to  be  asyrmetvia  about  the  boundary  of  the  geometrical- 
optics  shadow  and  It  extends  much  farther  Into  the  Illuminated  region 
than  the  dark  region.  This  point  Is  discussed  further  In  the  next  sub- 
section. 

Nevertheless,  It  Is  educational  to  compute  the  far-fleld  using 
the  aperture-field  distribution  synthesized  from  the  tapering  sinu- 
soidal roll-off  functions  In  the  penumbra-region  and  the  negative  of 
the  Incident  field  In  the  deep-shadow  region.  In  this  example,  the 
radius  of  the  cylinder  Is  3X;  the  closest  approximation  to  the  width 
of  the  transition  region  Is  IX  on  either  side  of  the  boundary  of  the 
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geometrical-optics  shadow,  which  means  that  a 2X-transltlon  region  Is 
assumed. 

The  far-fleld  pattern  obtained  by  this  method  Is  shown  as  a solid 
curve  In  Figure  6.11.  As  expected,  the  extraneous  oscillations  present 
In  the  far-fleld  pattern  obtained  by  Method  I are  suppressed.  However, 
there  are  still  quite  an  amount  of  discrepancies  from  the  exact  solu- 
tion. The  reason  for  these  discrepancies  can  be  traced  to  the  Ignorance 
of  the  phase  Information  In  the  course  of  solution.  So  far,  the  phase 
has  been  assumed  to  take  on  the  value  “ir  In  the  entire  aperture-field 
distribution.  Actually,  this  Is  not  the  case  In  the  transition  region. 

The  phase  Information  must  be  Included  In  the  process  of  solution; 
otherwise,  the  far-fleld  pattern  cannot  be  Improved  further  by  simply 
manipulating  the  magnitude.  In  the  next  subsection  we  discuss  another 
method  of  synthesizing  the  aperture  distribution  so  that  both  the  mag- 
nitude and  phase  variations  are  taken  Into  account  In  the  process  of 
solution;  hence,  further  Improvements  on  the  far-fleld  pattern  can  be 
achieved. 


6.4.3  Method  III 

In  the  last  two  subsections,  methods  of  synthesizing  the  scattered 
electric  aperture-field  distribution  in  the  geometrical-optics  deep- 
shadow  region  and  the  transition  or  the  penumbra  region  have  been  dis- 
cussed. Although  the  manipulations  for  magnitude  only  do  give  a sub- 
stantial Improvement  over  the  geometrical-optics  far  field,  the  failure 
to  account  for  the  phase  variation  In  the  transition  region  prevents 
us  from  further  Improving  the  far-fleld  pattern.  Recognizing  the  Im- 
portance of  the  phase-variation,  we  reexamine  the  basic  principles 


underlying  the  synthesis  of  the  aperture  distribution.  It  Is  found 
that  If  the  procedures  In  the  following  paragraph  are  followed  to 
synthesize  the  aperture  distribution*  the  proper  phase  Information  can 
be  Included  In  the  solution. 

As  alluded  to  earlier  In  Subsection  6.4.1  when  we  discussed  Method 
1,  the  scattered  electric  aperture  field  In  the  well-lllumlnated  region 
may  alternatively  be  approximated  by  the  geometrical-optics  scattered 
far-fleld  evaluated  on  the  aperture.  In  other  words,  we  stop  the  geo- 
metrical-optics rays  at  the  aperture  and  establish  an  aperture-distri- 
bution which  Is  nothing  more  than  the  equivalent  sources  In  the  Huygens' 
principle.  When  this  aperture-distribution  reradlates,  the  far-fleld 
beyond  the  aperture  can  be  recovered  as  If  the  rays  propagated  directly 
to  the  far-zone  undisturbed.  Then,  the  original  geometrical-optics 
far-fleld  pattern  should  be  recovered  Intact. 

However,  If  the  negative  of  the  Incident  field  is  assumed  In  the 
deep-shadow  region,  the  main  beam  of  the  scattered  far-fleld  pattern  is 
regained  as  shown  In  the  previous  two  subsections.  As  before,  for  the 
extraneous  oscillations  In  the  far-fleld  pattern  to  disappear,  the 
magnitude  of  the  aperture  distribution  must  not  contain  any  sharp  cor- 
ners or  discontinuities.  With  this  In  mind,  curve-fitting  techniques 
may  be  employed  In  the  transition-region  to  smoothly  Join  the  assumed 
approximate  aperture-distributions  In  the  well-lllumlnated  region  and 
the  deep-shadow  region. 

The  simplest  Interpolating  curve  that  fulfills  the  aforementioned 
requirement  for  suppressing  extraneous  oscillations  Is  the  cubic  spline 
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f(y)  “ . (6.29) 

The  derivative  of  f(y)  is 

f’(y)  “ 0 + C2  + Zc^y  + 3c^y^  . (6.30) 


The  coefficients  c^,C2»c^  and  c^  can  be  completely  determined  if  the 
values  and  the  slopes  of  the  interpolating  cubic  polynomial  are  knovm 
at  the  end  points  of  the  interpolation  Interval.  The  interpolation 
Interval  in  the  present  example  for  a circular  cylinder  with  radius  3X 
is  roughly  from  y * 2X  to  y * 6X,  which  is  obtained  by  observation  of 
the  width  of  the  transition  region  in  the  exact  aperture-field  distribution. 
The  values  and  slopes  of  the  magnitudes  and  phases  of  the  geometrical- 
optics  scattered  aperture-field  at  the  illuminated  end  and  the  negative 
of  the  incident  electric  field  at  the  dark  end  of  the  transition  region 
can  be  computed  numerically.  Then  the  coefficients  c^,C2»c^  and  c^ 
can  be  obtained  by  solving  the  following  matrix  equation: 


- 

— M 

* 

1 2X  (2X)^  (2X)^ 

*^1 

f(2X) 

0 1 2(2X)  3(2X)^ 

f’(2X) 

1 6X  (6X)^  (6X)^ 

*^3 

f(6X) 

0 1 2(6X)  3(6X)^ 

It 

f'(6X) 

(6.31) 


where 

f(2A)  is  the  value  of  the  magnitude  or  phase  of  the  negative  In- 
cident field  at  y ■ 2X, 

f'(2X)  is  the  slope  of  the  magnitude  or  phase  of  the  negative 
of  the  incident  field  at  y ® 2X, 
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f(6A)  is  the  value  of  the  magnitude  or  phase  of  the  geometrical- 

I 

optics  field  at  y > 6A, 

f'(6X)  Is  the  slope  of  the  magnitude  or  phase  of  the  geometrical- 
optics  field  at  y > 6X, 

and 

X Is  the  Incident  wavelength. 

t 

I Thus,  two  sets  of  coefficients,  one  for  the  magnitude,  and  the  other 

I for  the  phase,  are  obtained  from  (6.31)  by  first  substituting  the  mag- 

I nltude  entries  Into  the  right -hand- side  column-vector,  and  then  doing 

f 

i 

I the  same  thing  for  the  phase  entries. 

1 

I The  magnitude  and  phase  curves  are  chosen  to  be  Interpolated  be- 

cause  they  are  slowly  varying  smooth  curves.  The  curve  fitting  cannot 
j be  applied  to  the  real-  and  Imaginary-part  curves  of  the  same  complex 

quantity  because  they  are  highly  oscillatory  curves. 

To  summarize,  an  aperture-field  distribution  has  been  synthesized 
with  the  negative  of  the  Incident  electric  field  In  the  deep-shadow 

I 

j region,  the  geometrical-optics  scattered  field  In  the  well-lllumlnated 

> 

I region,  and  the  cublc-spllne  Interpolations  (applied  to  the  magnltude- 

I 

^ and  phase-curves)  In  the  transition-region.  The  synthesis  Is  sche- 

matically shown  In  Figure  6.12. 

The  scattered  far-fleld  pattern  Is  computed  In  two  steps:  (1)  For 

i 

I the  forward-scattering  hemisphere,  the  far-fleld  pattern  is  obtained  by 


Fourier  transforming  the  synthesized  aperture  distribution  via  the  well- 
known  relationship  between  the  far-zone  field  and  the  aperture  field. 

(2)  For  the  back-scattering  hemisphere,  the  far-fleld  pattern  Is  simply 
obtained  from  the  geometrical-optics  solution.  The  Fourier  transform 
can  be  performed  efficiently  by  a one-dimensional  fast  Fourier  trans- 
form (FFT). 
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However,  attention  must  be  called  to  the  fact  that  the  FFT  result 
Is  forced  to  zero  at  90°  and  270°  because  the  aperture  must  be  trun- 
cated to  a finite  width  (even  as  large  as  90X,  thirty  times  the  radius 
of  the  cylinder)  in  order  to  apply  the  FFT.  Nevertheless,  the  FFT  re- 
sult Is  extremely  good  In  the  forward-scattering  direction.  As  mentioned 
In  the  previous  section,  the  geometrical-optics  far-fleld  Is  well-de- 
fined and  Is  accurate  for  a wide  range  of  observation  angles  extending 
from  the  lit  side  Into  the  dark  side.  It  Is  found  that  the  FFT  result 
and  the  geometrical-optics  far-fleld  overlap  for  a relatively  wide  rar.je 
of  angles  (about  20  degrees)  In  the  neighborhood  of  70°  and  290°,  hence, 
the  two  far-fleld  patterns  can  be  Joined  smoothly  at  70°  and  290°  with- 
out any  discernible  discontinuities. 

The  far-fleld  pattern  obtained  In  the  above  manner  Is  shown  In 
Figure  6.12.  Except  for  the  slight  differences  In  the  directions  around 
50°  and  310°,  the  pattern  Is  In  excellent  agreement  with  the  exact  solu- 
tion, which  demonstrates  that  the  Inclusion  of  the  phase  Information  In 
the  transition  region  Is  crucial  In  obtaining  an  accurate  approximation 
to  the  exact  solution.  Due  to  the  slow  variation  of  the  magnitude  and 
phase  of  the  aperture-field  distribution,  the  phase  Information  Is  ade- 
quately obtained  by  a relatively  straightforward  cublc-spllne  Inter- 
polation. In  the  next  subsection,  we  discussed  another  method  of  syn- 
thesizing the  aperture-distribution  so  that  the  computation  of  the 
aperture  field  can  be  reduced  substantially  while  the  acciiracy  In  the 
far-fleld  pattern  is  still  preserved. 


6.4.4  Method  IV 


In  the  last  subsection,  an  aperture-field  distribution  has  been 
synthesized  with  the  negative  of  the  incident  electric  field  in  the 
deep-shadow  region,  the  geometrical-optics  scattered  field  in  the  well- 
illuminated  region,  and  the  cubic-spline  interpolations  (applied  to  the 
magnitude-  and  phase-curves)  in  the  transition  region.  The  scattered 
far-fleld  pattern  obtained  by  that  method  was  in  good  agreement  with 
the  exact  solution. 

However,  a reevaluation  of  the  synthesis  in  Method  III  reveals 
that  the  computation  of  the  geometrical-optics  scattered  field  on  the 
aperture  is  redundant.  Recall  that  in  Method  I or  Method  II,  the 
aperture-field  distribution  in  the  well-illuminated  region  was  assumed 
to  be  zero,  and  the  far-field  pattern  in  the  forward-scattering  hemi- 
sphere was  obtained  by  the  superposition  of  the  geometrical-optics 
far-field  and  the  far-field  radiated  by  the  assumed  aperture-field  dis- 
tribution in  the  geometrical  shadow  region.  The  superposition  of  the 
two  far-flelds  in  such  a manner  is  equivalent  to  an  alternative  approach 
that  was  mentioned  in  Subsection  6.4.3,  l.e.,  the  Huygens'  equivalent 
sources  for  the  geometrical-optics  far-fleld  are  created  on  the  well- 
illuminated  portion  of  the  aperture,  and  combined  with  the  assumed 
aperture-field  distribution  on  the  complementing  nonlllumlnated  portion 
of  the  aperture  to  obtain  a complete  scattered  electric  aperture-field 
distribution  which,  when  reradlated  into  the  forward-scattering  hemi- 
sphere, produces  the  scattered  far-fleld  pattern  in  that  hemisphere. 

Theoretically,  these  two  different  approaches  for  obtaining  the 
scattered  far-field  in  the  forward-scattering  hemisphere  give  the  same 
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result.  From  a computational  standpoint,  however,  the  approach  using 
the  geometrical-optics  scattered  field  directly  In  the  far-zone  for 
superposition  Is  more  efficient  than  the  other  approach.  The  reason  for 
Its  superiority  Is  two-fold:  First,  computational  time  can  be  reduced 

substantially  because  the  aperture-field  distribution  In  the  Illuminated 
portion  of  the  aperture  need  not  be  computed  but  Is  simply  set  equal  to 
zero;  and  second,  the  fast  Fourier  transform  can  be  applied  to  Its  full 
capacity,  meaning  no  Inaccuracies  near  90^  and  270*^,  because  the  aperture- 
field  distribution  Is  completely  confined  to  that  portion  of  the  aperture 
which  Includes  the  deep-shadow  plus  the  transition-region  only.  To 
further  emphasize  the  advantage  of  eliminating  the  calculations  of  the 
aperture-field  In  the  Illuminated  region,  we  discuss  the  labor  Involved 
In  setting  up  an  aperture-field  distribution  In  the  next  paragraph. 

In  order  to  obtain  a facsimile  of  the  Fourier  transform  of  the  In- 
finite aperture-distribution  using  the  fast  Fourier  transform  (FFT), 
the  aperture  must  be  truncated  at  points  beyond  which  the  aperture-field 
has  decreased  to  a negligible  quantity;  which  means  that  the  Illuminated 
portion  of  the  aperture  Is  typically  as  large  as  fifteen  times  the  geo- 
metrical-optics shadow  on  the  aperture.  As  mentioned  earlier  In  Section 
6.3  where  the  geometrical-optics  solution  was  briefly  derived,  there 
Is  no  simple  relation  that  can  be  found  to  explicitly  express  In  terms 
of  the  polar  coordinates  (p,<^)  the  parameters  6^,  1,  and  In  the  geo- 
metrical-optics scattered-field  expression.  The  far-fleld  approximation 
(p  ->  <»)  cannot  be  directly  employed  to  accurately  calculate  the  geo- 
metrical-optics aperture-field  distribution  because  the  aperture-plane 
Is  In  juxtaposition  with  the  cylinder;  hence,  for  a considerably  large 
portion  of  the  aperture,  the  far-fleld  approximation,  e.g.,  (6.27),  Is 


still  too  crude.  Therefore,  In  practice,  one  has  to  resort  to  numerical 
methods.  There  are  basically  two  routes  that  one  can  follow,  namely, 
either  searching  for  the  point  of  reflection  when  a point  of  observation 
on  the  aperture  Is  given,  or  Interpolating  for  the  field  when  a point 
of  reflection  Is  chosen.  While  the  accuracies  of  the  numerical  results 
obtained  by  these  two  different  approaches  are  comparable,  the  latter  Is 
more  straightforward  In  computer  programming  and  more  efficient  compu- 
tation-wise. But  In  either  case,  the  geometrical-optics  scattered 
aperture-field  has  to  be  computed. 

In  this  subsection,  we  discuss  a synthesizing  method  which  Is,  In 
part,  similar  to  Methods  I and  II,  l.e.,  the  calculations  of  the  aperture 
field  In  the  Illuminated  region  are  eliminated.  However,  the  major 
difference  Is  the  Inclusion  of  the  phase  Information  In  the  transition- 
region  In  the  present  synthesis.  As  In  all  of  the  previously  discussed 
methods,  the  scattered  electric  aperture-field  distribution  In  the  deep- 
shadow  region  Is  assumed  to  be  the  negative  of  the  Incident  field.  The 
field  distribution  In  the  well-lllumlnated  portion  of  the  aperture  Is  set 
equal  to  zero.  In  the  transition-region,  curve-fitting  techniques  are 
employed  to  smoothly  Join  the  assumed  approximate  aperture-distributions 
In  the  well-lllumlnated  and  the  deep-shadow  regions.  A cubic  spline  Is 
obtained  (for  the  magnitude  of  the  aperture-field)  by  simply  substituting 
f(6X)  « 0.0  and  f*(6A)  « 0.0  Into  the  rlght-hand-slde  of  the  column- 
vector  In  (6.31)  and  solving  the  resulting  matrix-equation  for  the  un- 
known coefficients.  The  Interpolated  curve  Is  the  gradually  tapered 
magnitude  curve  of  the  aperture-field  distribution.  It  seems  the  phase  . 
curve,  at  the  first  look,  cannot  be  obtained  as  straightforwardly  as  the 
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magnitude  curve,  because  the  assumed  aperture-field  at  the  Illuminated 
end  of  the  transition-region  is  zero  and  the  value  and  slope  of  the 
phase  curve  are  unknown.  However,  the  value  and  slope  of  the  phase  of 
the  geometrical-optics  scattered  field  at  the  illuminated  end  of  the 
transition  region  can  be  used  in  (6.31)  for  obtaining  the  phase-informa- 
tion in  the  transition  region.  This  point  is  discussed  further  in  the 
next  paragraph. 

In  particular,  the  aperture-field  being  synthesized  is  the  difference 
of  two  aperture-field  distributions,  namely,  the  aperture-field  synthe- 
sized by  Method  III  and  the  geometrical-optics  aperture-field,  i.e.. 


AP 


1^ 


IV 

AP 


- E 


III 

'AP 


1<P 


III 

AP 


G.O.  i(t> 


G.O. 


- E 


AP 


AP 


(6.32) 


where  E is  the  magnitude  and  Is  the  phase  of  the  complex  quantity. 
Superscripts  IV,  III,  and  G.O.  stand  for  Method  IV,  Method  III,  and 
geometrical  optics,  correspondingly.  Subscript  AP  stands  for  aperture. 
Note  that 


(6.33a) 


and 


III 
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(6.33b) 


which  are  direct  consequences  of  the  synthesis  in  Method  III.  Therefore, 
at  y ■ 6A,  we  have 
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(6.34b) 


This  confirms  the  remarks  made  In  the  last  paragraph  that  the  magnitude 
Is  zero  and  the  phase  Is  the  geometrical-optics  phase  at  the  Illuminated 
end  of  the  transition  region. 

Thus,  an  Interpolation  by  a cubic  spline  for  the  phase  In  the 
transition  region  gives  an  Identical  phase  variation  as  In  Method  III. 

Note,  however,  that  the  phase  curve  need  be  computed  only  In  the  deep- 
shadow  and  the  transition  region  where  the  magnitude  of  the  aperture 
distribution  Is  nonzero.  In  the  well-lllumlnated  region  of  the  aperture, 
the  phase  variation  becomes  Irrelevant  because  the  magnitude  of  the 
aperture  distribution  Is  zero. 

As  mentioned  previously,  the  Illuminated  portion  of  the  aperture 
Is  typically  as  large  as  fifteen  tiroes  the  geometrical-optics  shadow 
region  In  order  to  obtain  a facsimile  of  the  Fourier  transform  of  the 
Infinite  aperture  distribution  using  the  fast  Fourier  transform.  There- 
fore, In  setting  up  the  aperture  distribution  for  the  Fourler-transform 
operation,  the  present  method  Is  computationally  more  efficient  than 
Method  III.  In  addition,  the  fast  Fourier  transform  can  be  applied  to 
Its  full  capacity,  meaning  no  Inaccuracies  near  90°  and  270°,  because 
the  aperture-field  distribution  Is  completely  confined  to  that  portion 
of  the  aperture  which  Includes  the  deep-shadow  plus  the  transition  re- 
gion only.  This  renders  It  unnecessary  to  Join  the  two  far-fleld  patterns 

at  70°  and  290°  as  In  Method  III.  ? 

i 

i 

I 
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To  summarize,  an  aperture-field  distribution  has  been  synthesized 
with  the  negative  of  the  incident  electric  field  in  the  deep-shadow 
region,  the  zero  field  in  the  well-illuminated  region,  and  the  cubic- 
spline  interpolations  (one  of  which  is  between  the  magnitude  of  the 
incident  field  and  zero-field  for  the  magnitude-variation,  and  the  other 
is  between  the  phase  of  the  incident  field  and  the  phase  of  the  geometrical- 
optics  field  for  the  phase-variation)  in  the  transition  region.  The 
synthesis  is  schematically  shown  in  Figure  6.12. 

The  scattered  far-fleld  pattern  is  computed  in  two  steps: 

(1)  For  the  forward-scattering  hemisphere,  the  far-fleld  pattern  is  ob- 
tained by  Fourier  transforming  the  synthesized  aperture  distribution 
via  the  well-known  relationship  between  the  far-zone  field  and  the  aper- 
ture field.  The  Fourier  transformed  result  is  then  superposed  on  the 
geometrical-optics  far  field  to  obtain  the  scattered  far-fleld  pattern 
in  that  hemisphere.  (2)  For  the  back-scattering  hemisphere,  the  far- 
fleld  pattern  is  simply  obtained  from  the  geometrical-optics  solution. 

Note  that  the  geometrical  optics  far  field  is  used  throughout  the  en- 
tire range  of  observation  angles  (0  to  360*^)  because  in  the  synthesis 
of  the  aperture-field  distribution,  the  geometrical-optics  aperture 
field  has  not  been  taken  into  account  (see  (6.32));  hence,  in  step  (1) 
above,  a direct  superposition  of  the  far  fields  is  necessary  to  obtain 
the  scattered  far-fleld  pattern  in  the  forward-scattering  hemisphere. 

The  far-fleld  pattern  obtained  in  the  above  manner  is  shown  in 
Figure  6.12.  Observe  that  the  positions  of  the  maxima  and  minima  are 
in  good  agreement  with  the  exact  solution.  However,  the  oscillations 
are  not  as  pronounced  and  there  are  slight  discrepancies  (about  ±1  dB) 
in  the  levels  of  the  maxima  and  minima  in  comparison  with  the  exact 
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solution.  Nevertheless,  the  far-fleld  pattern  of  this  method  Is  In 
better  agreement  with  the  exact  pattern  In  the  directions  around  50° 
and  310°  than  the  far-fleld  pattern  of  Method  III.  Generally  speaking, 
the  overall  far-fleld  pattern  Is  In  reasonably  good  agreement  with  the 
exact  solution. 

The  slight  discrepancies  In  the  far-fleld  pattern  are  due  to  the 
failure  to  take  Into  account  the  shadowed  portion  of  the  cylindrical 
surface  In  the  solution  process.  As  a matter  of  fact,  the  shadowed 
portion  of  the  obstacle  can  be  of  an  arbitrary  shape,  but  as  long  as 
the  Illuminated  portion  of  the  obstacle  Is  circularly  cylindrical,  the 
far-fleld  pattern  obtained  using  the  present  method  (or  any  one  of  the 
methods  discussed  so  far)  will  be  the  same  as  that  of  the  circular 
cylinder.  In  the  next  subsection,  we  discuss  the  use  of  Galerkln's 
method  applied  In  the  spectral  domain  to  further  Improve  the  scattered 
far-fleld  pattern  by  taking  Into  account  the  Information  characterizing 
the  shadowed  portion  of  the  obstacle,  e.g.,  the  radius  of  curvature 
of  the  shadowed  surface.  In  the  process  of  solution. 

6.4.5  Method  V 

In  the  last  subsection,  we  presented  a method  of  synthesizing  the 
aperture-field  distribution  so  that  the  redundant  computations  of  the 
aperture-field  In  the  well-lllumlnated  portion  of  the  apertu  were 
completely  eliminated  while  the  accuracy  In  the  far-fleld  pattern  was 
still  preserved  to  an  extent  that  the  overall  far-fleld  pattern  was 
In  reasonably  good  agreement  with  the  exact  solution.  The  slight  dis- 
crepancies In  the  far-fleld  pattern  were  due  to  the  failure  to  take  Into 
account  the  shadowed  portion  of  the  surface  of  the  obstacle  In  the 
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process  of  solution.  In  this  subsection,  we  discuss  a method  of  synthe- 
sizing the  aperture-field  distribution  taking  Into  account  the  Informa- 
tion characterizing  the  shadowed  portion  of  the  obstacle  In  the  process 
of  solution  so  that  further  lnq>rovement  of  the  scattered  far-fleld  pattern 
can  be  achieved. 

As  alluded  to  In  Section  6.1,  the  Introduction  of  this  chapter, 
a marked  difference  between  the  present  approach  and  all  of  the  con- 
ventional methods  available  In  the  literature  Is  the  Introduction  of  an 
Infinite  planar  aperture  In  juxtaposition  with  the  obstacle.  The  scat- 
tered far-fleld  pattern  In  the  forward-scattering  hemisphere  Is  obtained 
by  a Fourier  transformation  of  the  scattered  field  distributions  on  this 
aperture.  Therefore,  operations  for  Improving  the  scattered  far-fleld 
pattern  can  be  performed  on  the  scattered  aperture-field  distribution. 

In  contrast  to  the  conventional  approaches  In  which  they  are  done  either 
directly  on  the  Induced  surface  current  or  directly  on  the  scattered 
far  field.  The  advantage  of  concentrating  our  efforts  on  Improving 
the  aperture-field  distribution  Is  two-fold.  Firstly,  the  magnitude 
and  the  phase  of  the  aperture-field  distribution  are  smoothly,  non- 
osclllatorlly  varying  functions  of  positions  In  the  plane  of  the  aper- 
ture; hence,  Improvement  can  be  achieved  In  a relatively  straightforward 
manner.  Secondly,  the  operations  for  Improvement  are  performed  on  a 
planar  surface  rather  than  on  a curved  surface  as  In  the  case  of  at- 
tempting to  directly  Improve  the  surface  current  on  the  obstacle  In  a 
conventional  approach,  hence,  the  present  method  Is  conceptually  simpler. 
An  Improved  scattered  far-fleld  pattern  In  the  forward  scattering  hemi- 
sphere of  an  n-dlmenslonal  obstacle  can  be  obtained  by  Fourier  trans- 
forming an  Improved  (n  - 1) -dimensional  scattered  aperture-field 


distribution,  where  n « 2,3.  This  is  a salient  feature  of  the  method 
which  manifest  Itself  to  be  numerically  efficient  by  the  circumvention 
of  the  unwieldy  three-dimensional  fast  Fourier  transform  (FFT) . These 
Important  points  have  been  demonstrated  in  the  previously  discussed 
methods  of  synthesizing  the  aperture-field  distribution. 

Because  of  the  aonosclllatory  nature  of  the  magnitude  and  the  phase 
variations  of  the  scattered  aperture-field  distribution,  a very  close 
approximation  to  the  exact  aperture-field  distribution  was  obtained  by 
Method  111  in  which  a simple  cubic-spline  interpolation-scheme  was  used 
in  the  transition  region  to  join  the  negative  of  the  incident  field 
in  the  deep-shadow  region  smoothly  to  the  geometrical  optics  scattered 
field  in  the  well-illuminated  region.  Again,  a similar  interpolation- 
scheme  was  employed  in  the  synthesis  of  the  aperture-field  distribution 
in  Method  IV  in  which  the  differenae  of  the  aperture-field  distribution 
of  Method  111  and  the  geometrical-optics  scattered  aperture-field  was 
closely  approximated  by  joining  the  negative  of  the  incident  field  in 
the  deep-shadow  region  smoothly  to  the  zero  "difference-field"  in  the 
well-illuminated  region.  However,  the  magnitude  and  the  phase  varia- 
tions in  the  transition  region  are  slightly  different  from  the  inter- 
polating cubic  polynomials.  Although  the  approximations  made  by  these 
cubic  polynomials  are  very  close  to  the  true  variations  of  the  magni- 
tude and  the  phase  of  the  aperture-field  distribution,  these  slight  de- 
viations manifest  themselves  in  the  resulting  scattered  far-fleld  pat- 
terns shown  in  Figure  6.12. 

.he  idea  to  further  improve  the  scattered  far-fleld  pattern,  or 
equivalently,  the  synthesis  of  the  aperture-field  distribution  of 
Method  IV,  is  to  expand  the  aforementioned  slight  deviations  of  the 
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aperture-field  distributions  In  the  transition  region  Into  a series  of 
known  basis  functions  with  unknown  coefficients.  The  unknown  coef- 
ficients can  be  determined  by  enforcing  the  boundary  conditions  on 
the  shadowed  portion  of  the  cylindrical  surface  of  the  obstacle.  This 
particular  portion  of  the  surface  of  the  obstacle  Is  chosen  because  It 
has  the  strongest  effect  on  the  variations  of  the  aperture-field  dis- 
tribution In  the  transition  region.  In  a sense,  this  Is  an  application 
of  the  Galerkln's  method  In  the  spectral  domain  to  solve  for  the  un- 
known deviations  of  the  first-order  synthesis  of  the  scattered  aperture- 
field  distribution  from  the  true  variations  of  the  scattered  aperture- 
field  distribution  to  obtain  an  Improved  second-order  synthesis.  Note, 
however,  that  the  unknown  quantities  are  the  deviations  from,  but  not 
the  true  variations  themselves.  Since  the  first-order  synthesis  of 
Method  IV  Is  already  a reasonably  good  approximation  to  the  true  solution, 
only  a refinement  of  the  synthesis  is  necessary;  and  the  Galerkln's 
method  applied  In  the  spectral  domain  Is  well-suited  for  this  purpose. 

A word  of  caution  Is  In  order  In  regard  of  the  choice  of  basis 
functions  In  the  Galerkln's  method.  Rectangular  pulses  are  not  recom- 
mended for  the  basis  because  they  Introduce  abrupt  discontinuities  Into 
the  aperture-field  distribution  and  give  rise  to  the  extraneous  oscil- 
lations In  the  scattered  far-fleld  pattern.  This  point  has  been  demon- 
strated In  Method  II  of  Subsection  6.4.2.  Therefore,  only  smoothly 
varying  and  gradually  tapering  off  functions  are  suited  for  the  use 
as  basis.  With  this  In  mind,  a Gaussian  pulse, as  shown  In  Figure  6.13, 

Is  used  as  the  basis  function,  l.e., 

- exp(-a\2) 


(6.35) 
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where  a - 2.75  and  Is  a local  coordinate.  This  particular  choice  of 
the  spread  of  the  Gaussian  pulse,  1/a,  Is  to  assure  that  the  Gaussian 
pulse  decreases  to  1 percent  of  Its  peak  value  at  x^  ■ 0.78,  and  It 
decreases  to  approximately  1/e  of  Its  peak  value  at  the  crossover  point 
of  two  adjacent  Gaussian  pulses,  which  means  that  the  gaps  between  ad- 
jacent Gaussian  pulses  are  filled  up  properly  and  the  Interferences 
with  Gaussian  pulses  beyond  the  adjacent  ones  are  practically  eliminated. 
Thus,  a superposition  of  these  basis  functions  results  In  a smoothly 
varying  curve,  which  Is  precisely  what  Is  wanted  In  the  refinement  of 
the  aperture-field  distribution. 

In  applying  the  Galerkln's  method  In  the  spectral  domain,  the 
Fourier  transform  of  the  basis  function  Is  needed.  For  the  particular 
choice  of  a Gaussian  pulse  In  (6.35)  as  a basis  function.  Its  Fourier 
transform  Is  simply  another  Gaussian  function,  namely. 


(6.36) 


which  Is  shown  In  Figure  6.14.  In  the  above  equation,  Is  the  Fourier 
transform  variable  corresponding  to  the  spatial  variable  x^. 

The  testing  function  In  the  Galerkln's  method  Is  the  same  as  the 
basis  function.  The  locations  of  the  basis  functions  on  the  aperture 
and  the  testing  functions  on  the  surface  of  the  obstacle  are  schematically 
shown  In  Figure  6.15.  Observe  that  the  basis  functions  are  situated 
within  the  transition  region  In  addition  to  the  cublc-spllne  interpolation 
of  Method  IV.  The  testing  functions  are  arranged  In  the  manner  shown 
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Figure  6.14.  Fourier  transform  of  the  Gaussian  pulse  in  Figure  6.13, 


I-'-. 


in  Figure  6.15  because  that  particular  portion  of  the  surface  of  the 
obstacle  has  the  . urongest  effect  on  the  variations  of  the  aperture- 
field  distribution  In  the  transition  region.  Due  to  the  symmetry  of  the 
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problem,  N basis  functions,  ****  used  but  only  N/2  un- 

known coefficients  need  be  solved,  hence  N/2  testing  functions, 

• • • » ^n/2‘ 

The  aperture- field  being  synthesized  Is 


n*l 


®n^^n  ^(N/2+n)^ 


(6.37) 


“iV 

where  Is  the  difference  aperture-field  distribution  of  Method  IV 

-IV 

(see  (6.32)),  and  the  deviations  of  E._  from  the  true  variations  are 

AP 

approximated  by  the  series  In  basis  function  of  the  form  given  In  (6.35), 
or  specifically, 

fn^y)  “ exp[-a^(y  - ygn)^]  » " “ ...,  N (6.38) 


where  yQ^  Is  the  position  at  which  the  peak  of  the  nth  basis  function 

Is  located  (see  Figure  6.15).  The  unknown  coefficients  c 's  are  to  be 

n 

determined  by  the  Galerkln's  method  In  the  spectral  domain. 

Consider  the  aperture  for  the  hemisphere  of  a testing  function 
as  shown  in  Figure  6.16.  The  portion  of  the  aperture  underneath  the 
Gaussian  pulse  Is  embedded  In  the  obstacle.  The  width  of  this  embedded 
portion  of  the  aperture  Is  determined  by  the  particular  choice  of  the 
spread  of  the  Gaussian  In  (6.35)  so  that  It  Is  subtended  by  an  angle 
of  30°  as  shown  In  Figure  6.16.  The  scattered  E-fleld  on  this  embedded 
portion  of  the  aperture  must  be  the  negative  of  the  incident  E-fleld, 


•-t»i5n«5vn»  ' 
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-E^,  to  satisfied  the  extended  boundary  condition  [31].  The  scattered 
E-fleld  outside  of  this  portion  of  the  aperture  Is  unknown  and  Is  desig- 
nated by  F as  shown  In  Figure  6.16.  Therefore,  the  scattered  E-fleld 
in 

on  this  aperture  can  be  written  as 


E®._  - -E^  + F 
mAP  m m 


(6.39) 


where  the  subscript  m signifies  the  aperture  corresponding  to  the  mth 

testing  function,  P . The  Fourier  transform  of  (6.39)  reads 

m 


-E  + F 
m m 


(6.40) 


where  " on  top  Indicates  the  Fourier  transform  of  the  corresponding  term 
in  (6.39). 

Now  we  take  the  Inner  product  of  the  Fourier  transform  of  P ,P  , 
with  (6.40),  which  results  In  the  following  equation: 


^^m’^mAP^ 


-<P  ,E^>  + <P  ,F  > 
m m m m 


(6.41) 


where  the  angle  brackets  Indicate  an  Inner  product  of  two  complex  func- 
tions of  e.g.,  A(|3)  and  B(3)*  defined  as 

“ * 

<A,B>  - / A (e)B(3)  d6  (6.42) 


where  the  asterisk  denotes  complex  conjugate.  By  Parseval's  theorem 
the  second  term  In  the  right-hand- side  member  of  (6.41)  can  be  written 
as 


<P 

m 


> 


F > 
m 


2tt<P  , 
m 


F > 
m 


0 


(6.43) 
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because  P and  F are  defined  In  complementary  domains  on  the  aperture  for 
m m 

the  hemisphere  of  (see  Figure  6.16).  Therefore,  (6.41)  becomes 


> 


(6.44) 


The  Fourier  transform  of  the  aperture  E-fleld,  Is  related  to  the 
fop- field  pattern  in  the  hemiaphere  of  following  re- 
lation: 


2 

k sin 


W 


(6.45) 


where  ~ being  the  local  polar  coordinate,  ((>  the  central 

polar  coordinate,  and  the  polar  coordinate  of  the  location  of  P^, 
all  of  which  are  shown  In  Figure  6.16.  T ((j».)  Is  obtained  by  truncating 

d Jw 

V 

the  complete  scattered  far-fleld  pattern,  T (4>j^)>  bo  the  hemisphere  of 
Pj^.  t'^((|))  Is  given  by 

V IV  f f 

r (<fr)  - T^\^)  + I c^(T^(^)  + 

n-1  ' 

where  T^^((|))  Is  the  scattered  far-fleld  pattern  of  Method  IV;  and  T^(<l>)  Is 
the  scattered  far-fleld  pattern  of  the  nth  basis  function,  f^.  In  Figure 
6.15.  More  specifically,  T^(<l>)  Is  expressed  In  terms  of  the  Fourier 
transform  of  the  nth  basis  function,  f jj»  in  the  following  manner: 


(♦)  . ^ .-lkAcos4  ^ 


n 


0 <.  <)  < 90°  and  270°  < ^ < 360° 


(6.47) 
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in  which  A is  the  radius  of  the  cylinder,  and  f^(4>)  1^  obtained  from 
(6.36)  using  the  shifting  relation  between  Fourier  pairs;  hence, 


, n -k^  sln^ 

*„w)  ■ - — — 


n 


1,2, 


N 


(6.48) 


4a 


where  Is  the  position  at  which  the  peak  of  the  nth  basis  function, 
f^.  Is  located  (see  Figure  6.15). 

Using  (6.48),  (6.47),  (6.46),  and  (6.45)  In  (6.44),  an  equation  of 
N/2  unknowns,  C2.*^2 %/2*  obtained.  Repeating  the  above  pro- 

cedure with  different  locations  of  the  testing  function  on  the  surface 
of  the  obstacle  results  In  a system  of  N/2  linear  equations  of  N/2  un- 
knowns, which  Is  then  solved  by  standard  matrix-inversion  techniques. 

The  matrix  elements  are  generated  by  the  Inner  products  of  (6.44) 

In  a straightforward  manner.  The  Inner  product  on  the  left-hand-side 
of  (6.44)  simply  becomes  a linear  superposition  of  the  Inner  products  of 

P with  the  Individual  far-fleld  patterns.  This  can  be  seen  from  (6.45) 
m 

and  (6.46).  Since  the  Invisible  spectra  of  the  far-fleld  patterns  are 
highly  attenuated  and  they  hardly  contribute  to  the  Inner  products, 
therefore,  they  can  be  Ignored  In  the  computation  of  the  Inner  products. 
Hence,  the  formula  for  computing  the  Inner  product  of  and  T Is 


*11"’^  lksln(()  Acosl5°  /-  j-k^  cos^  (f.  \ 

- I e ^ exp  5 ^ 2T(<^  ) A*  (6.49) 

♦,l-0  = \ 4a^  I “ 


IV  f 

where  T stands  for  T or  T . 

n 


Using  Parseval's  theorem,  the  left-hand-side 
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of  (6.44)  can  be  written  as 


-<P  ,E^>  - -2tt<P  ,E^> 
m m mm 


Since  P^  Is  assumed  nonzero  only  on  the  embedded  portion  of  the  aperture 
for  the  hemisphere  of  P^,  (Figure  6.16),  the  formula  for  computing  the 


Inner  product  In  (6.50)  Is 


<P  ,E^> 
m m 


Xjj^^«Aslnl5 


x^^«-Aslnl5 


exp(-a  x^^  + Ikx^^cos  (J)^  - IkA  cos  15°  sin  (J)^)4x^^.  (6.51) 


After  the  unknown  coefficients,  Cj^»C2 ^/2*  been  deter- 

mined by  the  Galerkln's  method,  we  can  substitute  them  Into  (6.37)  to 

obtain  the  Improved  synthesis  of  the  "difference"  scattered  aperture- 

-V 

field  distribution,  E^p,  and  then  proceed  as  In  Method  IV  to  compute 
the  far-fleld.  Alternatively,  the  coefficients,  °N/2’ 

be  substituted  Into  (6.46)  to  directly  obtain  the  Improved  far-fleld 
pattern,  T^(it)).  The  far-fleld  pattern  obtained  by  this  method  with 
N/2  ■ 7 and  ■ -15°  to  -45°  at  5°  Increments  Is  shown  In  Figures  6.17a 
and  6.17b  and  Is  In  excellent  agreement  with  the  exact  far-fleld  pattern. 


L6 


Figure  6.17a. 


Scattered  far-field  pattern  in  dB  of  a circular  cylinder 
with  radius  a - 3X  obtained  by  Method  V. 


Figure  6.17b.  Scattered  far-field  pattern  in  dB  of  a circular  cylinder 
with  radius  a “ 6A  obtained  by  Method  V. 


6.5  Computation  of  Induced  Surface  Current 


In  electromagnetic  wave  scattering,  two  quantities  are  of  vital 
Importance,  one  Is  the  scattered  far-fleld  pattern,  from  which  the  bl- 
statlc  radar-cross-section  (RCS)  can  be  determined,  and  the  other  Is  the 
Induced  surface  current  on  the  scatterer,  from  which  Information  can 
be  extracted  to  aid  In  the  design  of  the  mounting  of  antennas  and  sensor 
booms  on  a practical  scattering  obstacle,  e.g.,  a satellite.  A detailed 
discussion  of  the  computation  of  the  scattered  far-fleld  pattern  by 
systematically  synthesizing  the  aperture-field  distribution  has  been 
presented  In  the  last  section.  An  excellent  approximation  to  the  exact 
scattered  far-fleld  pattern  has  been  obtained  In  a relatively  straight- 
forward manner  by  the  introduction  of  an  aperture  on  which  the  Improve- 
ment operations  were  performed.  Such  an  accurate  scattered  far-fleld 
pattern  motivates  us  to  develop  a feasible  technique  for  the  computation 
of  the  Induced  surface  current  on  the  obstacle  to  a reasonable  degree 
of  accuracy. 

In  the  conventional  approaches  In  the  literature,  the  calculation 
of  the  scattered  far-fleld  pattern  from  an  Induced  surface  current  may 
be  formidable  because  of  the  numerical  multiple-integration  Involved; 
on  the  other  hand,  the  evaluation  of  the  Induced  surface  current  di- 
rectly from  an  approximate  scattered  far-fleld  pattern  may  not  be  fea- 
sible because  such  a far-fleld  pattern  Is  not  accurate  enough.  The 
development  of  the  present  approach  provides  a means  to  evade  such  a 
dilemma  so  that  the  Induced  surface  current  can  be  evaluated  with  ease 
when  the  approximate  scattered  far-fleld  pattern  has  been  Improved. 

This  method  of  computation  of  the  Induced  surface  current  Is  made 
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possible  by  the  extremely  close  approximation  of  the  true  scattered 
£ar-fleld  pattern  obtainable  by  the  method  of  synthetic  aperture  dis- 
tribution discussed  in  the  last  section.  The  detailed  development  of 
the  method  for  computing  the  Induced  surface  current  is  presented  in 
the  next  subsection. 


6.5.1  Method  of  computation 

Suppose  the  Induced  surface  current  at  a point  Q on  the  surface  of 
the  obstacle  as  shown  In  Figure  6.18  is  to  be  evaluated.  The  polar  co- 
ordinates of  the  point  Q are  simply  A being  the  radius  of  the 

circular  cylinder  and  ij)^  the  angular  displacement  from  the  y-axls  with 
positive  reference  direction  counterclockwise.  An  infinite  aperture 
with  Q as  the  point  of  tangency  to  the  cylinder  Is  established  for  the 
hemisphere  of  Q as  shown  schematically  In  Figure  6.18.  A local  rec- 
tangular coordinate  system  the  accompanying  local  polar 

coordinate  system  are  set  up  for  the  aperture.  The  Idea  Is  to 

use  that  portion  of  the  scattered  far-field  pattern  contained  In  the 
hemisphere  of  point  Q to  obtain  the  scattered  H-fleld  on  the  aperture. 

The  scattered  H-fleld  at  point  Q Is  subsequently  combined  with  the  In- 
cident H-fleld  at  point  Q to  obtain  the  total  surface  current  at  point 
Q.  For  other  points  on  the  surface  of  the  obstacle,  the  induced  surface 
currents  can  be  obtained  In  the  same  manner  by  using  their  corresponding 
hemispheres  and  apertures. 

The  relation  between  the  scattered  far-fleld  pattern  and  the  Fourier 
transform  of  the  scattered  magnetic  aperture  field,  H , Is  readily  ob- 
tained as 


T((j)j^)  exp(-lkA  sin  (p^) 


(6.52) 


.>  - w> 


HEMISPHERE 
OF  POINT  Q 


of  point  Q at  which  the  induced 
iated. 
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where  T(4>.)  is  the  scattered  £ar-£leld  pattern  In  the  proper  hemisphere, 

Zq  Is  the  free  space  characteristic  impedance,  and  is  the 

local  polar  coordinate.  The  phase  £actor  In  (6.52)  Is  to  assure  that 
the  aperture  at  ■ A but  not  at  ■■  0 Is  considered.  The  scattered 
magnetic  aperture  £ield,  H^,  Is  obtained  by  an  Inverse  Fourier  trans- 
£orm  o£  In  (6.52),  l.e., 

H®(Xj^)  - F'^{ii®(<^^(a))}  . (6.53) 


However,  only  the  value  o£  the  magnetic  scattered  £ield  at  ■ 0 Is 
o£  Interest.  There£ore,  Ignoring  the  negligible  contributions  £rom  the 
Invisible  range  o£  the  £ar-£leld  pattern,  H (x^  ■ 0)  can  be  expressed 
as  an  Integral  o£  the  £ollowlng  £orm: 


H®(x.  - 0)  - / T((Ji^)  sin  (j)^  d(t)^  . (6.54) 

0 

g 

The  numerical  evaluation  o£  H (x^  « 0)  Is  expedited  by  the  £ollowlng 
procedure:  the  Interval  o£  Integration  £rom  0 to  7r,[0,7r],  Is  divided 

Into  36  sub  Intervals,  each  o£  which  Is  of  length  tt/36.  Throughout  the 
entire  subinterval,  the  sinusoid,  sin  can  be  assumed  to  have  the 
constant  value  at  the  midpoint  of  the  subinterval.  Thus,  (6.54)  can 
be  written  as 


H®(x„ 


0) 


36 

« J sin(<f_)  ^ 
1-1 


A H- 
^Zl~72 


180 


(6.55) 
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where  Is  the  midpoint  of  the  1th  subinterval.  If  the  scattered 
far-fleld  pattern  Is  calculated  at  an  Increment  of  1°,  the  Index 

J "1»2»  •••»  5. 

Note  that  the  above  procedure  expedites  the  numerical  evaluation 
of  " 0)  without  sacrificing  the  fine  detail  of  the  scattered  far- 

fleld  pattern.  This  point  Is  worth-emphasizing  because  the  accuracy  of 
the  scattered  far-fleld  pattern  obtained  In  the  last  section  should  be 
preserved  as  much  as  possible.  Further  degradation  of  the  scattered 
far-fleld  pattern  may  give  Intolerable  numerical  results. 

The  Incident  H-fleld  at  the  point  Q Is  given  by 

H - ■=-  e ^ . (( 


Combining  (6.55)  with  (6.56),  and  using  the  relation  between  the  total 
surface  current,  J,  and  the  total  H-fleld  at  the  surface,  l.e., 

J ■ n X H 

- n X (H^  + H®)  (6.57) 

the  total  surface  current  at  the  point  Q(x^^  ” 0,y^  > A)  Is  obtained. 

This  completes  the  development  of  the  method  of  computation  of  the  induced 
surface  current. 


6.5.2  Results  and  comments 

The  magnitude  and  the  phase  of  the  total  surface  current  on  a per- 
fectly conducting  circular  cylinder  with  radius  ■ 3X  are  displayed  In 
Figure  6.19.  The  exact  eigenfunction  solution  of  the  total  surface  cur- 
rent Is  also  shown  In  the  figure  for  comparison.  It  Is  seen  that  the 


Figure  6.19.  Total  surface  current  on  a perfectly  conducting  circular  cylinder 
with  radius  a = 3X. 
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solution  of  this  study  deviates  slightly  from  the  exact  solution  In 
the  neighborhood  of  80°.  However,  the  small  deviation  occurs  at  the 
low  end  of  the  current  curve,  therefore.  It  is  tolerable  for  practical 
purposes.  The  rest  of  the  curve  Is  In  remarkable  agreement  with  the 
exact  solution. 

Before  closing  this  section.  It  Is  worthwhile  to  contrast  the 
present  approach  with  the  one  developed  by  Burnside  al.  [4  ].  Their 
method  uses  a combination  of  the  geometrical  theory  of  diffraction  and 
the  moment  method  to  solve  a circular  cylinder  scattering  problem, 
wherein  the  cylindrical  surface  Is  divided  Into  three  regions,  namely, 
the  physical  optics  region  In  the  well-lllumlnated  side,  the  6TD  region 
In  the  deep  shadow,  and  the  pulse  region  In  the  transition  region  be- 
tween the  lit  and  the  dark  sides.  It  Is  a typical  example  of  the  con- 
ventional approaches  of  attempting  to  directly  improve  the  surface  cur- 
rent. However,  one  finds  that  the  composite  surface  current  computed 
by  using  different  techniques  In  different  regions  on  the  cylindrical 
surface  as  Indicated  above  Is  not  smoothly  connected  at  the  bound- 
aries of  each  region.  These  discontinuities  are  clearly  shown  In 
Figure  6.20,  which  Is  Figure  4 In  [4  ] reproduced  here  for  comparison. 
Furthermore,  the  discontinuities  at  the  high  end  In  the  surface  cur- 
rent may  Introduce  extraneous  high-level  sldelobes  In  the  far-fleld 
pattern.  In  addition,  the  surface  current  of  Burnside's  method  de- 
viates from  the  exact  solution  in  both  the  high  and  the  low  ends  of 
the  current  curve.  Moreover,  the  pulse  region  In  Burnside's  method 
Is  typically  one  half  of  the  circumference  of  the  circular  cylinder; 
therefore,  many  current  samples  are  needed  for  an  electrically  large 
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cylinder,  especially  In  the  case  of  an  E-wave  Incidence  where  the 
surface  current  starts  to  deviate  from  the  physical  optics  current 
much  faster  than  It  does  In  the  H-wave  case. 

In  Burnside's  approach,  the  far-fleld  pattern  must  be  computed 
by  a numerical  double  Integration  of  the  surface  current.  Such  a com- 
putation may  be  very  time  consuming.  In  contrast,  we  have  presented  a 
conceptually  simple,  computationally  efficient  method  to  obtain  both 
the  scattered  far-fleld  pattern  and  the  Induced  surface  current, which 
are  In  remarkably  good  agreement  with  the  true  solution. 

6.6  Accuracy  Check 

In  the  previous  sections,  methods  have  been  discussed  for  obtaining 
scattered  far-fleld  patterns  and  Induced  surface  currents.  The  results 
have  been  compared  with  the  available  exact  solution  so  that  the  ac- 
curacies of  the  approximate  solutions  were  verified.  However,  In  many 
practical  situations,  the  exact  solution  Is  Impossible  to  obtain,  and 
quite  often,  the  solutions  arrived  at  by  other  approaches  may  not  be 
accurate  enough  to  compare  with.  Of  course,  experimental  results  pro- 
vide a good  evidence  to  validate  the  computed  results.  Nonetheless, 
experiments  are  expensive  to  set  up  and  also  can  be  time-consuming. 

Moreover,  special  techniques  and  Ingenuities  may  be  required  to  extract 
the  needed  Information  from  the  experimental  data,  viz.,  the  measurable 
quantities.  In  addition,  many  environmental  factors  may  Inflict  adverse 
effects  on  the  measurable  quantities  and  affect  the  results  of  the  ex- 
periment. Therefore,  an  economical  means  of  validating  a computed 
solution  Is  very  much  desired. 
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In  the  previous  chapters,  the  approximate  solutions  of  the  scattered 
far-fleld  patterns  or  Induced  surface  currents  were  validated  by  a self- 
checking procedure,  l.e.,  the  tangential  components  of  the  scattered 
E-flelds  on  the  surface  of  the  scatterers  were  computed  and  compared 
with  the  tangential  components  of  the  Incident  E-flelds  to  check  If 
the  boundary  conditions  were  satisfied.  For  scatterers  with  planar 
facets,  the  accuracy  checking  method  developed  In  the  previous  sections 
Is  feasible.  On  the  other  hand,  for  scatterers  with  smoothly  varying 
surface,  or  more  precisely,  surface  with  continuously  varying  curvatures, 
the  accuracy  checking  method  developed  In  the  previous  chapter  becomes 
Inefficient  because  a large  number  of  points  on  the  surface  of  the  scat- 
terer  need  be  checked  and  each  application  of  the  method  yields  rele- 
vant Information  at  a single  point,  viz.,  the  point  of  tangency  of  an 
aperture  and  the  obstacle.  In  order  to  provide  a more  efficient  means 
of  checking  the  accuracy  of  the  approximate  solutions,  the  self-checking 
procedure  of  the  previous  chapters  need  be  modified  to  apply  to  the 
curved  surface  scatterers. 

The  Idea  Is  to  use  the  concept  of  the  extended  boundary  condition 
[31]  Instead  of  the  ordinary  boundary  condition.  Very  briefly,  the 
extended  boundary  condition  requires  that  the  total  electric  field  be 
vanished  In  the  region  enclosed  by  the  surface  of  the  scatterer.  This 
concept  allows  a portion  of  the  aperture  be  embedded  In  the  obstacle 
as  Illustrated  In  Figure  6.16.  Thus,  the  accuracy  check  of  the  approx- 
imate solution  can  be  done  In  the  Galerkln's  sense,  which  means  that 
(6.44)  Is  employed  as  a criterion  for  validating  the  approximate  solu- 
tion. The  detailed  discussion  of  the  computations  of  the  Inner  products 
In  (6.44)  has  been  presented  in  Subsection  6.4.5  where  Method  V was 
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developed:  hence,  to  save  repetition,  the  analysis  Is  omitted  In  this 
section.  The  only  Important  point  need  be  emphasized  is  that  the  ap- 
proximate scattered  far-fleld  pattern,  T(i^),  whose  accuracy  Is  being 
examined.  Is  to  be  used  to  compute  In  the  Inner  product  of  the 
left-hand-side  member  of  (6.44). 

This  accuracy  checking  procedure  has  been  applied  to  the  approx- 
imate far-fleld  patterns  obtained  by  Method  IV  of  Subsection  6.4.4  for 
the  cases  of  A - 3X  and  A - 6X,  with  the  testing  Gaussian  pulse  located 
at  ■ -15°,  -30°,  and  -45°.  Results  of  these  accuracy  checks  are 
tabulated  In  Tables  6.1  and  6.2. 

Observe  that  the  dominant  parts  of  the  entries  In  the  second  and  the 
third  columns  of  these  tables  are  Identical  to  two  significant  digits. 

The  differences  of  the  entries  In  the  second  and  the  third  columns  are 
shown  In  the  fourth  columns  of  these  tables.  The  maximal  deviation  Is 
about  ±5  percent  from  the  expected  zero  value,  which  shows  that  the 
results  of  the  accuracy  checks  are  quite  satisfactory;  hence,  the  scat- 
tered far-fleld  patterns  of  Method  IV  can  be  concluded  to  be  quite  ac- 
curate without  comparisons  with  the  exact  solutions.  Of  course,  In 
the  case  of  A ■ 3X,  the  exact  scattered  far  field  has  also  been  com- 
puted and  a comparison  with  the  approximate  solution  of  Method  IV  not 
only  further  confirms  the  validity  of  the  approximate  solution,  but 
also  Indicates  that  the  accuracy  check  procedure  Is  effective.  For 
the  scattered  far-fleld  pattern  obtained  by  Method  V,  the  self-checking 
procedure  has  already  been  built  Into  that  method  as  part  of  the  solu- 
tion process;  hence,  a separate  accuracy  check  Is  repetitive. 
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TABLE  6.1 


ACCURACY  CHECK  RESULTS  FOR  THE  A - 3X  CASE 


m’  oAP 

-<P  ,E^> 
m m 

<P  *E®.r,>  + <P  .E^> 
m*  mAP  m m 

-15° 

0.0004793 
+ 11.1452 

- 0.000003465 
+ 11.19533 

0.0004828 
- 10.0501185 

-30° 

1.46487 
- 10.42989 

1.44525 
- 10.48077 

0.0196276 
+ 10.050886 

-45° 

- 2.05293 

- 10.59244 

- 2.016778 

- 10.64182 



- 0.036151 
+ 10.04938 

TABLE  6.2 


ACCURACY  CHECK  RESULTS  FOR  THE  A ■ 6X  CASE 


^ja 

<P  »E 

m mAP 

-<P  ,E^> 
m m 

<P  ,E®.„>  + <P  ,E  > 
m*  mAP  m’  m 

-15° 

1.18760 
+ 10.02917 

1.19917 
- 10.00001 

- 0.01565 
+ 10.029184 

-30° 

- 1.20338 
+ 10.943703 

- 1.21897 
+ 10.911953 

0.015594 
+ 10.03175 

-45° 

- 1.72229 

- 11.24086 

- 1.72114 

- 11.21891 

- 0.0011549 

- 10.0219545 
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6. 7 Sunnnary 

In  this  chapter,  a synthetlc-aperture-dlstrlbutlon  approach  to  the 
high-frequency  electromagnetic  scattering  of  obstacles  with  convexly 
curved  surface  has  been  presented.  The  Idea  was  to  transfer  Information 
characterizing  the  curved  surface  onto  a planar  aperture  where  operations 
for  Improving  the  geomoirlcal-optics  solution  were  carried  out.  The 
Introduction  of  a planar  aperture  in  the  present  approach  Is  a marked 
difference  from  all  of  the  conventional  methods.  A salient  feature  of 
the  method  Is  that  the  far  field  Is  obtained  by  Fourier  transforming 
the  aperture-field  distribution;  hence,  It  can  handle  problems  Involving 
an  n-dlmenslonal  obstacle  by  an  (n  - l)-dlmenslonal  fast  Fourier  trans- 
form (FFT),  where  n ■ 2,3.  Therefore,  the  unwieldy  three-dimensional 
FFT  Is  circumvented:  hence,  the  method  Is  numerically  efficient.  Another 
Important  feature  of  the  method  Is  that  the  aperture-field  distribution 
Is  slowly  varying,  l.e.,  not  rapidly  oscillatory.  In  magnitude  and  In 
phase 4 Hence,  operations  for  improvement  can  be  achieved  In  a relatively 
straightforward  manner.  Moreover,  the  method  Is  conceptually  simpler 
than  other  conventional  methods  In  the  literature,  because  operations 
for  Improvement  are  carried  out  on  a plane  rather  than  on  a curved  sur- 
face. 

The  method  has  been  applied  to  solve  the  scattering  problem  of  a 
plgve  wave  by  a perfectly  conducting  circular  cylinder.  The  underlying 
principles  Involved  have  been  Illustrated  by  a detailed  discussion  on  the 
syntheses  of  the  approximate  aperture- field  distribution. 

The  Galerkln's  method  applied  in  the  spectral  domain  not  only  results 
In  a self-checking  synthesis  of  Method  V,  as  described  In  Subsection  6.4.5, 


182 

but  also  provides  a convenient  accuracy  check  of  the  approximate  solution 
obtained  by  other  methods. 

The  remarkable  agreement  between  the  scattered  £ar**fleld  obtain- 
able by  the  method  with  the  true  solution  makes  It  possible  to  devise 
a means  of  obtaining  the  Induced  surface  current  on  the  obstacle  di- 
rectly from  the  knowledge  of  the  scattered  far  field.  The  computed  sur- 
face current  Is  In  reasonable  agreement  with  the  true  solution. 

Considering  the  efforts  Involved  In  the  formulations  of  the  solu- 
tions and  the  computations  of  numerical  results  In  conventional  methods, 
the  present  approach  Indeed  manifests  Itself  as  a conceptually  simple, 
computationally  efficient  method  for  solving  high-frequency  electro- 
magnetic scattering  problems  of  curved-surface  scatterers. 
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7.  CONCLUDING  REMARKS 

A new  approach  for  combining  the  Integral  equation  and  high  fre- 
quency asymptotic  techniques  has  been  demonstrated  with  four  Illustrative 
examples — diffraction  by  a strip,  a thin  plate,  a rectangular  cylinder, 
and  a circular  cylinder.  The  basic  Idea  Is  to  start  with  the  asymptotic 
high-frequency  solution  leading  to  the  zero-order  approximation  of  the 
scattered  far  field,  and  to  use  the  latter  In  the  Fourier-transformed 
version  of  the  extended  form  of  the  Integral  equation  to  derive  an  Improved 
result  for  the  Induced  surface  current  density  and  the  scattered  far 
field.  By  formulating  the  problem  In  the  spectral  domain,  the  spatial 
domain  Integral  equation  becomes  an  algebraic  equation,  which  can  be 
recast  In  an  Iteration  scheme  suitable  for  manipulations  on  the  computer. 

A salient  feature  of  the  method  Is  that  the  accuracy  of  the  solution  for 
the  surface  current  density  and  the  scattered  far  field  can  be  conve- 
niently checked  by  verifying  whether  the  scattered  field,  which  Is  also 
computed  In  the  process  of  generating  the  solution,  lndp-''.d  satisfies 
the  boundary  condition  at  the  surface  of  the  scatterer.  Therefore,  this 
approach  not  only  provides  a way  for  systematically  Improving  the  GTD 
solution  using  the  self-consistent.  Integral  equation  formulations,  but 
also  provides  a convenient  validity  check  of  the  ray  optics  solution. 
Furthermore,  the  method  of  solution  yields  both  the  Induced  surface 
current  density  and  the  far  field — an  important  feature  which  Is  not 
present  In  conventional  asymptotic  high-frequency  methods. 

The  convergence  of  the  Iterative  scheme  In  Chapter  2 has  not  been 
proven  rigorously  but  has  been  demonstrated  by  numerical  verification 
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only.  Note,  however,  that  the  Iteration  procedure  can  always  be  followed, 
after  a few  iterations,  by  an  application  of  the  Galerkin's  method  in 
the  transform  domain.  The  Iterated  solutions,  being  close  to  the  true 
solution,  provide  a very  convenient  choice  for  the  basis  set  that  may 
Include  some  other  functions  as  well.  In  fact,  for  curved  surface 
structures,  such  as  the  circular  cylinder  in  Chapter  6,  the  Galerkin's 
method  has  been  found  to  be  the  more  desirable  one  to  follow. 

We  would  like  to  mention  two  other  approaches  [14  ] , [4  ] that  are 
based  on  a combination  of  asymptotic  and  Integral  equation  techniques. 

The  one  developed  by  Thiele  [14  ] decomposes  a given  problem  Into  two 
parts,  one  of  which  Is  handled  using  the  GTD  method  and  the  other  using 
the  moment  method.  For  the  case  of  a wire  antenna  on  a finite  ground 
plane,  the  effect  of  the  edge  diffraction  from  the  ground  plane  Is  eval- 
uated using  GTD  and  the  result  Is  subsequently  used  to  augment  the  Im- 
pedance matrix  of  the  roonopole  antenna  over  an  Infinite  ground  plane. 
Although  the  method  works  rather  well  when  GTD  results  are  accurately 
known  for  the  ground  plane  problem,  e.g.,  a ground  plane  of  circular 
shape,  no  convenient  method  Is  available  for  Improving  the  solution  when 
there  are  corners  in  the  plane  that  contribute  substantially  to  the 
scattered  field.  The  latter  situation  arises  when  the  ground  plane  Is 
rectangular  shaped  and  Is  not  large  compared  with  the  wavelength,  or  when 
the  antenna  Is  mounted  close  to  one  or  more  of  the  edges. 

The  second  method  developed  by  Burnside  [4  ] tends  to  rectify  the 
situation  alluded  to  above  by  solving  for  surface  currents  via  the 
moment  method  In  the  regions  where  the  GTD  solution  Is  not  accurate,  and 
by  using  asymptotic  forms  for  the  surface  currents  In  regions  where  a 
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good  approximation  for  these  currents  can  be  enq>loyed.  However,  this 
method  cannot  be  conveniently  applied  to  either  the  strip  problem  with 
grazing  incidence,  or  to  the  large  plate  problem  discussed  in  this  thesis. 
For  the  strip  problem,  the  GTD  solution  is  quite  Inaccurate  when  the  in- 
cident angle  of  the  illuminating  wave  is  near  grazing.  For  the  plate 
problem,  the  current  does  not  settle  down  to  known  asymptotic  form  in 
the  center  region  of  the  plate  until  it  is  at  least  three  to  four  wave- 
lengths squared.  The  moment  method  is  Incapable  of  handling  the  number 
of  unknowns  required  to  accurately  solve  for  the  current  distribution 
on  plate  sizes  that  are  larger  than  2A  squared.  To  demonstrate  that 
the  present  method  is  equally  well  applicable  to  the  same  type  of  practical 
structures  which  have  been  solved  by  the  hybrid  technique  developed  by 
Burnside  [A  ],  the  geometries, including  a rectangular  cylinder,  possess- 
ing multiple  sharp  edges  and  a circular  cylinder  having  a smooth  curved 
surface,  have  been  investigated.  Both  the  relative  accuracy  and  the 
efficiency  of  the  transform  method  have  been  found  to  be  very  favorable. 

The  detailed  description  of  this  work  has  been  presented  in  Chapters 
5 and  6. 

One  other  method  developed  by  Bojarski  [32]  should  be  mentioned 
for  completeness  since  he  also  uses  the  transform  technique  to  convert 
the  Integral  equation  into  an  algebraic  form.  However,  he  uses  a 
three-dimensional  Fourier  transform  which  again  becomes  unwieldy,  both 
in  terms  of  computer  storage  and  time.  In  contrast,  the  present  method 
employs  (n  - 1) -dimensional  transforms  even  for  an  n-dlmenslonal  curved 
surface,  where  n - 2,3,  as  shown  in  Chapter  6,  thus  achieving  a saving 
in  the  storage  requirement  by  approximately  two  orders  of  magnitude 


186 


and  a corresponding  saving  In  CPU  and  I/O  times  as  well.  One  other 
point  with  respect  to  Bojarskl's  work  Is  that  no  advantageous  use  Is 
made  of  the  available  analytic  form  of  the  GTD  solution  for  the  far  field 
In  the  visible  and  Invisible  ranges,  which,  In  many  cases,  forms  an  ex- 
cellent starting  point  of  the  Iteration  procedure. 

In  conclusion,  a conceptually  simple,  computationally  efficient 
method  for  solving  high-frequency  electromagnetic  scattering  problems 
has  been  illustrated.  The  method  has  devised  ways  and  means  to  system- 
atically Improve  the  far  field  as  well  as  the  Induced  surface  current. 

It  has  been  shown  by  examples  that  the  method  Is  applicable  to  scat- 
terers  with  multiple  edges  as  well  as  those  with  smooth  curved  surface. 

The  computational  efficiency  of  the  method  has  been  demonstrated.  Both 
the  far  field  and  the  surface  current  are  obtainable  In  the  process  of 
solution.  Furthermore,  accuracy  check  Is  a built-in  feature  of  the 
method. 

It  was  the  purpose  of  this  study  to  demonstrate  that  a new  approach 
to  the  high-frequency  scattering  is  feasible.  However,  to  develop  the 
method  further,  It  will  be  desirable  to  investigate  Its  applications  to 
more  practical  problems.  The  pioneering  work  done  and  documented  herein 
Is  hoped  to  stimulate  strong  Interests  in  future  research  activities  along 
similar  lines  based  on  the  present  approach.  The  ultimate  goal  of  the  ap- 
proach as  well  as  any  other  approaches  is  to  solve  scattering  problems  In- 
volving arbitrarily  shaped  obstacles.  Therefore,  the  next  stage  of  de- 
velopment of  the  method  should  be  an  extension  of  It  to  cope  with  el- 
liptical cylinders, spheres,  and  ellipsoids.  All  of  these  scatterers 
are  of  fundamental  Interests  In  high-frequency  scattering.  The  fact 
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that  the  present  approach  uses  two-dimensional  FFT  even  for  three- 
dimensional  scatterers  should  make  it  feasible  in  obtaining  reasonable 
results  for  these  scatterers,  which,  in  turn,  can  considerably  augment 
the  scope  of  future  research. 
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APPENDIX  A 

UNDERSTANDING  AND  APPLYING  HIE  FAST  FOURIER  TRANSFORM 

This  appendix  is  written  to  serve  as  a tutorial  presentation  of 
the  Cooley-Tukey  Fast  Fourier  Transform  for  rapid  computation  of  the 
Discrete  Fourier  Transform  of  a function  represented  by  a discrete  set 
of  data  points. 

The  Fourier  transform  of  a function  f(t)  may  be  defined  by 
00 

F(ui)  “ F{f(t)}  ■ / f(t)  e dt  . (A.l) 

•>00 

Both  F(oj)  and  f(t)  are  complex  functions  of  teal  variables. 

Suppose  that  f(t)  Is  not  available  analytically,  but  Instead  Is 
expressed  as  a sequence  of  samples,  f(nT).  If  these  samples  are  taken 
to  be  the  sequence  produced  by  sampling  f(t)  with  a Dirac  delta  func- 
tion sampler,  we  may  write 

00 

f*(t)  - I f(nT)  6(t  - nT)  (A. 2) 

n*-» 

where  the  asterisk  Implied  a sampled  function. 

The  equation  (A. 2)  has  a Fourier  transform 

00 

F*(u)  * f(nT)  e . (A. 3) 

n«-» 

If  f(nT)  - 0 for  n < 0,  (A. 3)  reduces  to  the  familiar  z-transform 
00 

F(z)  - ^ f(nT)  z 
n-0 


(A.4) 


Observation  of  (A. 3)  shows  that  the  Fourier  transform  of  any  train  of 
Impulses  Is  periodic  In  u with  a period  2it/T. 
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A.  1 The  Discrete  Fourier  Transform 

Suppose  that  f(t)  Is  represented  by  a finite  number  of  samples, 
say  N,  such  that  the  samples  are  f(nT),  0 < n < N - 1.  Further,  let 
the  spectral  function  F(u)  be  represented  by  F(kO),  0 < k < N - 1. 

The  variables  T and  0 represent  the  spacing  of  samples  in  time  and  fre- 
quency, respectively.  Then  In  (A. 2)  and  (A. 3)  there  are  only  a finite 
number  of  samples  and  we  can  write 

N-1 

f*(t)  - I f<nT)6(t  - nT)  (A.  5) 

n*0 


and 


F*(kfi) 


■T  f(nl) 

n*0 


(A.  6) 


The  equation  (A. 6)  represents  a discretizing  of  the  Integral  transform 
(A.l)  In  both  time  and  frequency.  Thus,  the  expression  (A. 6)  can  be 
readily  computed.  From  (A. 6)  we  drop  the  asterisk  notation  and  define 
the  Discrete  Fourier  Transform  (DFT)  of  f(t)  as 


F(kD)  - I f(nT) 
n»0 


(A.  7a) 


where 


(A.  7b) 
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We  may  show  that  for  this  choice  of  there  are  only  N distinct  values 
of  F(kn)  which  can  be  computed. 

Consider  an  arbitrary  Integer  k.  This  Integer  can  be  expressed 
as 


k - rN  + kg 


where 


k^  ■ k modulo  N 


and 


r - (k  - kQ)/N 


We  can  thus  write 


N-1  -jfiTn(rN+k.) 

F(kfl)  - I f(nT)  e 
n-0 


f(nl)  . 

n«0 


Now 

OTN  - 2tt  from  (A.  7b),  so  - 1 

r- 

and 


N-1 

F(kfl)  - \ f(nT) 
n«0 


-jnTnkg 

e 
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So  for  k > N - 1,  F(kn)  reduced  Co  F[(k  mod  N)(2].  That  Is  to  say  Chat 
F(kn)  Is  periodic  In  k with  period  N. 

A.  2 Properties  of  the  DFT 

The  DFT  can  be  shown  to  be  frequency  selective.  Suppose  f(nT) 

■ where  we  have  let  u ■ qD  be  a multiple  of  D.  (It  is  not 

necessary  that  q be  an  Integer.)  Then 

F(kD)  - ^ gj(q-k)nfiT 
n«0 


1 . gj  (q-k)NDT 
1 _ J (q-k)DT 


(A. 8b) 


Consider  (A. 8)  for  integral  q.  If  k ■ q mod  N,  gJ  ^9~^)NnT  _ (q-k)DT 
1,  and  (A.  8b)  Is  seen  to  be  a 0/0  form.  However,  it  Is  clearly  seen 
from  (A. 8a)  that  F(kO)  ■ N In  this  case.  For  k q mod  N, the  exponential 
In  the  numerator  of  (A.  8b)  still  has  a value  of  1,  but  the  exponential 
In  the  denominator  does  not,  so  F(kD)  ■ 0. 


F(kfi) 


k ■ q mod  N 
k q mod  N 


(A.9) 


If  q Is  not  an  Integer,  we  may  deduce  from  (A. 8b)  that  F(kD)  ^ 0 for  all 
k.  Thus,  from  (A.9)  we  see  that  the  DFT  Is  frequency  selective  for  fre- 
quency which  Is  an  Integral  multiple  of  0;— a property  which  corresponds 
closely  to  the  frequency  selectivity  of  the  continuous  Fourier  transform. 

The  DFT  Is  a linear  operator  directly  from  the  definition  of  (A.  7). 
This  may  be  expressed  by  the  statement 


DFT{af(nT)  + bg(nT)}  - aDFT{f(nT)}  + bDFT{g(nT)} 


(A. 10) 
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An  Inverse  DFT  exists  which  will  map  a DFT  back  Into  Its  generating 
sequence.  This  Inverse  Is  expressed  by 


fttD-i  I F(kO)  . 

“ k-0 


This  differs  from  (A.  7)  In  the  scale  factor  and  the  sign  of  the  rota- 
tlon  of  the  exponential. 

We  may  show  (A. 11)  to  be  valid  by  substituting  (A. 7)  for  F(kn). 


nm  - i Y Y £(■>« 

k*0  n»0 


iYfCnDY^^”^'”^''  • 

n*0  k»0 


YeJt!T(t-n)k  J” 

k-0  jo  n 1 mod  N 


f (AT) 


f[(A  mod  N)T]N 
N 


f[(A  mod  N)T] 


We  observe  that  the  Inverse  transform  reproduces  the  original  function 
with  a periodic  extension  of  period  N. 
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A.  3 Fast  Fourier  Transform  Algorithms 

The  DFT  of  a function  can  be  calculated  throu{;h  a direct  application 

of  (A.  7).  The  compute  time  required  for  this  calculation  would  be  pro- 
2 

portlonal  to  N If  all  N frequency  points  were  extracted. 

The  Cooley-Tukey  algorithm  uses  the  presence  of  redundant  calcu- 
lations In  (A.  7)  and  of  factorable  multiplications  In  order  to  reduce 
the  compute  time  of  the  computation.  For  N» which  Is  a power  of  two, 

the  computation  Is  proportional  to  N log2  N.  This  provides  substantial 

2 

time  savings.  For  example.  If  N - 1024,  N 1,048,576,  while  N log2  N 

2 

- 10,240,  or  Just  under  one  percent  of  N . 

In  describing  the  Cooley-Tukey  algorithm,  for  notatlonal  simplicity, 
we  define 


(A.  12) 


This  definition  requires  that  we  not  neglect  the  dependence  of  W on  N. 

For  Instance,  when  we  transform  a sequence  with  N/2  points  W would  be 
2 

replaced  by  W . 

Consider  a set  of  N samples  f^,  where  N Is  divisible  by  2.  Let  us 
decompose  the  sequence  Into  a sequence  of  even  numbered  samples  and  a 
sequence  of  odd  numbered  samples: 

I - 0,1,2,  ...,  N/2  - 1 . (A. 13) 

“ ^2£+l 


The  DFT's  of  these  sequences  may  be  written  as 


We  can  also  write  the  DFT  of  the  entire  sequence  as 


n-0 

-“I'' 

1-0 

£-0  £-0  ^ 

-Gj^  + w\  . (A.  15) 

Hence,  the  DFT  of  the  sequence  can  be  written  as  a combination  of  the 

DFT's  of  two  shorter  sequences.  Thus,  through  separating  the  sequences 

the  computation  time  of  (A.  15)  can  be  reduced  to  something  proportional 
2 

to  N + (N/2)  . Further,  since  G and  H have  a period  N/2,  only  half  of  the 
N points  Is  needed  to  reproduce  F. 

If  N/2  Is  also  divisible  by  2,  each  of  G and  H can  be  reduced  to 
the  DFT  of  two  shorter  sequences.  This  procedure,  termed  decimation  In 
time,  can  be  continued  until  the  sequence  Is  reduced  to  subsequences 
of  two  points.  These  subsequences  can  be  DFT'ed,  and  from  these  DFT's 
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the  DFT  of  the  total  sequence  can  be  constructed.  The  computation  time 
In  this  case  can  be  derived  to  be  proportional  to  N log2  N. 

In  the  case  that  N Is  not  divisible  by  two  but  Is  factorable,  the 
above  procedure  may  be  repeated  for  other  divisions  than  two.  The  number 
of  points  involved  In  the  largest  DFT  subsequence  will  be  the  largest 
prime  factor  of  N.  The  decimation  In  "time"  algorithm  may  be  used  In 
exactly  the  same  fashion  for  the  Inverse  transform  by  changing  the  sign 
of  the  exponent  In  the  definition  of  W. 

A similar  algorithm  termed  decimation  In  frequency  has  been  derived 
by  Sande,  Cooley,  Stockham  and  others.  This  algorithm  calculates  even 
and  odd  subscripted  points  In  the  frequency  domain  from  half  sequences 
In  the  time  domain.  The  computational  efficiency  of  this  method  Is 
equivalent  to  that  of  the  decimation  In  time  method. 

This  procedure  can  be  employed  In  an  N-dlmensional  Fourier  trans- 
form. The  algorithms,  being  equivalent,  will  not  be  discussed  here. 

A.  4 Employing  FFT  Subroutines 

There  are  five  FFT  routines  available  through  the  IBM  Common  Computer 
User  Group  Library.  The  routines  are  written  In  USASI  Basic  Fortran 
and  hence  may  be  used  with  most  third  generation  Fortran  processors. 

The  routines  are  designed  such  that  they  take  maximum  advantage  of  data 
form  for  the  problem  at  hand.  Substantial  time  penalties  will  result 
from  using  a program  general  beyond  the  requirements  of  a given  problem. 
Table  A.l  lists  the  routines  and  their  functional  characteristics. 

The  user  Is  referred  to  IBM  documentation  for  the  details  of  using  these 


routines . 
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TABLE  A.l 
FFT  SUBROUTINES 


FOURI 


FOURG 


FOURZ 


FOURT 


FOR2D 


One  dimensional 
Power  of  two  data  points 

General  one  dimensional 

Multidimensional 

Power  of  two  data  points  in 
all  dimensions 

General  multidimensional 

Multidimensional 
Power  of  two  data  points 


Data  stored  on  direct  access  file 


